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| THERMAL-EXPANSION CHARACTERISTICS OF SOME 
GROUND-COAT ENAMEL FRITS 


By William N. Harrison, Benjamin J. Sweo, and Stephen M. Shelton 


ABSTRACT 


Thermal-expansion determinations were made on 11 sheet-iron ground-coat 
frits of the single-frit type, using the interferometer and also a quartz-tube dial- 
indicator apparatus. “Critical”? temperatures for both heating and cooling were 
determined. Mean coefficients of expansion from 25° to 400° C, and total contrac- 
tion from the critical temperatures to 25° C during cooling were correlated with 
variations in composition. Total differential contractions between enameling- 
type iron and the frits were also computed and correlated with variations in com- 
position. Constants were determined which, when substituted in simple first- 
order equations, permit the computation of the expansivity of each frit at any 
temperature or temperature interval between 25° C and the critical temperatures. 

The observed difference in the critical temperature on heating and that on 
cooling, the difference in the temperature-expansion curve of highly strained 
specimens below the critical temperature in heating, and the shrinkage of annealed 
specimens when held 30° C below the critical temperature for heating, are ex- 
plained by the work of Tool, Lloyd, and Merritt on glass. 
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I. INTRODUCTION 


A series of 11 ground coats for sheet iron of the single-frit type 
were studied at this Bureau, and reports on their fusion properties,’ 
strength, and elasticity ? have been published. Subsequently, deter- 
iinations on their thermal expansions were completed and consider- 
able work done in comparing the results obtained by different methods 
of test and of analyzing the data. Further, some experiments were 


Say ae eccoee 

‘W.N. Harrison and B. J. Sweo, Some fusion properties of ground coat enamels as influenced by composition. 
BS J. Research 10, 189 (1933) RP524. 

W. N. Harrison, 8. M. Shelton, and W. H. Wadleigh, Strength and Young’s modulus of some ground 
Coat enamels for sheet iron. J. Am. Ceram. Soc. 18, 100 (1935). 
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conducted to obtain confirmation of the logical presumption that the 


temperature-time-volume relations, which, 


according to the work of 


Tool, Lloyd, and Merritt * apply to glasses studied by them, als, 


apply to the enamel frits used in this investigation, which were 
a 


essen. 


ly homogeneous glasses. This paper is a report of these studies 
The method of preparation of the frits, and also their compositions 
both as computed from batch weights and as determined by chemicaj 
analysis, have been published.* The scheme of variation in compos 
tions is shown in table 1. 


TABLE 1.—Scheme of variation in composition of frits and data on critical lempera- 


tures and thermal-expansion characteristics 
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! Each frit contained, in addition to the tabulated constituents, 7.7 percent of fluorspar, 1.7 percent of 


manganese oxide, and 0.6 percent of cobalt oxide. ; 
1 The values for cooling may be obtained by subtracting 60° C from those determined for heating, an ap- 


proximation which is explained in the text. 


II. METHODS OF TESTING 
1. DIAL METHOD 


The apparatus used was similar to that described by Hidnert and 
Sweeney,° with the following refinements. A sensitive dial indicator, 
which could be read to 0.0001 inch and estimated to an additional 
decimal, was used. To minimize temperature gradients, the furnace 


winding was tapped to divide it into three sections. The middle sec- 
tion was the same length as the specimen and inclosed it during a test. 
The two end sections were each about half the length of the specimen. 
Current could be shunted through variable resistances in parallel 
with the respective sections of the winding. Correction for the ex- 


3 Dimensional changes caused in glass by heating cycles. BS J. Research 5, 632 (1930) RP219. Also J. Am. 


Ceram. Soc. 13, 4 (Sept. 1930). ; i 
4 W.N. Harrison and B. J. Sweo, Some fusion properties of ground coat enamels as influenced by composition 

BS J. Research 10, 189 (1933) RP524. ‘ ’ Re 
§ Peter Hidnert and W. T. Sweeney, Thermal expansion of magnesium and some of its alloys. BS J. 


search 1, 771 (1928) RP29. 
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ansivity of the particular fused quartz used, and for temperature 
cradients in the quartz tubes during test, was made by calibration 
against a specimen of porcelain, the expansivity of which had been 
accurately determined at this Bureau. It was established that the 
expansivity - this specimen did not change in repeated determinations 
to 400° C. 
EThe dial was read when the furnace was at room temperature and 
at 100°, 200°, 300°, and 400° C, respectively. At each temperature 
the dial reading used in computing expansion had remained constant 
over a 20-minute period, during which the temperature of the furnace 
was held constant, the resistance being adjusted so that the tempera- 
ture was uniform within 1° C from top to bottom of the specimen. 
The requisite 8-inch specimens were cast in molds made from blocks 
of electrode carbon. Side-opening molds were used so that the vessel 
from which the molten frit was poured could be moved from one end 
of the mold to the other during the pouring process. Each mold was 
preheated to about 600° C and, after the specimen was cast, was re- 
tuned to the furnace, held at 550° C for 1 hour and allowed to cool 
with the furnace. While in the furnace the molds were encased in 
sheet-iron jackets buried in mixtures of sand and graphite to retard 
oxidation. Expansion curves on specimens prepared in this way 
indicated a satisfactory quality of annealing. 


2. INTERFEROMETER METHOD 


The construction and operation of the apparatus used for measuring 
expansion by the interferometer method are described by Merritt.® 
The test specimens were made from rods approximately 3 mm in 
diameter drawn from the molten frit. A blunt point was formed at 
one end of each specimen, and three short legs formed on the opposite 
end, by grinding. These specimens were approximately 5 mm long, 
over-all. The temperature was increased continuously in these tests 
it the rate of 3° C per minute, as indicated by the average of two 
Chromel-Alumel thermocouples, the junction of one being placed just 
above the cover plate, and the junction of the other just below the 
vase plate. The cup in which the specimens were placed during test 
was made of aluminum. Tests were made on unannealed specimens 
cooled from the molten state in air to determine their “softening tem- 
peratures.”’ 7 They were then annealed by heating for 3 hours at a 
temperature about 5° C below their respective softening temperatures, 
and cooling in the furnace. These annealed specimens were used for 
expansivity determinations. 


III. RESULTS AND DISCUSSION 
1. EXPANSIVITY TO 400° C 


The coefficients of linear thermal expansion from 25° to 400° C, as 
determined by both the dial and the interferometer methods, are given 
in table 1. The initial temperature in each case was approximately 
23° C and the values were all corrected to 25° C in order that they 


mgt be comparable. The coefficient of expansion values given in 
saecinnnan 
au. Merritt, The interference method of measuring thermal expansion. BS J. Research 10, 59 (1933) 


‘W.N, Harrison and B. J. Sweo, Some fusion properties of ground coat enamels as influenced by composition. 
DS J, Research 10, 189 (1933) RP524. - _— os . ” sie 
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this table represent averages of determinations on three separately 
prepared batches of frit. Individual deviations from the report; 
values averaged less than 0.1X10~°, the maximum deviation being 
0.3 10-°. 4 

The values obtained with the interferometer are 0.2 to 0.5% ~ 
higher than those obtained with the dial apparatus. A number of 
experiments made to determine the source of this difference failed to 
reveal the cause of the discrepancy. 

From table 1 it may be seen that the assumption commonly used jy 
calculating expansivity from factors, namely, that the expansivity js 
an additive function of the constituent parts, held approximately 
true for these frits. Each of the six substitutions of 3 percent of boric 
oxide for 3 percent of sodium oxide, with other constituents unchanged. 
resulted in a decrease of about 1.0 10~° in the coefficient of expansion, 
the extremes being 0.8 and 1.110-°. The substitution of 10 percent 
of flint for 10 percent of feldspar produced a decrease of about 
0.5X10-°8 The J-K change in composition resulted in a substan 
tial reduction in the coefficient of expansion. 


2. “CRITICAL” TEMPERATURES 


It is recognized that the real importance of the thermal expansivyity 
of enamels lies in the differential contraction between the enamel and 
the underlying metal as an enameled article cools after firing, and this 
Dioemdial damande not only upon the difference in expansivity but 


also upon the temperature at which the enamel becomes sufficiently 
stiff to set up strains in cooling. Hence it is obvious that the coefl- 


cient of expansion alone, whether computed or determined experi- 
mentally, does not give all of the required information. 

The temperatures at which these strains were initiated were not 
determined directly. However, the ‘critical’? temperatures °® of the 
frits were determined because it was thought that a knowledge of the 
total contraction from these temperatures to room temperature would 
be helpful in estimating the relative amounts of residual strain after 
firing. In locating these critical temperatures, only the interferom- 
eter method was used. No attempt was made to carry the observa- 
tions with the dial apparatus up to such temperatures, as the use of 
the large specimens, and the accompanying necessity for holding the 
temperature constant at each observation level, render this method 
unsuitable for the purpose. Data obtained during cocling were com- 
pared with those obtained during heating and, for greater accuracy, 
new determinations were made, using a more sensitive potentiometer, 
the precision of the new pyrometer being better than 0.01°C. The 
cooling rate of the furnace was approximately 2°C per minute through 
that temperature range which included the critical temperatures, 
becoming slower as the temperature decreased. 

Figure 1 is indicative of the results that were obtained. In the 
upper curves, expansion and contraction with change of temperature 
are plotted in the usual manner. In the lower curves the rates (00- 
efficients) of expansion (or contraction) for the small temperature 
intervals between pairs of adjacent interference bands are plotted 
against the corresponding mean temperatures. Since the lower curves 

* This 10-percent change in batch composition entailed a much smaller change in chemical composition. 


°C. G. Peters and C. H. Cragoe, Measurements on the thermal dilatation of glass at high temperatures, BS 
Sci. Pap. 16, 445 (1920) S393. 
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are the first derivatives of the upper ones, they accentuate changes in 
rate of expansion or contraction. Straight lines fitted the points on 
the derivative curves sufficiently well up to approximately the critical 
temperature, which may be defined for the purpose of this work as the 
temperature at which a smooth curve drawn through the plotted 
points diverges or “breaks away” from the straight line.” That this 
“break”’ was sharper and that it occurred at a substantially bigher 
temperature on heating than on cooling, as shown in figure 1, held 
true for every frit tested. 

Even with the breaks accentuated by this system of plotting the 
data, it was not possible to determine the exact temperature at which 
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FicurE 1.—Curves showing the expansion and contraction of frit B. 


The upper curves show the expansion from 350° to 500° C and contraction from 500° to 50°C. Thelower 

rves were obtained from the same data and show the change in coefficient of expansion (or contraction) 
with temperature under the conditions of test. To avoid confusion, not all of the data are plotted, except 
at the high-temperature end of the lower curves. 


each break occurred, but it is thought that the uncertainty of the 
values obtained from the heating curves is within +5° C. However, 
itis believed that the differences between the critical temperatures of 
different enamels were obtained with greater accuracy, since all were 
judged in the same way. The critical temperatures given in table 1 
ae those for heating data. Because of the greater difficulty in 
loceting the critical temperatures for cooling from the data as plotted 
in figure 1, an estimate was made of the average difference, and this 
value (60° C) was subtracted from the critical temperature obtained 
fom a heating curve to obtain the critical temperature for the 
corresponding cooling curve. 
_ One correlation between change in composition and resultant change 
in critical temperature is very striking. The A-B, D-E, and G- 
Although the points of divergence will not be accentuated to the same degree as in the described pro- 
ure, the same principle can be used by fitting a curve to expansion (or contraction) data from the con- 


os referred to further on in this paper, and noting the temperatures at which smooth lines drawn through 
# data points diverge from such curves. 
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changes in critical temperatures (see table 1) are about four times os 
large as the B-C, E-F, and H-I changes, although the changes }, 
composition are in the same direction and are of the same magnitude. 
The first three pairs all involve an increase of boric oxide from 13 to 
16 percent and a reduction of sodium oxide from 17 to 14 percent 
The second three pairs involve a further increase in the boric oxide 
to 19 percent, with corresponding decrease in sodium oxide to }} 
percent. The average effect of increasing flint from 25 to 35 percent 
at the same time reducing feldspar from 35 to 25 percent (the 4-( 
B-H,and C-I changes), was to increase the critical temperature on|y 
about 7° CC." The J-K change in composition resulted in an increas 
of 38° C in critical temperature. ‘ 


3. EXPANSION CONSTANTS 


The following equation has frequently been used to express changes 
of length with temperature, as graphically illustrated in the upper 


curves of figure 1. 
L/In=1+at+b0? (1) 


in which: 

L=length of specimen at the temperature, t. 

Iyn=length of specimen at the temperature, 0° C. (No 
appreciable error is introduced if this length js 
measured at room temperature.) 

a and b are constants. Bx 

The instantaneous coefficient of expansion, C,, at any temperature, 
t, is therefore given by the derivative of eq 1. 


=. 
C= Hh Xz = e+ 2bt, (2) 


from which it follows that the mean coefficient between two tempers- 
tures, ft, and fg, is: 


a+6(t,+t) 


Onqa at". 


The constants, a and 6, for the enamels studied are given in table |. 
Although these constants are empirical, and small variations in the 
data may make relatively large percentage differences in their values, 
they are, nevertheless, useful in several ways. Various investigators 
have reported thermal expansions of enamels and similar materials 
between different limits of temperature, and hence it is sometimes 
difficult to compare the results reported by different authors. The 
constants given in table 1, when used in eq 3, make it possible to 
compute the coefficient of expansion of the enamels within any 
desired temperature range up to the critical temperature of a frit by 
a very simple calculation. Furthermore, unless elaborate precautions 
are taken, it is difficult to obtain accurate data on thermal expansion 
near room temperature when a schedule of continuous increase 
temperature is used."2 This difficulty may be largely remedied by 
use of the constants, as indicated in the following discussion. 


11 See footnote 8, p. 130. (1933) 
11G. E. Merritt, The interference method of measuring thermal erpansion. BS J. Research 10, 59 (100s 


RP515. 
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In this study determinations of total expansion from 25° to 400° C 
and hence mean coefficient of expansion) were in satisfactory agree- 
ment whether obtained in heating or in cooling. However, the 
qreement was not satisfactory at the lower temperatures, notably 
irom 50° to 150° C. Below 50° C no readings were taken in cooling. 
Since temperature measurements were made from two thermocouples, 
one above and one below the quartz plates, and a quartz thermometer 
was not used, it was concluded that the temperature readings may 
pave differed significantly from the true temperature of the specimen 
juring heating. In the cooling tests the average rate of change of 
temperature in the 50° to 150° C interval was about 0.3° C per 
minute, as compared with a rate of 3° C per minute in heating, and 
because of this slower rate, the temperature gradients were presumably 
snaller. When a very slow heating schedule was tried, a higher 
apparent rate of expansion was observed, as in cooling. Hence, in 
computing the constants given in table 1, the data obtained in cooling 
were used in establishing the lower end of the expansion curves down 
» 50°C. Extrapolations from 50° C down to 25° C by means of the 
constants used in eq 3, are thought to be more reliable than the 
observed data in that range. 

In table 1, in addition to the constants for eq 3, the mean linear 
coeflicients of expansion computed from these constants are given 
for the following temperature ranges: 25° to 200° C, 25° to 400° C, 
25° C to the eritical temperature, as determined from heating data, 
and 25° C to the critical temperature, as determined from cooling data. 

Concerning the total expansions to the critical temperatures in 
heating, or total contraction from critical temperatures in cooling, a 
reversal of trend takes place as sodium oxide is substituted for boric 
oxide. Increasing the sodium oxide content from 11 to 14 percent 
and decreasing boric oxide from 19 to 16 percent (the C-B, F-E, and 
I-H changes shown in table 1) caused a substantial increase in total 
expansion to (or contraction from) the critical temperatures, as well 
asin the average coefficient of expansion. A further 3-percent increase 
of sodium oxide at the expense of boric oxide (the B-A, H-D, and 
H-G changes shown in table 1) increased the coefficient of expansion, 
as in the first substitution, but the accompanying reduction in critical 
temperatures more than compensated the effect of this increase in 
average coefficient of expansion on the total length change between 
25°C and the critical temperatures." 

It would be a mistake, however, to deduce that the total differential 
contraction between iron and enamel in the temperature ranges under 
discussion undergoes a reversal, as boric oxide and sodium oxide are 
varied between the indicated limits. From values given in table 1 
it is apparent that the differential contraction increases without 
ieversal as the coefficient of expansion of the frits decreases. The 
total contraction of the iron to 25° C from the respective critical tem- 
peratures of the frits in cooling, given in table 1, was computed from 
constants which were experimentally established for one specimen of 
ehameling-type iron in the range 25° to 400° C, namely, a=11.4 107° 
and b=0.0055 X 1078. 


_' A nonlinear effect on the properties of glasses caused by intersubstitution of boric oxide for silica was 
doted by W. E. S. Turner, J. Glass Tech. 7, 155 (1923). 
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4. APPLICATION OF CONCLUSIONS OF TOOL, LLOYD, AND MERRITT? 
TO FRITS STUDIED 


[Vol, 9 


The described shift in the “critical” temperature obtained on cool. 
ing, as compared with that obtained on heating (fig. 1), has beey 
explained by Tool, Lloyd, and Merritt." The critical temperature 
observed on heating indicates a behavior analogous to that occurring 
in superheating and the critical temperature observed on cooling 
bears a similar analogy to under cooling. In either case the high rates 
of change of volume with respect to temperature illustrated by the 
lower curves in figure 1, between 400° and 500° C indicate drifts 
toward the equilibrium volumes, for the temperatures involved, which 
proceed at insignificant rates in lower temperature ranges, because of 
the higher viscosity of the specimens. The familiar dip in the ey. 
pansion curve of highly strained specimens of enamel frit (or glass) 
which begins as much as 100° C below the critical temperature in 
heating, is also explained by those authors in the same way. Figure? 
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Figure 2.—Thermal expansion curves for frit H. 


Drawn from molten frit and cooled in air (solid circles) and the same specimens after annealing (open 
circles). 

















se 





shows typical heating curves for specimens of frit H in the annealed 
and the unannealed conditions. The relatively flat portion of the 
latter curve between 400° and 500° C reflects a drift toward the 
equilibrium volume-temperature line. The forces tending toward 
equilibrium volume are larger in the highly strained specimens and 
hence can operate more rapidly against high viscosities. 


14 Dimension changes caused in glass by heating cycles. BS J. Research &, 627 (1930) RP219. Also J, Am. 
Ceram. Soc. 18, 4 (Sept. 1930). 
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The fact that even annealed specimens will undergo such drifts if 
riven sufficient tame was shown by Tool, Lloyd, and Merritt. In the 
present study, annealed specimens of enamel frit E were held for 5 
hours at 450° C, which is 30° C below the critical temperature obtained 
on heating and 30° C above that on cooling. The drift toward equi- 
ibrium volume at that temperature resulted in a decrease in length 
equal to 13 percent of the expansion from 25° to 450° C (see fig. 3). 
Upon subsequent reheating no significant change in coefficient of 
expansion was detected. 
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Figure 3.—Curve showing the expansion of annealed specimens of frit E heated at 
8° C per minute to 450° C and subsequent contraction on holding at that temper- 


ature for 6 hours. 
IV. CONCLUSIONS 


Aconsistent difference of about 3 percent was observed between the 
expansivities from 25° to 400° C, as determined with the interferome- 
ter, using specimens 5 mm tall and 3 mm in diameter, and those 
determined with a quartz-tube and indicating-dial apparatus, using 
specimens 8 inches long and %inch square in cross section. Since inter- 
ferometer specimens cut from the larger specimens gave the same 
difference in results, it was concluded that this difference was inherent 
in the methods of test and may have been related to the difference in 
sie of the specimens. Because the smaller specimens are better 
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adapted for determining change of length with continuous change j; 
temperature, especially through the “critical” range, the interferon. 
eter method is preferred. F 

Equal substitutions of boric oxide for sodium oxide in the rance 
studied caused approximately equal decreases in coefficient of expan. 
sion, but caused unequal increases in “critical’”’ temperature, the first 
3-percent substitution causing about four times as much change in this 
property as a further 3-percent substitution in the same directioy 
The total expansion to critical temperature underwent an actya| 
reversal in the same range of substitutions, that for the intermediate 
composition being greater than that for either higher or lower sodiyr 
oxide content. Computation from data of effective differential cop. 
tractions between the frits and a typical enameling-type iron indicated, 
however, that this important relation followed a continuous trend 
increasing with decreasing sodium oxide content. 

From expansion constants and a simple first-order equation whicl 
are given, the coefficients of expansion of the frits studied can be easily 
computed for any desired temperature or any temperature intery,| 
between 25° C, and the critical temperatures. 

The report of Tool, Lloyd, and Merritt on the dimensional changes 
resulting from heating cycles in glasses explains why the data obtained 
in the present study indicate that the critical temperature for cooling 
was lower than that for heating and also explains the shrinkage of 
annealed specimens of enamel frit held for 5 hours at a temperature 
30° C below the critical temperature for heating. The difference 
between the critical temperatures was about 60° C and the shrinkage 
was about 6 u/em. Both of these phenomena, as well as the familiar 
dip in the expansion curve of a highly strained specimen, which occurs 
well below the critical temperature observed in heating, apparently 
result from drifts toward the equilibrium temperature-volume condi- 
tions which those authors reported for glasses. 


WasuinetTon, October 16, 1938. 
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SEPARATION OF 1,2,4-TRIMETHYLCYCLOHEXANE AND 
AN ISONONANE FROM A MIDCONTINENT PETRO- 


LEUM ? 
By Joseph D. White,? and Augustus R. Glasgow, Jr. 


ABSTRACT 


A fraction of Oklahoma petroleum boiling at 141° C, from which the aromatic 
hydrocarbons previously had been removed, was separated by distillation at 215 
mm Hg into a fraction containing the bulk of a naphthene constituent and one 
enriched in paraffins. From the former, nearly pure 1,2,4-trimethyleyclohexane 
was isolated by erystallization from solution in liquid dichlorodifluoromethane. 
(Continued distillation of the paraffinic fraction at normal pressure, alternated 
with distillation at 215 mm Hg, yielded a fraction containing 85 mole percent of 
an isononane, probably 2,3-dimethylheptane. The 1,2,4-trimethyleyclohexane 
constitutes about 0.1 percent, and the isononane about 0.05 percent of the original 
petroleum. 

Effective separation of a paraffin-naphthene mixture, constant-boiling at normal 
pressure, by distillation at a different pressure arises from the greater change in 

boiling point with pressure of the naphthene component. For hydrocarbons 

normally boiling near 140° C, the mean interval between boiling points at 760 mm 
Hgand 215mm Hg is about 42° C for paraffins, and about 43.2° C for naphthenes. 
The boiling points at the two pressures were determined for six hydrocarbons, 
2.6-dimethylheptane, a nonanaphthene (boiling at 136.7° C), 1,2,4-trimethyl- 
cyclohexane, 4~ methyloctane, 3-methyloctane, and 2- -methyloctane. The boiling 
points at these pressures for the ‘“2,3-dimethylheptane” were estimated from 
values for an impure sample. 

The boiling point, freezing point, density, refractive index, and critical solution 
temperature in aniline have been determined for the 1,2,4- trimethyleyclohexane 
and the isononane. Comparison of the properties of the 1,2,4-trimethylcyclo- 
hexane with those reported for its synthetic stereoisomers indicates that the pe- 
troleum hydrocarbon is either leis, 2¢raney4cie-trimethylcyclohexane or l¢ie2traney 
4anetrimethyleyclohexane. 
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I. INTRODUCTION 


In 1883 Markownikoff and Ogloblin [1]* identified 1,2,4-trimethy). 
cyclohexane (hexahydropseudocumene) in Caucasian petroleum from 
the product of its nitration. Since then, frequent allusions have been 
made to the presence of nonanaphthenes (C,H) in petroleum without 
specific reference to 1,2,4-trimethyleyclohexane. Recently, however. 
Kasansky and Eliséeva [2] dehydrogenated the nonanaphthene frae- 
tion, distilling between 135° and 138° C, of Baku petroleum into 
aromatic material, and, by brominating the portion of the product 
which boiled at 170° C, obtained some _ tribromopseudocumene. 
From this experiment, they concluded that 1,2,4-trimethylcyclo. 
hexane was a constituent of the original petroleum fraction. — 

The present investigation has resulted in the isolation of 1,24- 
trimethylcyclohexane and the partial separation of a nonane, proba- 
bly 2,3-dimethylheptane, from the methyloctane fraction [3] of an 
Oklahoma petroleum. Until now, the presence in petroleum of the 
isononane (boiling point, 140.7° C) appears to have escaped detection, 


II. PREVIOUS FRACTIONATION OF THE PETROLEUM 
DISTILLATE 


After a systematic distillation of the whole naphtha fraction of the 
petroleum, the material boiling between 130° and 145° C was treated 
with liquid sulfur dioxide to extract the bulk of the aromatic hydro- 
carbons [4]. To complete the removal of aromatic material, the 
immiscible portion was redistilled, and the large-volume fraction 
boiling between 140° and 145° C was treated with silica gel [3]. 
The residual oil consisting of paraffins and naphthenes was then 
distilled at a pressure of 215 mm Hg through a fractionating column, 
packed with jeweler’s locket chain, having the equivalent of about 60 
theoretical plates [5]. The distilled material was recovered in 50-ml 
fractions. After the measurement of their boiling points and refrac- 
tive indices, fractions boiling within a range of 0.2° C were combined. 
The distillation separated the material into two main portions, each 
with a boiling range (at 215 mm Hg) of about 1° C. The larger, 
boiling between 101° C and 102.2° C (143° and 144.2° C at 760 mm 
Hg), contained the three methyloctanes. The isolation of these 
hydrocarbons from the petroleum has already been described [3]. 
The other portion distilled chiefly between 98° and 99° C. The 
volumes and refractive indices of the combined fractions comprising 
this portion are given in figure 1. 


4 Figures in brackets indicate the literature references at the end of this paper. 





nie | 1,2,4-Trimethyleyclohezane and Isononane from Petroleum 139 


Glasgow 


ITERS 
@ 
° 
° 


' 
ee 


= 
- 
= 


IN 





VOLUME 





960 966 976 98.4 ~ 100.0 100.8 
BOILING RANGE OF FRACTION IN °C (215 mm) 

















Figure 1.—Volumes and refractive indices of the petroleum fractions, boiling between 
96.0° and 100.8° C at 215 mm Hg, after removal of aromatic hydrocarbons and 
subsequent distillation. 


Lower graph, volume distribution with respect to boiling range; upper graph, refractive indices. 


II], COMPOSITION AND PROPERTIES OF THE FRACTION 
BOILING BETWEEN 98° AND 99° C AT 215 mm Hg 


In figure 2 the boiling points and refractive indices of the material 
represented in figure 1 are plotted against the volume distilled. The 
refractive indices are those for the 50-ml fractions which marked a 
rise of 0.2° C in the boiling point of the distillate. The plateau in the 
boiling-point curve between 98° and 99° C represents nearly 3 liters 
of material whose gradual change in properties indicated it to be a 
mixture. The values of the boiling points and refractive indices signi- 





140 Journal of Research of the National Bureau of Standards {yx » 


fied that the material was composed of nonanaphthenes and igono. 
nanes. As disclosed by the high refractive index of the lower-boijlin: 
material, the naphthenes were more volatile than the paraffins at ; 
pressure of 215 mm Hg. 

To determine the degree of separation of the material into its con. 
ponents when distilled at normal pressure, a 1,200-ml portion boiling 
between 98.4° and 98.8° C at 215 mm Hg was distilled at atmospheric 
pressure in the column used before. The results of the distillatioy 
are given in figure 3. The boiling-point curve (labeled 760 mm) 
shows that fully half of the material distilled between 141.0° and 
141.1° C. Somewhat more surprising was the fact, as indicated by the 
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Figure 2.—Fractionation by distillation of the material represented by the graphs in 
jigure 1. 


Boiling points and refractive indices are plotted against the volume of material distilled. The individua 
points are those for 50-m] fractions which marked a rise of 0.2° C in the boiling point of the material. 


refractive-index curve drawn through the open circles, that at 760 mm 
Hg the order of distillation of the constituents was the reverse of 
that at 215 mm Hg. At normal pressure the paraffinic material was 
more volatile than the naphthene and tended to concentrate in the 
lower-boiling fractions.’ It was hoped that the normal boiling points 
of the component hydrocarbons could be located by the atmospheric 
distillation and compared with those of known hydrocarbons to 
obtain a clue as to what compounds might be present. Little knowl- 
edge was obtained concerning the boiling points of the constituents. 
What fractionation there was indicated that the naphthene present 
boiled above 141.0° C and the paraffin below. 


5 Prior to this experiment, S. T. Schicktanz had observed similar behavior of material from the ‘‘isodecane” 
fraction of this petroleum. 
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Since it was difficult to fractionate the material at normal pressure, 
the fractions were combined and distilled again at 215 mm Hg. The 
change in composition of the material as the distillation progressed 
isindicated by the curves passing through the solid circles in figure 3. 
The continual rise in boiling point and the marked change in refractive 
index indicated that better separation took place at the reduced 
ressure. 
i combustion analysis the composition of a fraction of the dis- 
tillate boiling at 98.45° C (n%=1.4151) was 48 mole percent of naph- 
thenes and 52 mole percent of isononanes. The cooling curve of a 
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Ficure 3.—Fractionation of the paraffin-naphthene mixture by distillation at nor- 
mal pressure contrasted to that obtained by distillation at 215 mm Hg. 


Curves passing through the open circles represent the change in boiling point and refractive index of the 
{actions obtained from the distillation at 760 mm Hg. Curves passing through the solid circles denote the 
fractionation at 215 mm He. 


ilml fraction boiling at 98.45° C, (n?3=1.4154) is shown in figure 4. 
By repeatedly plunging a metal rod, chilled with liquid air, into the 
material as it was cooled, freezing was finally induced at —120° C. 
Qn warming the partially frozen mass, a rapid change from a higher 
toa lower viscosity took place around —115° C, and at —104° C the 
lsst of the crystals melted. The cooling and warming behavior indi- 
tated that the material was substantially a binary mixture. As de- 
11827339 2 
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termined later, the portion of the melting curve, figure 4, extendine 
from —115° C to —104° C represents the melting of the naphthene 
component, while the break in the curve at —115° C is attributed 
to a change in state of the isomeric nonane present (see section VII) 
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Figure 4.—Cooling and warming curves of a distillation fraction boiling at 98./° 
C at 215 mm He. 


S, stirred with a cold rod to induce crystallization. 


IV. ISOLATION AND IDENTIFICATION OF 1,2,4-TRI- 
METHYLCYCLOHEXANE 


When the material boiling at 98.4° C was mixed with liquid dichlo- 
rodifluoromethane, Freon 12, and cooled, a crystal fraction was 
produced which was found to be enriched in a naphthene having 4 
normal boiling point of approximately 141.3° C. Systematic crystal- 
lization from this solvent, alternated with distillation, yielded from 
the material boiling between 97.0° and 98.8° C at 215 mm Hg some- 
what more than a liter of the impure naphthene.® This was subjected 
to further fractional crystallization to yield a 95-ml fraction of “best” 


¢ The lowest-boiling mother liquor contained another naphthene, probably 1,3,5-trimethylcyclohexane. 
The higher-boiling mother liquors were rich in paraffins. 
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naphthene. After distillation through a small fractionating column 
and filtration through silica gel to remove remaining traces of solvent, 
this fraction had the following properties: ” 

Normal boiling point, °C 

Pressure coefficient, °C/mm Hg at 760 mm 

Melting point, °C 

Density at 20° C, g/ml 

Refractive index, n? 

Critical solution temperature in aniline, °C 

The freezing range of 1° C in the cooling curve, figure 5, showed 

‘hat the compound was sufficiently pure for use in identification. 
(ombustion analysis indicated that its chief impurity was an isomeric 
nonane. As a result of three determinations by the procedure de- 
«ribed by Rossini [6], the mean ratio, moles H,O/moles CO,, was 
found to be 1.0060 +0.0004. Compared with the theoretical value 
of 1.0000 for a naphthene, C,H2,, and 1.1111 for a nonane, CyH., the 
combustion analysis indicated that the material was composed of 
46 mole percent of a naphthene and 5.4 mole percent of an iso- 
nonane. 
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ficure 5.—Time-temperature cooling and warming curves of the 1,2,4-trimethyl- 
cyclohexane isolated from petroleum. 


The molecular weight of the fraction, determined by the lowering 
freezing point of cyclohexane,’ was 128.5, which agrees within the 
mor of the method with the theoretical value of 126.2 for a nona- 
naphthene, CyHyjs. 
The nonanaphthene was identified as 1,2,4-trimethylcyclohexane 
y dehydrogenation over palladinized asbestos (20 percent of palla- 
lum). The catalyst was prepared by the method described by 
Husky and Borisoff [7]. A 5-ml portion of the petroleum hydro- 
whon, having n?*= 1.4229, was vaporized at the rate of six drops a 
lunute and passed through a tube containing a 50-cm section of the 
italyst held at a temperature of 340° to 380° C. The passage 
‘quired 25 minutes, during which time hydrogen was constantly 
volved. A hydrocarbon product of 3.2 ml was obtained which had 
telractive index at 25° C of 1.4827, and smelled strongly of pseudo- 
nene (1,2,4-trimethylbenzene). This aromatic material was sulfon- 
"The boiling point was measured in a Swietoslawski ebulliometer. The difference between the boiling 


# condensing points was found to be 0.013° C. The density was determined by the Division of Weights 
Measures of this Bureau. 


“ Yalue of 7.4 cal/g was used for the heat of fusion of cyclohexane. See Parks and Huffman, Ind. Eng. 
mM. %, 1138 (1931). 





144 Journal of Research of the National Bureau of Standards {vu 


ated at 45° C with a mixture of 20 ml of concentrated sulfuric aciq 
and 1 ml of fuming sulfuric acid (30 percent of free SO.) and they 
hydrolyzed by steam distillation. After rejecting the oil distilling 
near 100° C as unsulfonated material, 1.3 ml of hydrocarbon wa; 
produced by hydrolysis of the sulfonated material at 145° C. The 
product had a refractive index at 25° C of 1.5022, as compared wit) 
an observed index of 1.5017 for a reference sample of pseudocumerne 
[14]. Confirmatory evidence of the identity of the dehydrogenated 
naphthene with pseudocumene was obtained by comparing the 
melting point of the former’s nitro derivative with that of pure 
trinitropseudocumene, and by determining the melting point of g 
mixture of the two nitro compounds. Comparative values for the 
three melting points were as follows: 

(1) From the dehydrogenated naphthene 

(2) From pseudocumene 

(3) For a mixture of (1) and (2) 

The results of the dehydrogenation and nitration experiments 
leave little doubt that the petroleum naphthene is hexahydropseudo- 
cumene, namely, 1,2,4-trimethylcyclohexane. However, further in. 
formation regarding the naphthene’s identity is obtained by com- 
paring its physical properties with those reported for the synthetic 
stereoisomers of 1,2,4-trimethylcyclohexane. Of the four stereoj- 
somers possible, three have been prepared and two identified {3}. 
The properties of these and the properties of the petroleum compound 
are given in table 1. The comparison shows that the 1,2,4-trimethyl- 
cyclohexane from petroleum is not the “cis” stereoisomer boiling at 
146° C and, also, that it differs appreciably in properties from the 
l ciss2cis,4trans-trimethyleyclohexane prepared by Skita and Schneck 
[8]. The compound therefore is probably either 1 15,2 trons 4citti- 
methyleyclohexane or 1 ¢4s,2 trans:4trane-trimethyleyclohexane. In the 
absence of critical data, it is not possible to state which of the two 
it Is. 

TaBLE 1.—Comparison of the properties of 1,2,4-trimethylcycloherane from petro- 

— with those reported for three of the four stereoisomers of 1,2,4-trimethyleyclo- 

exane* 





Boiling 
Density | Refractive 


Compound at 20° C index 
Hg) 





1,2,4-Trimethylcyclohexane from Oklahoma petroleum; White and g/ml 

Glasgow é 0. 7720 
leis,2trene,42-Trimethylcyclohexane; Von Auwers [9] 777 
leis, 2eie,4trane” TTimethylcyclohexane; Skita and Schneck [8] . 786 
leis,2cie,4ete-TTimethylcyclohexane; Skita and Schneck [8] 79 














« The data for the properties of the synthetic compounds have been taken at their face value. Since 
the values were determined on small samples (5 to 10 g) whose purity was not specified, they are not neces: 
sarily accurate values. 

» If the compound prepared by Von Auwers is leis,2trenes4transtrimethylcyclohexane, the fourth, and 
unknown, stereoisomer is 1 cis, 2irane,4eis-trimethylcyclohexane. 
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y. COMPOSITION AND PROPERTIES OF THE MATERIAL 
DISTILLING BETWEEN 98.8° AND 100.8° C AT 215 mm 
He, AFTER PARTIAL REMOVAL OF 1,2,4-TRIMETHYL- 
CYCLOHEXANE 


Figure 6 shows the boiling-point and refractive-index curves of 2.3 
jiters of material obtained by redistilling a mixture of 0.6 liter of the 
higher-boiling mother liquors from the crystallization of 1,2,4-trimeth- 
yevelohexane and the 2 liters of petroleum distillate boiling between 
9g,8° and 100.8° C.2 The comparatively low refractive index indi- 
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FicurE 6.—Fractionation by distillation of crystallization residues and petroleum 
fractions boiling directly above 1,2,4-trimethylcyclohexane. 


The circles on the two curves represent, respectively, the boiling points and refractive indices of succes- 
sve fractions of distillate. 

cated that the redistilled material, while still containing some 1,2,4- 
tnmethyleyclohexane, was largely composed of paraffinic hydrocar- 
bons. The shape of the boiling-point curve suggested the presence 
of at least two paraffins. The indications were that the bulk of the 
paraffin previously encountered as a binary mixture with the isolated 
tnmethyleyclohexane had been concentrated in the distillate boiling 
between 98.8° and 99.8° C, and that another paraffin had accumu- 
lated in the fractions boiling near 100.7° C. The presence of two 
paraffins could not be clearly established, however, until the amount 
of naphthene present had been determined with accuracy. Naph- 


wenn atari boiling above 100.8° C had been set aside for the separation of its chief constituent, 2-methyl- 
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thenes differ distinctly from paraffins in the magnitude of their refra¢. 
tive indices, but the amount of naphthene in the apparently threo. 
component system under consideration could not be determined py 
measuring the refractive index alone, because of a possible differenc, 
in refractive index of the paraffinic material in the successive fractions 
of the distillate. Since, however, the two paraffins had the same 
molecular weight,” the content of the naphthene in the fractions 
could be determined by combustion analysis, and from the knoyy 
refractive index of the distillate and the assumption that the naph. 
thene was 1,2,4-trimethylcyclohexane, a refractive index of the residyg] 
paraffinic material could be computed. Combustion analyses of 
representative fractions are given in table 2. The figures in the last 
column indicate clearly that the paraffinic material in the fractions 
boiling at 99.05° and 99.5° C differed in composition from that boiline 
at 100.4° C. The analysis was confirmed when 4-methyloctane was 
later separated from the material distilling between 100° and 101°C 


TABLE 2.—Properties and composition of fractions of ‘cr, distilling between 99° 
and 100.5° C at 215 mm Hg 





Estimated mole} ,, ,. 
Boiling percentage 
Boiling point of | point of At Refrac- Moles H20 
fraction (215 fraction |(boiling| tive aoa 
mm Hg) (760 mm |point) *| index Moles COs Paraffin] 
Hg) CeHao 





°C °c nit 
141.3 42.3 | 1.4115 1.0740 (1 experiment) 
141.6 42.1 | 1.4101 (1) 1.0862 
(2) 1.0861 
ae 142.4 42.1 | 1.4070 1. 0890 (1 experiment) 


























® Difference between the boiling point at 760 and 215mm Hg. 
» Assuming ni for the naphthene to be 1.4250. 


Attempts were made to separate the lower-boiling paraffin from the 
fractions boiling near 99.5° C by pe Anse agg from dichlorodifluoro- 
a 


methane, and from propane, but, although the distillate contained 75 
mole percent of the paraffinic material, only a tacky mass was obtained 
when the solutions were cooled to low temperatures. At —170° C 
the glue-like material appeared to be insoluble in the solvents used. 
An effort to fractionate the material which boiled at 98.8° C by dis- 
tilling it with excess glacial acetic acid [10, 11] was also unsuccessful. 
Apparently the composition of the fraction was not as simple as wes 
supposed or else the azeotropic properties of the component hydr- 
carbons were too much alike. A system of fractionation by distilla- 
tion at normal pressure, alternated with distillation at a pressure of 
215 mm Hg, was then tried and found to be fairly successful in sep- 
arating the chief paraffinic constituent from the fraction boiling be- 
tween 98.8° and 99.8° C. 


10 It is almost axiomatic that paraffins from petroleum which boil near 140° C are nonanes. 
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2 ,6-Dimethyiheptane 93.4 135.2 7 
Nonanaphthene (see text) 93.6 136.7 
1,3,5-Trimethylcyclohexane (95.0) (138.0) 
2,3- Dimethytheptane 98.8 1407 
1,2,4-Trimethylcyclohexane 97.8 141.2 
4- Methyloctane 100.4 1424 
2- Methyloctane 101.3 143.3 _J 
3-Methyloctane 102.2 144.2 
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Figure 7.—Calculated vapor pressure curves of certain hydrocarbons found in Oklahoma petroleum. 


The numbers in-cireles refer to the curves of the numbered hydrocarbons in the table at the lower right of the figure. The table gives the boiling points, at 760 mm and at 215 mm Hg, of the 
listed hydrocarbons. Those without parentheses are experimentally determined values, those in parentheses are values estimated from the boiling points of impure fractions. For the mode of 
construction of the curves, see text. 
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vy]. SEPARATION OF PARAFFIN-NAPHTHENE MIXTURES 
BY ALTERNATE DISTILLATION AT NORMAL AND RE- 
DUCED PRESSURES 


As pointed out in section III, 1, 2, 4-trimethylcyclohexane has a 
higher volatility at a pressure of 215 mm Hg than the paraffin asso- 
ciated with it in the petroleum, and a correspondingly lower volatility 


at 760 mm Hg. ‘This behavior illustrates the general principle that 
with paraffins, naphthenes, and aromatic hydrocarbons which boil 
at nearly the same temperature, members of one class differ materially 
from those of another class in the change of their vapor pressures 
with temperature. Because of this, a binary mixture consisting of a 
paraffin and a naphthene boiling close together can have the partial 
ressures of its components so altered as to enhance their separation 
c fractional distillation, or to reverse the order in which they distil 
(is). Figure 7 shows calculated vapor-pressure curves of the paraffin 
and naphthene hydrocarbons in the fraction of this petroleum which 
normally boils between 135° and 145° C. In the lower right of the 
feure the boiling points at 215 mm and 760 mm Hg are given for the 
hydrocarbons." The values for boiling points were experimentally 
determined, except those in parentheses. These were estimated from 
the boiling points of impure fractions by correcting for the effect of 
the major impurity. To construct the vapor-pressure curves of 
figure 7, the reciprocals of the absolute boiling points of 1,2,4-tri- 
methyleyclohexane and of 2-methyloctane were first plotted with respect 
to the logarithm of the corresponding pressures for the two experi- 
mental points (at 215 and 760 mm Hg). Intermediate values for 
boiling points and pressures were then obtained from the straight 
line connecting the experimental points. With these additional 
values, the vapor-pressure curves of the two hydrocarbons were 
drawn. Lines for the other paraffins were drawn with the same 
curvature as that for 2-methyloctane, and those for the naphthenes 
wih the same curvature as that for 1,2,4-trimethylcyclohexane. 
This was permissible because the difference between the boiling 
points at 760 and 215 mm Hg for each of the paraffins was quite 
close to the mean value of 42° C, and for each of the naphthenes, 
close to the mean value of 43.2° C. It should be noted that the 
curves as constructed may not correspond exactly with the actual 
vapor pressure curves of the individual hydrocarbons, but they 
nevertheless illustrate the relative change in vapor pressure with 
temperature of the two classes of hydrocarbons. All of the curves 
have been useful in the analysis of the distillate boiling between 135° 
and 145° C. Curves 1 and 2, for example, indicate the substantial 
difference between the boiling points at normal pressure of 2,6- 
dimethylheptane and a nonanapthene, and the close proximity of 
,|A practical demonstration of the possibilities of the principle is afforded by the experiments of S. T. 
Schicktanz mentioned in footnote 4. The results of the experiments were kindly called to the attention of 
the authors by Schicktanz, prior to the present investigation. ‘ : 
Ne he hydrocarbons used for the boiling-point determinations were, with the exception of 4-methyloctane, 
the best fractions prepared from petroleum. The sample of 4-methyloctane was supplied by George Cal- 
= hag eae) Concuine Corporation. Properties of the samples used are listed in the following papers 
‘ng to 2,6-dimethylheptane, J. Research NBS 17, 952 (1936) RP955; nonanaphthene, 134.7° C, J. Re 
sarch NBS 18, 808 (1934) RP745; isononane (probably 2, 3-dimethylheptane), 140.7° C, described in this 
paper; 1,2,4-trimethyleyclohexane, described in this paper; and 4-methyloctane (Calingaert), 2-methyl- 
octane, and 3-methyloctane, J. Research NBS 19, 423 (1937) RP1033. The fractions containing the nona- 
naphthene (probably 1,3,5-trimethyleyclohexane) boiling at 138° C, and the nonanaphthene boiling at 


i C, had refractive indices at 25° C of 1.4233 and 1.4263, respectively. The isolation of these compounds 
tom the petroleum has not been completed. 
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their boiling points at 215 mm Hg. It would be expected that , 
mixture of these two hydrocarbons could be separated by fractiong| 
distillation at 760 mm Hg, but not at 215 mm, at which pressure g 
constant-boiling mixture would result. This was found to be the 
case when a fraction of the petroleum composed of these two hydro. 
carbons was resolved into its constituents [12]. Similarly, use was 
made of the pressure-temperature relationships depicted in curves ¢ 
and 9, in separating 3-methyloctane from the petroleum [3]. Special 
attention, however, is called to curves 4 and 5, which concern the 
present work. Curve 4 crosses curve 5 at a pressure of about 650 mm 
In this region of pressure the pair of hydrocarbons form a mixture 
which is constant-boiling at about 135° C. At higher or lower pres. 
sures, the mixture would distil over a definite temperature range, 
It follows that constant-boiling mixtures of naphthene and parafliy 
hydrocarbons are only incidental, and that such mixtures may be 
resolved into their components by distillation at a pressure somewhat 
removed from that causing the constant-boiling behavior. Closer 
inspection of curves 4 and 5 indicates that a binary mixture of 1,24. 
trimethyleyclohexane and the “2,3-dimethylheptane” could be sepa- 
rated more readily by distillation at 215 mm Hg than at 760 mn, 
The present problem, however, was to separate a three-component 
mixture which contained, besides these two compounds, 4-methyl- 
octane (curve 6). It was believed that by fractionally distilling the 
mixture at atmospheric pressure, the 4-methyloctane could be removed 
by rejecting the end fractions and residue which, from the order of 
distillation of the constituents, would contain the 4-methyloctane 
and some of the naphthene. The bulk of the lower-boiling paraffin 
would be concentrated in the initial fractions as a binary mixture 
with the naphthene which afterwards could be removed by fractiona- 
tion at a pressure of 215 mm Hg. 


VII. SEPARATION OF THE PARAFFIN-NAPHTHENE/IMIx. 
TURE BOILING BETWEEN 98.8° AND 99.8° C AT 215 mm Hg 


Using the procedure outlined above, the 1600-ml fraction shown in 
figure 6 to boil between 98.9° and 100° C was mixed with 400 ml of 
paraffinic mother liquors obtained from further crystallization of the 
1,2,4-trimethylcyclohexane stock, and redistilled at 760 mm. The 
boiling-point and refractive-index curves of the distillate are shown in 
figure 8. The individual points denote the boiling points and refrac- 
tive indices of successive 50-ml fractions. The curves show that the 
first liter of distillate was almost constant-boiling at 141.0° to 141.1°C, 
and changed but little in refractive index, 1.409 to 1.410. Both prop- 
erties indicate a constancy in composition of the material. The 50-nl 
fraction taken after 900 ml had distilled (bp 141.1° C, n?=1.4101) 
contained 20 mole percent of a naphthene and 80 mole percent of 8 
nonane. From the composition, the nonane was calculated to have 
an index of 1.4065 at 25° C. From this point in the distillation, the 
remaining material distilled at increasingly higher temperatures aud 
exhibited a steady increase in refractive index. For the 50-m! 
fraction which marked the 1,750-ml point in the distillation (bp 141.9 
C, n?=1. 4130), the naphthene content was 40 mole percent. The 
calculated refractive index of the paraffinic material in this fraction 





white | 1,2,4-Trimethyleyclohexane and Isononane from Petroleum 149 
was 1.4046 at 25° C as compared with 1.4036 for pure 4-methyloctane. 
From the analysis it was evident that the naphthene and 4-methyloc- 
‘ane had collected in the end fractions of the distillation, and for this 
reason the last 400 ml of distillate and the residue were rejected. 
The first 1,400 ml to distil, which was rich in the lower-boiling paraffin, 
was redistilled at 215 mm Hg. 

As a result of this second distillation at reduced pressure, much of 
the naphthene still present collected in the first portion to distil, and 
most of the 4-methyloctane which remained was held in the residue. 
The 850 ml of oil representing the middle portion of the distillate 
consisted chiefly of the desired isononane. The individual fractions 
comprising the middle portion ranged in boiling point from 98.7° 
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VOLUME OF DISTILLATE IN LITERS 


Ficure 8.—F ractionation by distillation at normal pressure of the petroleum fraction 
enriched in 2,3-dimethylheptane. 


The points on the curves indicate the boiling points and refractive indices of successive fractions, 


to 98.9° C and varied in refractive index (in the order distilled) from 
14125 to 1.4087. By combustion analysis, the naphthene content 
was 30 mole percent for a fraction with nj=1.4122 and 15 mole 
percent for one with n?3=1.4087. The paraffinic material in most 
of the distillate had a calculated refractive index of 1.406 at 25° C. 
In the last fractions, however, the index of the paraffinic material had 
decreased to 1.405. 

As a final step in the separation of the isononane, the middle 
portion (850 ml) of the distillate described above was redistilled in 
the same column, at 760 mm Hg. The boiling points of the successive 
fractions from this final distillation were measured at 760 mm and 
at 215 mm Hg in order to determine the nature of the distillate. 
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The boiling-point curves for the two sets of data are shown in the 
lower part of figure 9. The curve at the top of the figure indicatgs 
the range in refractive index of the distillate. Barring the firs 
100 ml of oil distilled, the gradual rise in refractive index of the 
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distillate, together with the divergence of the boiling-point curves, 
shown by the steady increase in boiling point at 760 mm Hg (middle 
curve) and the steady decrease in boiling point of the fractions when 
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mite 1,2,4-Trimethyleycloherane and Isononane from Petroleum 151 
measured at 215 mm Hg (bottom curve), is regarded as substantial 
evidence that the major portion of the distillate was composed of a 
binary mixture of the unidentified nonane and 1,2,4-trimethylcyclo- 
hexane. Ae 

Although the mean refractive index of the material, 1.4095, signified 
that the isononane was the chief constituent, none of the distillation 
fractions would freeze. A 240-ml portion, characterized by a con- 
stant boiling point of 140.9° C and a constant refractive index of 
1.4090, was subjected to a fractional crystallization from its solution 
in propane and methane [13]. Apparently a crystalline phase formed, 
but when separated and freed from solvent it had the same physical 
properties as the oil in the mother liquid, indicating that no fractiona- 
tion had taken place." 


VIII. PROPERTIES OF THE ISOMERIC NONANE FRACTION 
AND THE PROBABLE IDENTITY OF THE NONANE 
WITH 2,3-DIMETHYLHEPTANE 


The properties of the recombined material from the crystallization 
experiment are as follows: Boiling point at 760 mm Hg, 140.80° C; 
congeals to a glass at —119° C; density at 20° C, 0.7340 g/ml; re- 
fractive index, n=1.4121; critical solution temperature in aniline, 
71°C. From two determinations by combustion, the ratio of moles 
of H,O to moles of CO, was found to be 1.0950+0.0002. From this it 
was estimated that the petroleum fraction contained 85.5 mole percent 
of the nonane and 14.5 mole percent of 1,2,4-trimethyleyclohexane. 

The properties of the constituent isononane, caheglatad from those 
of the petroleum fraction by correcting for the effect of the 1,2,4- 
trimethyleyclohexane present, are given in table 3. That the com- 
pound is related to, if not identical with 2,3-dimethylheptane, is 
indicated in the table by the similarity in their physical properties. 
Of the 16 isomeric nonanes which have been prepared (chiefly the 
substituted octanes and heptanes) 2,3-dimethylheptane is the only 
one which has properties similar to those of the petroleum hydro- 
carbon. However, in the absence of specific information concerning 
ihe properties of other possible isomers of nonane, it cannot be stated 
definitely that the two compounds are identical. General knowledge 
relating the properties of paraffin hydrocarbons to their molecular 
structure [15] indicates that only a few of the 19 unreported nonanes 
can be closely related to the nonane in question. From such a correla- 
tion, all of the trimethylhexanes and tetramethylpentanes may be ex- 
pected to boil below 140°C and all of the dimethylethylpentanes may be 
expected to have densities greater than that of the isolated hydrocarbon. 
Some of the methylethylhexanes, however, as well as some of the di- 
nethylheptanes, can be expected to have properties similar in magni- 
tude to those of the nonane from petroleum. These include, 3,5-, 3,4-, 
and 4,4-dimethylheptane and 3,3-, and 3,4- methylethylhexane. 

‘A method was considered for purifying the material by dehydrogenation of the naphthene present, 


lollowed by adsorption of the aromatic product, but time was not available to make an adequate study of its 
possiblities, 
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TABLE 3.—Comparison of the properties of an isomeric nonane from petrol 
with those of 2,3-dimethylheptane 


eum 





Boiling 
point Behavior Density 
(760mm | on cooling | at 20°C 
Hg) 


Nonane, CoH 





°C , g/ml 

Petroleum fraction, corrected for naphthene content-__- 140.7 0. 727 
glass at 
2,3-Dimethylheptane, Whitmore and Southgate ® 6140. 65 —116.7 . 7235 

















® Private communication from F. C. Whitmore and H. A. Southgate, Pennsylvania State Co!] 
» Initial boiling point. dp/di=0.049. 


IX. CONTENT OF 1,2,4-TRIMETHYLCYCLOHEXANE ANp 
THE ISONONANE IN THE PETROLEUM 


Of the 3 liters of petroleum shown in figure 2 to boil between 
98° and 99° C at 215 mm Hg, nearly 2 liters was 1,2,4- trimethy!- 
cyclohexane and somewhat less than a liter was the isomeric nonane. 
It is estimated that the percentage of 1,2,4-trimethyleyclohexane in 
the petroleum is about 0.1, and that of the isononane about 0.05, 
by volume." 
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\ See references [16, 17] for a discussion as to the accuracy of these estimates. 
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SUMMARY OF AN INVESTIGATION OF THE COMPOSITION 
OF A MIDCONTINENT PETROLEUM DISTILLATE, 
BOILING BETWEEN 100° AND 130° C.? 


By Robert T. Leslie ? 


ABSTRACT 


The results of 7 years of investigation of the composition of the naphtha fraction 
f petroleum boiling between 100° and 180° C are summarized. Of the 130 
liters found to boil in this range after the first laboratory distillation, the constitu- 
ents of about one-sixth of the volume were identified or characterized by physical 
properties, isolated or concentrated, and their concentration in the petroleum 
roughly estimated. About one-fourth of the volume remains as unidentified 

stillate. This is composed in part of material of which the constitution is only 

dicted (about one-twelfth of the volume) and of additional quantities of the 
nstituents which were isolated previously (about one-sixth of the volume). 
The remainder was lost during the investigation or distilled out of the range on 
ibsequent distillations. Photomicrographic studies were made of those con- 
centrations of distillate which could be crystallized. There is given a bibliography 
of the physical constants of hydrocarbons reported to boil within this range. 


CONTENTS 


I. Introduction 
I] ae ly of the composition of the distillate. 
. Distillate between 98° and 102° C 
. Distillate between 103° and 106° 
3. Distillate between 108.5° and 110.5° ¢ 
4. Distillate between 112° and 116° C_- 
5. Distillate between 117° and 124° - es 
§. Distillate between 125° and 126° 
7. Distillate between 127° and 130° oC ; 
III. Probs ible composition of a midcontinent petroleum boiling between 100° 
and 130° C 


I. INTRODUCTION 


This paper summarizes the work done to date since 1929 with the 
portion of midcontinent petroleum distilling between 100° and 130° C 
and includes some new information obtained since the publication of 
‘previous summary paper * and two subsequent specialized papers.‘ ® 


eer es 
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About 6 percent of the total volume of the crude petroleum js. 
tilled between 100° and 130° C and, of this, by far the largest yolyn, 
distilled between 117° and 128° C. Curves A; and A, of figure | 
show, respectively, the refractive indices and the distribution 9 
volumes after four distillations in the laboratory, before the remoyy| 
of any of the constituents. Curves B, and B, of the same figure shoy 
the distillate still remaining after the removal of seven hydrocarhoy: 
and the distillation of the residual material. This is composed of the 
unidentified constituents and of unisolated residues of the identified 
compounds. Although unidentified constituents which boil chiefy 
between 101.5° and 116° C still remain, the constituents present in thp 
largest proportions have been removed. The paper includes, besides 
a summary of the compounds actually isolated, all the informatio, 
concerning the material represented by curves B, and B, of figure |. 
In addition, the values of physical constants for the hydrocarboys 
reported in the literature to boil within the range are given in table |, 
This information was obtained by surveying the literature up to 1937 
including standard works of reference. The values were selected from 
publications dealing with the preparation of pure compounds or ap. 
parently reliable measurements of the properties of nearly pure con. 
pounds. Where necessary, the only available data were included. 
The best values, in the author’s opinion, have been plotted in figure |, 
and are numbered for convenient reference to the table.® 


¢ Figures in brackets throughout the paper are the reference numbers given in figure 1 and table 1. 
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II. STUDY OF THE COMPOSITION OF THE DISTILLATE 
1. DISTILLATE BETWEEN 98° AND 102° C 


In earlier work on the portion of petroleum boiling between gy: 
and 100° C, a large quantity of a nearly constant-boiling mixture o; 
n-heptane and methylcyclohexane was found just below 100° ¢ 
The material also contained, as did all the distillate boiling betwee, 
95° and 108° C, a considerable quantity of toluene. 

Further systematic distillation, after the removal of the major por- 
tion of these constituents,’ * probably resulted in the reappearance a} 
about 98° and 102° C of the concentration shown by curve B, of 
figure 1. The maximum volume was found to distil between 100° to 
100.5° C, and this coincides with the temperature at which a maximuy 
in the refractive-index curve (1.4227) appears. The boiling range an¢ 
refractive index correspond very nearly to the properties of methyl- 
cyclohexane [3]. The photomicrographs of crystals of the material 
boiling at 101.2° C (fig. 1) support this conclusion.® 

The material distillmg just below this peak shows a rapidly de- 
creasing refractive index, the minimum index (1.3925?) occurring at 
97.5° to 98° C. This indicates that the lower-bo'ling constituent of 
the large-volume fraction was chiefly a paraffin hydrocarbon. Table 1 
lists two hydrocarbons, n-heptane [1] and 2,2,4-trimethylpentane 
[2], which are reported to have approximately the properties of this 
fraction. Comparison of photomicrographs of the crystals from the 
material at 98.3° C (fig. 1) with photomicrographs of pure n-heptane 
and 2,2,4-trimethylpentane show that the chief constituent was 
n-heptane.” 


2. DISTILLATE BETWEEN 103° AND 106° C 


No great effort was made to separate the constituents other than 
toluene from the material boiling between 102° and 106° C, because 
of the relatively small quantities present and because the fractions of 
distillate did not crystallize readily. 

Curve B, of figure 1 shows about 4 liters of material boiling between 
102.5° and 105.5° C after the toluene was removed. The fact that 
the refractive-index curve slopes down more gradually above 101° C 
than it does below that temperature indicates either the presence o/ 
one or more naphthenic hydrocarbons boiling just above methyleyclo- 
hexane, such as ethylcyclopentane [4], 1,1,2-trimethyl-2-ethyleyclo- 
propane [5], or 1,2-dimethyl-1-ethyleyclopropane [6], or the presence 
of a paraffin, such as 2,2,3,3-tetramethylbutane [7] or 2,2-dimethyl- 
hexane [8], which tended to form a constant-boiling mixture with 
methyleyclohexane. At 105° to 105.5° C the refractive index has 
dropped to a value which shows that the distillate was largely parafiine. 
Attempts to crystallize fractions in this range by cooling without 
dilution with solvents resulted in clear glasses. When a fraction dis- 
tilling at 105.2° to 105.3° C was diluted with propane and poured into 
liquid methane, cooled by liquid nitrogen, two liquid layers appeared. 
These layers were separated by using a separatory funnel constructed 

7M. M. Hicks-Bruun and J. H. Bruun, BS J. Research 8, 525 (1932) RP432. 


8 J. H. Bruun, R. T. Leslie, and S. T. Schicktanz, BS J. Research 6, 363 (1931) RP280 
#10 R. T. Leslie, and W. W. Heuer, J. Research NBS 18, 639 (1937) RP1000. 
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Leslie} 


as shown in figure 2. The residue obtained from each layer, after 
poiling off the solvents, showed no difference in refractive index or the 
tendency to crystallize. Temperatures lower than that of liquid 
nitrogen might have caused crystallization of one of the layers. 

From a consideration of the properties of the 5 hydrocarbons known 
to boil between 102° and 108° C, ethyleyclopentane [4] and 2,2-di- 
methylhexane [8] are most likely to occur in this small volume. 





}, DISTILLATE BETWEEN 108.5° AND 110.5° C => 


Though the boiling point of toluene is given 
i, the literature as about 110.8° C, very little 
of it was found in this fraction. The maxi- 
mum concentration of the aromatic hydro- 
carbon was actually found to correspond with 
‘he maximum in the refractive-index curve, 
4,, which oceurs at about 105° C. The ma- 
terial represented by curve B, contained prac- 
tically no toluene, as evidenced by the fact 
that the refractive index was not changed by 
treatment with nitrating agents or filtration 
through silica gel. The fact that the refrac- 
tive-index curve in this range rises with boil- 
ing point indicates, therefore, that the dis- 
tillate contained paraffinic mixed with higher- 
polling naphthenic constituents. 

When cooled, none of the fractions could be 
crystallized without solvents, and distillation 
with acetic acid failed to separate the naph- 
thene, as shown by figure 3. Curves A, and 
4, of this figure show the results obtained 
ifter about 15 distillations without acetic acid 
in stills ranging in type of construction from 
\2 bubbling plates to 3-meter chain-packed 
columns and curves B, and B, show the results 
ptained after two successive fractional dis- 
ir with acetic acid. The meg was 
distilled in three batches of about 1.5 liters , - 
ech, and the distillate was again divided gg gn yen 
according to refractive index for the second liquid layers formed 
distillation. It is to be observed that curves when hydrocarbons are 
indicate very little improvement in separation. — @/uled with liquid meth- 

: i ane and propane. 

Attempts to crystallize the material from 
methane and propane as solvents, as de- 4, <rainorifice with ground stop- 
bimhad ¢ . : * Sena per. Funnel removed from 
scribed above, resulted in two liquid layers, bath of liquid nitrogen toa cold 
oue of which was chiefly methane and the {st tubeand lower layer drain- 
dther propane, each containing petroleum ; 
lydrocarbon. The layers were separated but no differences in the 
voperties of the fractions were observed. 

The 2,5- [9], 2,4- [11], 3,3- [16], 3,5-dimethylhexanes [15], and 2,2,3- 
imethylpentane [14] are shown in figure 1 as possible paraffinic 
foustituents. The fact that the refractive-index curve rises with 
boiling point indicated that the naphthenic constituents were the 
ugher-boiling components of the mixture. It appears that 1,2,3- 
7], 1,1,3- [24], 1,1,2- [23], and 1,2,4-trimethyleyclopentane [19] 
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were the most likely possibilities, though their boiling points gp, 
higher than those of 1-methyl-1, 2-diethyleyclopropane [10] ala 
methyl-2-isobutyleyclopropane [12]. Since the total volume of 
material distilling between 108° and 111° C was less than 4 liters, it on 
estimated that the constituents of this fraction originally constityte; 
not more than 0.1 percent of the crude petroleum. Of this, the me 
constituents appeared to be paraffinic. It is possible that the dif 
culty in separating the constituents sufficiently to obtain crystals may 
have been caused by the presence of a number of the isomeric octanes 
with nearly the same boiling points. , 


4. DISTILLATE BETWEEN 112° AND 116° C 


This small volume of distillate, less than a liter, had a range of 
refractive indices which indicated that it was chiefly naphthenic in 
composition. The most likely constituents are the trimethyleyelo. 
pentanes [17, 19, 23, 24] mentioned above. Attempts to crystallize 
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BOILING RANGE IN °C 
Figure 3.—Volume and refractive index of material distilling between 108° and 
110° C before and after distillation with acetic acid. 


A, April 1935, before distillation with acetic acid; B, December 1935, after distillation with acetic scid. 


this material with methane and propane resulted in two liquid layers. 
Additional careful distillation should concentrate this fraction to the 
point where it will crystallize, unless it contains small quantities of 
many of the possible hydrocarbons. 


5. DISTILLATE BETWEEN 117° AND 124° C 


Because the proportion of the material which distilled between 11” 
and 124° C was large, most of the effort was devoted to separatilg 
constituents in this range. 

The refractive-index and volume curves, A; and A, of figure |, 
show that the fraction distilling between 114° and 128° C contained 
originally a large quantity of paraffinic material which had a slight 
tendency to separate into two large volumes. An amount of n-octane 
[44] equivalent to about 1 percent of the crude petroleum was remove 
by crystallization from the distillate between 123° and 126° C." Sub 


u R, T. Leslie and 8. T. Schicktanz, BS J. Research 6, 377 (1931) RP282, 
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sequent distillations resolved the fraction into two more or less dis- 
tinct ones boiling between 116° and 124° C and between 124° and 
128° C. The lower-boiling fraction appeared to have the more com- 
plex composition. In previous publications there has been described 
in detail the isolation of a relatively large quantity of 2-methylhep- 
‘ane ” [26] (about 0.5 percent of the petroleum) boiling at about 117° 
(. a smaller quantity of a mixture of naphthenes * (about 0.2 per- 
cent) boiling between 119.2° and 120.8° C, and of a small quantity 
of a naphthene (about 0.04 percent), boiling at about 123.4° C. 

The photomicrographs at the bottom of figure 1 show the crystal 

behavior of fractions from this range. The photograph of the mate- 
rial boiling at 116.9° C shows material identified as 2-methylheptane 
26]. Those of the fraction at 119.8° C probably show 1,1-dimethyl- 
cyclohexane [34]. The three photographs of the same fraction were 
made successively and show two crystal forms.'® At 120.2° C the 
crystals still resemble those at 119.8° C. At 120.3° C, however, the 
crystals resemble m-dimethyleyclohexane * [36]. The naphthenic 
material boiling at 123.4° C, thought to be trans-o-dimethyleyclohex- 
ane [39], differs in appearance from the crystals of the lower-boiling 
naphthenes. The fine crystals which appear in one of the photo- 
craphs may be a small amount of n-octane. The result of adding 
more octane from the higher-boiling fractions is shown by the photo- 
graph of the ‘‘mixture.”’ 
' The residue shown by curve B, between 116° and 122° C repre- 
sents the relatively small quantity of material which concentrated 
after removing, by crystallization as completely as practicable, the con- 
stituents discussed above. The refractive-index curve B,, figure 1, 
indicates that the material contained lower-boiling paraffinic constitu- 
ents and higher-boiling naphthenes, which were probably an addi- 
tional small quantity of the hydrocarbons mentioned above. 

In addition to the 2-methylheptane [26], it is possible that this frac- 
tion also contained small quantities of 3-ethylhexane [29], 3,4-di- 
methylhexane [27], 3-methylheptane [28], 4-methylheptane [31], and 
3-methyl-3-ethylpentane [32]. If they were present they would be 
found in the residue shown by curve B, of figure 1. Of the seven 
naphthenes which are reported to boil in this range [30, 33, 34, 35, 36, 
39, and 40], three were found and tentatively identified [34, 35, 39]. 
Cycloheptane [30], because of its high freezing point, and p-dimethyl- 
cyclohexane [33], because of the absence of the characteristic crystals 
which it forms, appeared to be absent. 1-methyl-3-ethyleyclopentane 
35] and isopropyleyclopentane [40] may be present in small quantities 
but were not detected. 


6.“DISTILLATE BETWEEN 125° AND 126°FC 


The approximately 7 liters of distillate obtained in the narrow 
range, 125° to 126° C, as shown by the volume curve B, of figure 1, 
was probably n-octane, as indicated by the minimum in the refractive- 
index curve B,;. The photomicrographs of crystals from these frac- 
tions show that the chief constituent behaves like n-octane "7 [44]. It 
is possible, however, that small quantities of 2,2,4,4-tetramethylpen- 


TT 
- Leslie, BS J. Research 10, 609 (1933) RP552. 
oe eslie, J. Research NBS 18, 41 (1935) RP808. 
wi aL: Leslie, J. Research NBS 17, 761 (1936) RP943. 
"4. T, Leslie and W. W. Heuer, J. Research NBS 18, 639 (1937) RP1000. 
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tane [37], and 2,2,5-trimethylhexane [38], were also present since their 
boiling points are reported in this neighborhood. From the manner 
in which the refractive-index curve rises abruptly on either side of 
this narrow boiling range, it is apparent that the material also eo) 
tained naphthenic constituents in small quantities, which were prob. 
ably some of the cis-dimethylcyclohexanes [43, 46], 1-me thyl-2. 
ethyleyclopentane [41], or isopropyleyclopentane [40]. 

125° to 126° C distillate must have contained only a oe percentare re 
of naphthenic material, since the refractive- index of the large-volume 
fraction is very low and there are only small quantities of the frac. 
tions on either side which have high refractive indices. 


7. DISTILLATE BETWEEN 127° AND 130° C 


The material represented as ‘“‘residue”’ in figure 1 was composed of 
“tailings” from several systematic fractionations and probably con- 
tained considerable cracked material. Treatment with silica gel, or 
other reagent for removing olefins, and distillation would return most 
of this material to lower-boiling fractions. Very little material acty- 
ally distilled between 127° and 130°C. This might be expected since, 
as shown by figure 1, very few hydrocarbons are reported to boil iy 
this range. If 2,3,5-trimethylhexane [45], cis-o-dimethyleyclohexane 
[46], and n-propyleyclopentane [47] were present, they must have been 
minor constituents. The next higher-boiling concentration of a con- 
stituent occurred at about 130° C, and the isolation of ethyleyelo- 
hexane from it has been reported by Rose and White.'® 


III. PROBABLE COMPOSITION OF A MIDCONTINENT 
PETROLEUM BOILING BETWEEN 100° AND 130° C 


Most of the conclusions which have been reached in the investiga- 
tion are summarized in table 2. This table presents, in addition to 
the relative quantities of constituents actually isolated, conclusions as 
to the probable presence or absence of each of the other hydrocarbons 
boiling between 98° and 131° C on which the author has found data 
in the literature. A column of the table also gives the approximate 
boiling ranges of the distillate from which constituents were isolated 
or in which constituents whose presence is uncertain might be ex- 
pected to occur. 


18 F, W. Rose, Jr., and J. D. White, J. Research NBS 15, 151 (1935) RP817. 
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TaBLE 2.—Known hydrocarbons boiling between 100° and 130° C (approximately) 
“with relation to the corresponding naphtha fraction of a midcontinent petroleum 





Refer- Relative 
ence Pee tage — 
num-} ; +3 ; oiling range ‘ound, 
Compound ber ee in which ma- | referring 
to Pp I terial occurred, | to n-oc- 
table if present tane as 
1 unity! 


Occur- 
rence 2 





°C °C 
98. 4 95 to 190 
99. 3 
101.2 98 to 102 
103.0 +0. 5 98 to 102 
103. 5 to 104.5 


n-Heptane 

2, 2, 4-Trimethylpentane- - - - 
Methyleyclohexane 

Ethyleyclopentane 

1,1, 2-Trimethyl-2-ethyleyclopropane - 


|, Dimethyl-1-ethyleyclopropane-.-_- 
9 3, 3-Tetramethyl butane 


104 to 105. 2 
9. 9.-Dimethylhexane.........------ 106 to 107 102 to 106 

2, §-Dimethylhexane 109. 2 +0.1 108. 5 to 110.5 
L-Methyl-1, 2-diethyleyclopropane- --- 108 to 109 


COND AP wrre 


2, 4-Dimethylhexane 108 (750 mm) 
|-Methyl-2-isobutyleyclopropane 109. 5 to 110.5 
Toluene... ‘ ; 98 to 111 

2, 2, 3 Trimethylpentane 108. 5 to 110.5 
q-3, 5-Dimethylhexane 108. 5 to 110.5 
3, 3-Dimethylhexane 111 to 112 108. 5 to 110. 5 
Trimethyleyclopentane- - ._..-. 7 | 111.5 to 114 112 to 116 
Trimethylpentane.-....----... 112.8 108. 5 to 110.5 
Trimethyeyclopentane 112. 5 to 113 112 to 116 
Trimethylpentane 113.6 112 to 116 


oe 
2, 3 

2, 3, 4 
2, 4 
, By 3° 


2, 3s Dimethylhexane 113.9 112 to 116 
2-Methyl-3-eth ylpentane 2/114 112 to 116 
1,1, 2Trimethyleyclopentane. --...-- 113 to 114 (749 mm) 112 to 116 
1,1, 3-Trimethyleyclopentane---__... 115 to 116 112 to 116 
3, 3-Dimethylbicyclo-(0, 1, 3) hexane__ 115 


2Methylheptane 117.2 115 to 118 
3, 4 Dimethylhexane 118.7 116 to 122 
}-Methylheptane 116. 5 (46 mm) 116 to 122 
$-Ethylhexane__.._- 118. 8 to 119 (766 mm) 116 to 122 
118 to 120 


118 116 to 122 
118. 5 to 119 (750 mm) 116 to 122 
119. 63 116 to 122 
1, 1-Dimethyleyclohexane........... 119.8 116 to 122 
1-Methyl-3-ethyleyclopentane 120. 5 to 121 (756 mm) 116 to 122 


trans-1, 3-Dimethyleyclohexane 120. 40 118 to 1223 
2, 2, 4, 4-Tetramethylpentane 4 120 to 126 
2, 2, 5-Trimethylpentane 120 to 126 
trans-1, 2-Dimethyleyclohexane 122 to 124 
Isopropyleyclopentane 122 to 126 


1-Methyl-2-ethylcyclopentane...___.. 124 122 to 126 
1, 1, 3, 3-Tetramethyl-2, 4-diethylcy- 

clobutane 124 to 125 
tis-1, 4 Dimethyleyclohexane 124. 59 122 to 126 
n-Octane 125. 59 122 to 126 
129 


cit-l, 2-Dimethyleyclohexane--_...._- 130. 04 
n-Propyleyclopentane 130. 6 +0. 2 




















'The numbers in this column give the estimated relative amounts by volume of the given hydrocarbons 
in the petroleum referred to normal octane as unity. To obtain the approximate magnitude of the per- 
‘entage of hydrocarbon in the crude, these numbers should be multiplied by a factor which is roughly 
‘timated, on the basis of the ‘‘over-all’’ losses, to lie between 1 and 134 (see abstract). Similarly, the 
approximate magnitude of the percentage in the 100° to 130° C fraction may be obtained by multiplying 
he numbers by a factor estimated to lie between 17 and 30. (Estimates in previous papers were made by 
foubling the actual quantities found in an attempt to compensate for losses. This factor is probably low.) 
+ Letters in this column have the following significanes: P, present; A, presence in more than small 
amounts considered unlikely because of the properties of the compound; S, possibly present in amounts 
too small to be detected by the method used; U, presence considered unlikely because of structure or small 
quantity of distillate at its boiling point. ; 

This fraction also contained an unidentified naphthene boiling at about 119.8° C (possibly No. 34). 
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Of the distillate from the crude petroleum boiling between 100° anq 
130° C, n-octane [44] was by far the chief constituent. Reasonably 
careful distillation caused most of the n-octane to concentrate betwee) 
122° and 126° C. From this distillate a quantity of nearly pure 
hydrocarbon was isolated equal to about one-tenth of the volume 
with which the investigation was begun. Toluene [13] was the solo 
aromatic constituent and was present in a relatively large proportion 
Distillation caused it to appear in a wide temperature range (98° {, 
111° C), and a quantity was isolated equal to about one-third the 
volume of n-octane. Methylcyclohexane [3] about equal in quantity 
to toluene was found in the distillate between 98° and 102° GC. \ 
quantity of 2-methylheptane [26] amounting to about one-fifth of the 
n-octane occurred in the material distilling between 115° and 118° ( 


At least three cycloparaffinic hydrocarbons were present in appreciable 
quantities in the distillate between the boiling range of the materia! 
containing the 2-methylheptane and that of the n-octane (trans-13- 
dimethylcyclohexane [36], trans-1,2-dimethyleyclohexane [39], and an 
unidentified hydrocarbon*which#may be _ 1,1-dimethylceyclohexane 
[34]). The other known hydrocarbons boiling in this range were in 
relatively small amounts, if at all. Between 126° and 130° C there 
was only asmall quantity of distillate, indicating insignificant amounts 
of the few hydrocarbons whose boiling points are within the range, 
There was, however, distillate other than toluene between 102° and 
115° C which though small, showed the presence of di- and trimethy! 
isomers of octane [8, 9, 11, 14, 15, 16, 18, 20, 21, 22] and possibly one 
or more of the cyclopentanes [4, 17, 19, 23, 24]. 


WasaineotTon, October 21, 1938. 
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GASOMETRIC METHOD AND APPARATUS FOR THE ANAL- 
YSIS OF MIXTURES OF ETHYLENE OXIDE AND 


CARBON DIOXIDE 
By Joseph R. Branham and Martin Shepherd 


ABSTRACT 


A gasometric method for the analysis of ethylene oxide mixed with carbon 
dioxide is reported. Ethylene oxide is swiftly and quantitatively removed from 
the gas phase by a relatively small amount of sulfuric acid. This reagent absorbs 
over 5,000 times its own volume of ethylene oxide, and the quantity used (only 
0.2 ml) dissolves no significant amount of carbon dioxide. The carbon dioxide is 
absorbed in a concentrated solution of potassium hydroxide. A reproducibility 
of +0.05 percent is attained; the accuracy is commensurate. All gases are 
measured dry in & new apparatus designed to perform this type of analysis 
rapidly. 
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I. INTRODUCTION 


_ The use of ethylene oxide as a fumigant has increased considerably 
in this country and abroad during the last 8 — This substance 


is generally used with about 90 percent of carbon dioxide to produce 
noninflammable mixture. The analysis of this simple mixture has 
caused some difficulty, in spite of the remarkable chemical reactivity 
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of the ethylene oxide. The methods so far reported have bee, 
limited in number and are not satisfactorily accurate. They haya 
involved a combination of gravimetric or gasometric methods, always 
with a final estimation of the ethylene oxide alone by titrimetric 
procedure. The method here developed involves only direct gasp. 
metric measurements made with a standard burette-manometer-com. 
pensator unit. Its accuracy is commensurate with that obtained by 
the best of the former methods, and both components of the mixtyre 
may be determined in less time and with less effort than were required 
previously for the determinations of ethylene oxide alone. 


II. PREVIOUS INVESTIGATIONS 


Deckert ' employed the reaction of ethylene oxide with hydrochloric 
acid and sodium chloride to form glycol chlorohydrin, titrating the 
excess acid with sodium hydroxide to determine ethylene oxide. 
Later, Deckert? reported a colorimetric method based upon the 
reaction of ethylene oxide with a concentrated solution of sodium 
chioride to form glycol chlorohydrin, phenolphthalein and bromothy- 
mol blue serving as indicators. He reported a third method (colori- 
metric) * for the examination of air containing traces of ethylene oxide, 
based on the reaction with potassium thiocyanate, with phenol- 
phthalein as indicator. 

According to Lubatti,* Deckert’s titrimetric method indicated only 
90 percent of the ethylene oxide present. Lubatti employed a nearly 
saturated acid solution of magnesium chloride as an analytical reagent 
for the examination of samples of liquid ethylene oxide, titrating with 
sodium hydroxide. The reaction was over 99 percent complete. 

More recently the apparatus described by Page ® was utilized by 
Lubatti,® who modified his original modification of Deckert’s method 
to the extent of adding sodium, lithium, and magnesium chlorides to 
0.1 N hydrochloric and sulfuric acids, titrating the excess acids with 
sodium hydroxide. 


III. DEVELOPMENT OF A GASOMETRIC METHOD 


1. PRELIMINARY CONSIDERATIONS 


The reduction in gas volume which occurs when ethylene oxide 
reacts with reagents of these types was suggested by Deckert ‘ as an 
obvious means of determining ethylene oxide. The gasometric pro- 
cedure has been employed in several laboratories, but the results 
have not encouraged publication. The reagents used to absorb 
ethylene oxide have been aqueous solutions of various salts or mix- 
tures of these prepared indiscriminately. It is not unlikely that these 
reagents reacted quantitatively with ethylene oxide, since water 
alone will do so, and with great rapidity if the gas is bubbled througli 
it.’ The major mistake apparent in these efforts has been the use 

1 W. Deckert, Z. anal. Chem. 82, 297-307 (1930). 
1W. Deckert, Z. angew. Chem. 45, 559-62 (1932). 
3W. Deckert, Z. angew. Chem. 45, 758 (1932). 
40. F,. Lubatti, J. Soc. Chem. Ind, 51, 361 (1932). 
8A. B. P. Page, J. Soc. Chem, Ind. §4, 421 (1935). 
¢ 0. F. Lubatti, J: Soc. Chem: Ind. 54, 424 (1935). 
T See footnote 1. 


* Several experimenters state that the reaction is slow. This is not true if an efficient absorbing pipette 
is employed. 
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of 100- to 200-ml portions of aqueous reagents in ordinary forms of 
absorption pipettes. — Thus speed of reaction has been attained, but 
with a grievous sacrifice of accuracy, since carbon dioxide dissolves 
abundantly in such aqueous solutions. Elaborate precautions have 
been taken to saturate the reagent (at controlled temperature) with 
carbon dioxide. This prerequisite to actual analysis not only re- 
guires much time, but unfortunately does not accomplish its objective 
ance the dissolved carbon dioxide wanders in two directions: Into the 
ir through the atmospheric arm of the pipette, and into the gas 
smple during the absorption of ethylene oxide. The amount lost 
by the reagent to the gas sample may be tediously replaced by re- 
peated passage into the pipette, but the amount lost to the air must 
pe replaced by an equal amount from the gas sample itself, thus 
causing a proportionately high estimation of ethylene oxide. 

4 brief outline of a successful gasometric method for the analysis 
of a mixture of ethylene oxide and carbon dioxide should serve to 
make this problem clear. 

To begin with, the sample taken for analysis must be measured in a 
iry burette. If the mojst burette generally used were employed, 
water present would react with the ethylene oxide (hereinafter des- 
mated C,H,O) and alter the composition of the sample before its 
analysis. 

Next, the dry sample must be transferred to a reagent which will 

quantitatively remove C,H,O, without removing CO). 
“Since the reagent which removes C,H,O will also add H,O vapor to 
the residual gas, this gas must be returned through a suitable desic- 
cant in order that it may reach the burette in a dry state. The 
liference between the measured dry sample and dry residue should 
then represent the C,H,O in the sample. 

Finally, the CO, must be removed, and the residual gas returned 
after drying) to the burette for measurement. The contraction 
observed will then represent the CO, in the sample. 

Since the sample may consist entirely of C,H,O and CO,, a measured 
volume of dried air must be added during the initial steps of the 
analysis. This will act as a transfer or ‘‘push” gas, and must be inert 
with respect to the reagents employed. The difference between its 
volume and the volume of the residual gas observed at the end of 
the analysis represents the inert content of the sample. 

From this outline, it can be seen that the development of a successful 
gasometric method for the analysis of this mixture will resolve na- 
turally into three separate problems, as follows: 

1. Selection of a reagent which is chemically and physically selective 
with respect to C,H,O in the presence of CO, and air. 

2. Selection of a desiccant which is chemically and physically pas- 
sive with respect to CO, and air. 

Since a satisfactory reagent for CO, (in the presence of air) already 
exists, the final problem becomes: 

3. Design of a volumetric apparatus capable of performing the 
lecessary procedures indicated. 

These problems will be discussed in their order. 

11827339 4 
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2. REAGENT FOR ETHYLENE OXIDE 


A suitable reagent for the determination of C,H,O must possess i, 
following qualifications: ” 

1, The reagent must react quantitatively with C,H,O, completely 
removing this substance from the gas phase, or replacing it ey 
stoichiometrically definite volume not equal to its own. 7 

2. The product or products of reaction must possess negligible yayo: 
pressures (or vapors which can be removed chemically without affect. 
ing CO, or air). 2 

3. The reaction must produce no gaseous by-product (unless this 
can be removed). 

4. The reagent or products of reaction must not react with, or dis. 
solve, significant amounts of CO, or air. . 

5. Since any reagent likely to be conceived will dissolve CO, to 
some extent, the volume of the reagent actually employed must be 
very small in order to reduce the attendant error to the order of 
magnitude of the required volumetric precision. 

6. The above qualification at once imposes several additional and 
severe ones. Since the volume of the reagent must be reduced to 
0.2 to 0.3 ml, it can most conveniently be confined over mercury. 
Therefore, it must not react with mercury to produce significant 
amounts of gas. 

7. A pipette of the type mentioned above is at best an ineflicien: 
absorber, presenting little contact between gas and reagent. Ae. 
cordingly, the rate of reaction between the reagent and C,H,O must 
be very rapid. 

8. Finally, because of its small volume, and for the sake of con- 
venience, the reagent should have an extremely large capacity for 
C,H,O, so that its frequent replacement may be avoided. 

This list of qualifications would appear to be somewhat appalling 
to the gas analyst; but actually they were easily realized. The 
behavior of C,H,O with various substances was roughly studied ina 
preliminary set of experiments. C,H,O was observed to react readily 
with the contents of the reagent bottles at hand on the laboratory 
shelf, such as HNO;, HCl, H;PO,, H,SO,, CrO3, H,Cr,0;, HgNO,, etc, 
in various dilutions with water. As a matter of fact, since C,H,0 
reacts readily with H,O alone, it seems likely that few inorganic 
reagents would prevent the formation of the glycol in the presence 
of water, although some accelerate the reaction. The apparent 
tendency of C,H,O to react with many such mixtures suggests the 
possibility of a number of satisfactory reagents. 

Only one from this list of reagents was selected for extensive study. 
The reagent selected was a mixture of equal parts of H,O and con- 
centrated H,SO,. This met every one of the qualifications previously 
delineated. In the event other reagents develop, it would not seem 
worth-while to issue additional papers for each in turn. 

Apparently, the reaction between the diluted H,SO, andC,H0 
goes through an intermediate state to form ethylsulfuric acid, which 
breaks down to H,SOQ, and the first glycol. This in turn reacts with 
additional C,H,O to produce higher glycols, increasing enormously 
the capacity of the reagent. 

The behavior of this reagent with C,H,O is interesting. If the gas 
is confined over mercury within a glass tube, and the reagent is adult 
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Figure 1.—Reaction unit. 
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‘ed in successive small portions of 0.2 ml each, the initial reaction is 
very rapid and the final reaction very slow with the addition of each 
portion of reagent. However, if the wall of the tube is coated with 
‘he reagent by shaking, the reaction thereafter is extremely rapid. 
The glycols which form drain slowly, producing a film or droplets of 
honey-like consistency over the entire wall of the tube. 

Over 2,500 ml of C,H,O was passed (at a rate exceeding 200 ml per 
minute), into a tube containing 0.5 ml of the H,SO, solution. At the 
ond of this experiment, the reaction was still proceeding rapidly and 
completely. Approximately 2 ml of the glycols had been collected. 
The capacity of the reagent Is accordingly over 1:5000, expressed as a 
ratio by volume of reagent:C,H,O, a novel behavior from the view- 
point of the gas analyst. 

' There was no direct evidence of the formation of interfering by- 
products of the reaction. Acetaldehyde was suspected, but only 
‘races of it were found. 


3. SUITABLE DESICCANTS 


To qualify for occupancy of the drying tube within the analytical 
apparatus, a desiccant must possess the following characteristics: 

|. It must leave somewhat less water vapor in the dried gas than 
corresponds to the volumetric accuracy of the apparatus. This limit 
with respect to the desiccant can be set at 0.02 ml. 

2. It should dry the gas rapidly, preferably at the rate of 100 ml per 
30 seconds or less. 

3. It must not react with CO, or air, nor adsorb significant quanti- 
ties of these gases. : 

4. It should not react objectionably with C,H,O in the event that 
this gas is accidently brought into contact with it through a mistake in 
operation. 

5. It should be convenient to use. 

Of the available desiccants initially suggesting themselves, P.Os, 
Meg(ClO,)2, CaSO,, and CaCl, (fused, granular, anhydrous) seemed best 
suited. These were studied experimentally, and the following facts 
were ascertained: 

1, All four met the first qualification, CaCl, in its best available form 
just passing. 

2. All four met the second qualification. 

3. POs, Mg(ClO,)2 and CaCl, met the third qualification. (CaCl, 
must, of course, be of reagent quality. Traces of alkaline impurities 
ust be neutralized by treatment with CO,-+air before the gas 
aualysis.) CaSO, adsorbed large amounts of CO, and was therefore 
discarded. 

4. Of the remaining three desiccants, CaCl, alone met the fourth 
qualification completely. This is a “should” and not a ‘‘must” quali- 
fication, however, since C,H,O would not come into contact with the 
desiccant except by mistake. 

If C,H,O is inadvertently passed over POs, the reaction is marked, 
as is the P,O; with a visible deposit of carbon, which leaves no doubt 
as to the necessity of its replacement. 

Mg(ClO,). is open to a curious objection. A double mistake might 
toncelvably result in (a) bringing H,SO, into contact with this desic- 
rant, and later (b) bringing C,H,O into the same tube. It is not 
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probable that such mistakes would }, 
made; but if they were, a serious ba 
plosion might occur. Meg(ClO,), was 
actually used in the analytical work 
which follows. Since it demonstrate) 
no marked superiority, it had probably 
best be laid aside in this instance. * 

5. Of the two remaining desiccants 
CaCl, is by far more convenient to use. 
However, P,0; may, without too great 
difficulty, be sifted onto glass woo! 
within the drying tube of the apparatus 
This is best done by plugging a tube 
(smaller than the drying tube) at one 
end with glass wool, drawing (by 
vacuum) P.O, from the reagent bottle 
into this tube, and finally inserting the 
small tube into the desiccating tube of 
the apparatus, discharging its contents 
therein by means of a thin glass rod. 
Once this tube is filled, prolonged 
service is to be expected. 

The analyst will probably make his 
own choice between PO; and CaC); 
With respect to H,O and CO,, P.O, 
offers the best chemical behavior. Both 
desiccants gave good analytical results. 
In any event, the same desiccant 
should be used in both drying tubes 
attached to the apparatus. 


4. DESIGN OF THE VOLUMETRIC 
APPARATUS 


In measuring the mixture of C,H,0+ 
CO, gasometrically, it is obvious that 
the burette must be dry. The dry 
burette accordingly should be com- 
plemented with a dry compensator. 
The burette and compensator-manom- 
eter unit of our standard gas-analysis 
apparatus ® served satisfactorily. 

Aside from measurement, it is neces- 
sary only to transfer the measured 
dry sample to the reagent which re- 
moves C,H,O and subsequently remove 
the water vapor which this reagent 
adds to the residual gas prior to return- 
ing it to the burette. Thereafter, the 
residue is to be passed into a reagenl 
which removes CO;, and the residual 
gas dried and measured. 


* M. Shepherd, BS J. Research 6, 122-180 (1931) RP266. 
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An apparatus ’° by means of which these procedures may be easily 

achieved is shown in figure 1. 
" This unit is sealed directly onto the horizontal outlet of stopcock 4 
of the burette-manometer-compensator assembly previously men- 
‘ioned, which is reproduced here for convenience (fig. 2). A thermo- 
plastic cement (pine tar and shellac) serves to make the joint at the 
Jeeve, D, (fig. 1). Aceording to tests made, rubber reacted with 
(,H,0 too rapidly to permit its use. The cement reacted too slowly 
to cause significant error. Actually, the main concern is at the 
sampling inlet. : ; 

An ordinary bubbling pipette connected at the outlet, A, of T-stop- 
cock 1 is filled with a concentrated solution of KOH to remove the COs. 
A concentrated solution is used because of the low solubility of 
atmospheric inerts. )" 

A reservoir, B, connected to T-stopeock 2 serves as the reaction 
pipette for the C,H,O. It is filled with mercury, and a leveling bulb 
is attached to its lower outlet. The small volume of reagent (0.2 to 
).3 ml of equal parts by volume of H.O and concentrated H,SO,) 
may be conveniently introduced into B by means of a medicine 
dropper drawn out into a tip which is bent downward at a right angle. 
This tip may be inserted in the capillary inlet of B when the key of 
cock 2 is removed. One such charge will serve for a hundred or 
more analyses before it must be removed in the same manner. 

The capillary inlet of B has a reference or zero mark etched around 
the stem just under cock 2. The lower outlet of B connects to a stop- 
cock by means of nitrometer tubing, around which is placed a com- 
pression screw for adjusting the reagent meniscus to this zero mark. 

The tube, C, contains a desiccant for drying the gases returned from 
Aor B. Stopeock 3 is drilled with three bores, as indicated on the 
handle. When turned 45° counterclockwise to the position shown in 
figure 1, tube Cis by-passed, and gas from the burette flows directly 
within the horizontal distributor to either _A or B. When the cock is 
turned in the position shown, gas flowing from A or B to the burette 
must pass through C, entering at the bottom of this tube and leaving 
the top, as indicated by the directional arrows. 

The tube, C, is closed by means of plates, H—-E, which are ground flat 
to fit corresponding ground seats at the ends of C. A heavy lubricant 
serves to seal them. They are held in position by two coil springs 
fastened to saddles of sheet metal, which, in turn, fit over the glass 
knobs of the end plates, H-E. Thus, the tube, C, may be easily opened 
in position on the apparatus for supplying or removing desiccants. 
Plugs of glass wool at the ends of C keep the desiccants in place. 

The gas volume of C, when filled with desiccant, is small enough so 
that changes of 0.2° C. in room temperature will produce no significant 
error in comparing volumes during an analysis. The time required 
for an analysis is so short that trouble from changing temperature 
should never occur.!? 


* This apparatus was designed by Shepherd. The unit comprising cock 3 and tube C is to be used in a 
multip e assembly in a new type of gas-analysis apparatus now under test. 
"J. R. Branham and M. Sucher, J. Research NBS 21, 63 (1938) RP1113. 
_ '! One of the most embarrassing questions asked the authors with respect to any type of analysis is, ‘How 
2hg Coes it take to do it?” There are so many factors involved—the accuracy desired, the composition of 
ture, and, ofcourse, the analyst—that definite answerscannot begiven. Butinthe present instance, 
y be of interest to know that it is possible to complete an analysis in 15 minutes, operating carefully. 
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To complete the apparatus, a drying tube is connected to the lef. 
hand outlet of stopcock 33 of the burette (fig. 2), and the source of 
sample to the right-hand outlet of the same cock. These connections 
were made with sleeve joints and thermoplastic cement in the presen; 
work, but ground joints would be preferable. 

The apparatus described in figure 1, as previously noted, may bp 
used in conjunction with the regular volumetric unit and with the 
supporting frame of our standard gas analysis apparatus. If, how. 
ever, the number of samples to be analyzed would justify a separate 
apparatus for the purpose, the following modifications would improve 
the unit somewhat: 

1. The burette should be blown with a 20-ml bulb at the top, with 
intervals on the stem from approximately 20 to 100 ml. This would 
(a) shorten the apparatus or (b) permit estimations of intervals of 
smaller volume if the apparatus were not shortened. 

2. The two outlets of stopcock 33 of the burette (fig. 2) should 
terminate in interchangeable ground-glass joints (male members), 
to permit easy connection to an auxiliary drying tube and source of 
sample. 

3. The sleeve, D, (fig. 1), should be tapered (female) to fit a cor- 
responding taper (male) on the outlet of cock 4 (fig. 2). 


5. VOLUMETRIC BEHAVIOR OF THE DRY APPARATUS 


Before developing the analytical procedure, the volumetric be- 
havior of the dry burette and compensator was studied experi- 
mentally. The following facts were ascertained: 

1. With a desiccant in the tube, C, and with the reservoir, B, the 
burette and the compensator dried, a volume of dried air could be 
repeatedly measured, transferred from the burette to the reservoir 
and back, and even allowed to stand overnight, with no significant 
change in volume (+0.02 ml). Indeed, there are indications that 
the dry system might prove more accurate than the saturated (wet) 
system ordinarily employed. 

2. With water in the reservoir, B, dry air transferred to B and re- 
turned to the burette through C (containing any desiccant later 
employed) gained approximately 0.05 to 0.08 ml if the rate of return 
exceeded 70 ml per 20 seconds. However, no change of volume was 
noted at the lesser rates of gas flow employed. Returning a gas 
from B through C in 30 seconds insures drying. 

3. If dried air is passed directly over water in B and returned to 
the burette, by-passing C, the volume is increased to correspond to 
the existing vapor pressure of water. Thereafter, if the moist air \s 
passed back and forth between B and the burette through the desiccant, 
the volume gradually returns to its original dry value. Approximately 
4 or 5 passages (2 minutes each) through the desiccant are required 
to restore the dry condition, since adsorbed water is apparently given 
up reluctantly by the burette wall. This illustrates the importance 
of always returning gas to the burette through the tube, C. 


6. ANALYTICAL PROCEDURE 


The following detailed procedure was developed during the course 
of a few preliminary analyses, and of the analyses reported in table |. 
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Initial drying of apparatus.—The compensator was dried by alternate evacuation 
and filling with air dried over P;Os. (A small amount of desiccant left permanently 
in the bottom of the compensator would probably be desirable in an apparatus to 
be used for this type of analysis only.) The burette was dried by repeatedly 
passing air in and out through the desiccant connected to cock 33 (fig. 2). The 
distributor was dried by passing air through C to A and back. About 10 passages 
of 2 minutes each are required for drying, if no water was originally visible. 
Onee dry, the apparatus should remain so. 

Analytical technique.—The distributor, manometer, and C are filled with dried 
air and adjusted to the fixed pressure of the compensator. 

Approximately 60 ml of the C,H,O+ CO, mixture to be analyzed is drawn into 
the burette. With cock 22 closed the pressure is balanced against atmospheric 
with the leveling bulb, so that a minimum amount of gas will later flow between 
the burette and manometer. Cock 4 (fig. 2) is then turned to the position - to 
connect the manometer with the burette; the pressure is balanced against the 
compensator, cock 4 is turned to the position y to disconnect the burette from the 


manometer and train, and the volume of the sample is read. Cock 8 (fig. 1) is 
turned to the position X, the leveling bulb attached to B is lowered slightly to 


reduce the pressure therein, cock 2 to the position Fk, and cock 4 to the position L. 
The sample is transferred to B, the leveling bulb being dropped to maintain 
approximately atmospheric pressure in B. When the mercury in the burette 
reaches cock 22, it is turned to connect to cock 33 and the line between them is 
flushed with mereury. Cock 33 is then turned to connect the burette to the source 
of dried air, 25 to 30 ml of which is taken into the burette and measured exactly 
as was the sample. The gas in the arm of the manometer connecting to the 
distributor is then displaced to the burette, cock 4 is turned to the position y¥, the 
pressure of the gas in the burette adjusted to slightly in excess of atmospheric by 
the leveling bulb, and the air is then transferred to B with cock 4 in the position L. 
The dry air has now swept the C,H,O+ CO, into B, where the mixture is allowed 
to stand 1 minute. Agitation by means of mercury flow is not necessary after 
the first passage of the first analysis. Cock 3 is now turned to the position 
and the residual gas (minus C,H,O) is returned through C to the burette, 30 
seconds being taken for passage through C. The volume is measured as before, 
the contraction noted being computed as C,H,O. 

The gas in the distributor arm of the manometer is again displaced to the burette, 
and thence transferred to the pipette, A, containing a solution of KOH, being 
routed through C, both back and forth. During the eighth and tenth passages 


into A, cock 3 is turned for a second to the position 4 to displace CO, trapped in the 


straight bore of this cock. At no time must gas by-pass C on the way back to the 
burette. After the tenth passage, the residual gas (minus CQ.) is measured and 
the contraction computed as CQ:. (Ten passages are used to insure complete 
removal of CO, from the tube, EH, and the distributor, by dilution.) The difference 
in volume between this residue and the dry air originally taken for a transfer gas 
is computed as the inert content of the sample. 


IV. ANALYSES OF MIXTURES OF ETHYLENE OXIDE AND 
CARBON DIOXIDE 


1. PREPARATION OF MIXTURES FOR ANALYSIS 


To test the method, apparatus, and procedure just discussed, two 
mixtures of C,H,O and CO, were prepared and analyzed. Both of 
these mixtures were made to contain approximately 10 percent of 
C;H,0 by volume, which corresponds to the mixture used for fumi- 
gating. ‘The mixtures were prepared in a 5-liter glass flask to which a 
barometric manometer was attached. The flask was evacuated, 
C,H,O was admitted until the manometer registered a pressure of 
approximately 75 mm, and CO, was then admitted until the pressure 
reached approximately 750 mm. The room temperature was ob- 
served to remain constant (+0.1° C). The two gases, of almost 
identical molecular weights, were allowed to mix for several days. 
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Convection was occasionally induced by local warming of the 4 
The flask was then joined to the burette as previously described. an: 
samples were withdrawn using the mercury within the burette as . 
pump. ; 

The composition of the samples prepared‘in' this manner was yo; 
necessarily expected to check the composition as determined by thy 
analysis, even assuming absence of significant errors in both Cases 
Since the pressure-volume relationships existing during the prepurg. 
tion and during the analysis of these mixtures are different, the devia. 
tions of the gases and their mixtures from the behavior prescribed {y; 
ideal gases could account for differences in the composition as measured 
in the two cases. 


] 
ask, 


a 


f2. PURITY OFiTHE'ETHYLENE OXIDE 


The C,H,O used in preparing the two mixtures was purified jin , 
large commercial rectifying column. Approximately 10 liters of this 
was stored as a liquid in a cylinder. About one-third of the contents 
of this cylinder had been withdrawn before the present use of the gas 

The completeness of reaction between this C,H,O and the H,0. 
H,SO, reagent employed during the analyses was determined as 
follows: 

1. 30.08 ml of dried air was measured in the burette and transferred 
to the reservoir, B, containing the reagent. 

2. 274 ml of C,H,O was measured in three portions and transferred 
to B. (The reaction was so rapid that only with great difficulty 
could the mercury in B be prevented from being drawn into the 
horizontal distributor.) 43 4% « 

3. The residual volume was passed through CaCl, in C, and meas- 
ured in the burette as 30.53 ml. 

4. The difference between the initial and final volumes corresponds 
to 0.017 ml of “‘inert’”’ per 10 ml of C,H,O. Since 6 ml of C,H,0 is 
the amount in the average sample, approximately 0.01 ml of the 
C,H,O used did not react or was inert during each analysis. This is 
less than the expected volumetric error. 

The above test does not preclude the possibility that a compound 
other than C,H,O was present and reacted with the solution of 
H.SO,. The purification of the C,H,O by rectification makes such a 
possibility remote; and analyses of C,H,O+CO, mixtures obtained 
from the same source as that supplying the present C,H,O, made by 
isothermal fractional distillation in our laboratory, failed to disclose 
a notable impurity in the C,H,O fraction. 


3. PURITY OF THE CARBON DIOXIDE 


The carbon dioxide used was drawn from the last half of a cylinder 
of liquid obtained from a commercial source. Its “purity,’”’ so far 
as this work is concerned, was tested as follows: 

1. Four volumes totalling 384 ml of CO, were measured; 30.43 ml 
of air was measured separately. 

2. These gases were transferred to the pipette containing a saturated 
solution of KOH, first without bubbling, and finally with repeated 
bubbling passages. 

3. The residual air was dried in C and measured as 30.54 ml. 
This corresponds to 0.03 ml of “‘inert’’ per 100 ml of COs, or about 
0.016 ml per analysis. (A small amount of this inert may have beer 
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air displaced from the solution of KOH during the absorption of CO).) 
This is of the same order of magnitude as the volumetric accuracy at 
its best. 

4. ANALYSIS OF MIXTURE 1 


According to the pressures observed in preparing this mixture, its 
composition was calculated to be: 
(%, by volume; or mole %, 
assuming ideality) 
of _—_—_———— _ 9.59 +0.03 
CO; 90.15 +.03 
i 0.26 +.03 
Correcting these for the amount of inert found in the C,H,O and 
CO;: 
. 9.57 +0.03 


90.12 +.03 
0.31 +.03 


TABLE 1.—Analyses of mizture 1 





| 

Percentage, by volume 
: A=deviation from average 
Analysis 


Yasiccan 
Desiccant number 





C:H,O A CO, A Inert 





9. 91 0.05 | 89.91 0.12 0.19 
9. 81 -05 | 89.98 05 . 22 
9. 84 -02] 90.14 | . 02 
9. 87 .O1 90. 07 . 04 . 06 


9.86 | +0.03 90.03 | +0. 08 0. 12 


9.66 | 0.04 | 90.08 |” 0.08 
. 58 04 | 90.12 . 04 
. 63 01] 90.16 . 00 
61 -O1 | 90,27 I 
.62 | +0.03 90. 16 _£0. 06 

0.04 | 90.17 | 0.02 7 
02 | 90. 20 . 05 
01 | 90.09 . 06 
-05 | 90.15 .00 

+0.03 | 90.15 | +0. 03 | 

0.06 | 89.99 | 0.17 
. 08 01 
03 | 90, . 05 
eT M4 . 23 

+0.07 | 90.16 | +0. 12 a 


63 | +0.05 | 90.16 | £0. 07 | 
‘aa | 
6 H ) | ) | (*) 


a 











| Mg(C1O4)3-_--- 
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_ The analyses of this mixture are given in table 1. Determinations 
) to 16, inclusive, made with 0.3 ml of reagent and with three different 
desiccants, give the following average composition: 


(%, by volume, as 
observed) 


9.63 +0. 05 

90.16 +.07 

0.22 +.07 

This agreement is closer than was to be expected. 
A study of the data noted in table 1 discloses the following facts: 
|. The value for C,H,O is apparently 0.2 percent too high, and that 
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for CO, correspondingly too low, if 2 ml of the H,SO, reagent is use, 
instead of 0.3 ml. (Compare the average of determinations 1 to 4 
with the average of 5 to 8.) This illustrates the effect of the solubility 
of CO, in the aqueous reagent, and the importance of reducing the 
volume of the reagent to a working minimum. 

2. The results obtained when the desiccants Mg (ClO,):, Ca(), 
and P.O; are used are concordant, since no significant differences 
occur between the averages for the determinations 5 to 8, 9 to 19 
and 13 to 16, inclusive. When CaSO, is employed, the results ay: 
aberrant. (Note determinations 17 to 19.) 

3. The reproducibility attained is adequate. 


5. ANALYSIS OF MIXTURE 2 


The analysis of this mixture was a slightly more extended effor 
than that of the first one. According to the pressures observed when 
preparing the mixture, the composition (corrected for inert found jn 
C,H,O and CO,) was as follows: 

%, by volume 

9.97 +0. 03 
89.98 + .03 

0.05 + .03 

The average of 28 determinations reported in table 2 shows the 
composition of this mixture to have been: 

%, by volume 
+0. 04 
+ .05 
+ .04 


TABLE 2.—Analuses of mixture 2 








Percentage, by volume 

: Analy- A=deviation from average 
Desiccant sis 
number 





C2H,O CO: Inert 








9.96 0.03 89. R4 0. 20 
9.86 07 90, 01 : mL 
9. 82 .11 | -90.08 09 
9.93 . 00 &9. 99 . 0S 
10. 08 15 89. 92 . 00 


9.91 .02 | 90.07 . 02 
9. 83 -10 | 90.02 5 15 
9.98 .05 | 89.99 04 
9. 98 .05 | 89.88 : 13 
9. 91 .02 | 90.02 0 07 


9. 93 -00 | 90.05 
9. 95 .02 | 90.00 





Average biel tial 9.93 | 40.05 | 89.99 .08 | +0. 05 














9. 93 0.01 90. 01 0.03 | 100.00 
9. 99 . 07 89.95 : . . 03 oa : 
} , . 90 . 02 90. 06 B 01 W).( 
RFR RIOD osiks acaminacavevecnveueds ce 00 90, 05 a 05 00 100.0 
9. 85 . 07 90. 15 ‘ ; . 03 100 W 
9.95 . 03 90. 05 P ‘ .03 | 100.0 


Average oetenta 9.92 | +0.03 90. 05 +0. 02 | 


——E 
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TABLE 2.—Analyses of mixture 2—Continued 





Percentage, by volume 
Analy- A=deviation from average 
sis 
number 





C2H4O A CO; Inert Total 


10. 0. 07 89. 99 . 06 0. 00 100, 00 
.05 90, O1 ¢ . 00 st 100. 00 
Ol 90. 05 ‘ . 02 . 00 100. 00 
. 03 90. 03 ‘ . 00 . 02 100. 00 
.02 90. 04 ‘ ‘ . 02 100. 00 


03 90, 07 ‘ ° . 00 100. 00 
«i} 90. 17 oi . . 02 100. 00 
. 04 90. 10 OF . . 02 100. 00 
Ol 90. 03 ‘ . 04 .C2 100. 00 
01 89. 97 ‘ . OF . 06 100. 00 





0. 04 90.05 | +0. y +0. 02 


























| 9.93 | +0. 04 | 90.02 | +0. 05 





The agreement in this case is about as close as before. 

For the convenience of study, the data in table 2 include the devia- 
tion of each analysis from the average, the average and the average 
deviation of each of the three groups of analyses, and the total aver- 
aves and deviations for the entire series. A study of these data gives 
the following information: 

1. So far as the C,H,O is concerned, P,O;, CaCl,, and a mixture of 
the two serving as the desiccant are equally good. 

2. The CO, tends to be slightly lower when CaCl, alone is employed. 
However, the difference between the amount obtained when using 
this desiccant, and the amounts obtained when using P.O; or P,O;+ 
CaCl, is a matter of no concern. 

3. The results on the whole are concordant, and the reproducibility 
attained should satisfy the general need (+0.05 percent). The 
maximum deviation noted was 0.15 percent, which should serve to 
indicate the desirability of check determinations. 

It may be noted that the remarkable totals obtained (last column) 
have little significance other than to testify to the fact that no meas- 
urable amount of air leaked out of the apparatus during the analysis 


V. ACCURACY OF THE METHOD 


1. PRESSURE-VOLUME CONSIDERATIONS 


So far we have made mention of reproducibility only. This is as 
faras the gas analyst has usually found it convenient to go. Repro- 
ducibility, often under the pseudonym of precision, has long mas- 
queraded under the cloak of accuracy. It is a demonstrable entity, 
and can be captured by anyone with suflicient patience. Accuracy is 
not so easily demonstrated. This is unfortunately true in the present 
Instance, 

If two gases of known and satisfactory purity, which deviate by no 
significant amount from the formulas laid down for ideal gases either 
separately or mixed, had been mixed according to the pressure method 
employed in preparing the samples previously analyzed, then, indeed, 
4 definite figure for the accuracy of this method could be expressed. 





184 Journal of Research of the National Bureau of Standards (Volt 


But these two gases, C,H,O and CO,, are known to deviate from 
ideality as separate gases, and the extent of the deviation of C.H0 
has not been measured. Worse yet, their behavior when together 
and when mixed with air, has not been reported. It seems unreasop, 
able to suppose that these two, noticeably imperfect with respect to 
Boyle’s law as indiwiduals, would strictly observe only one of the 
perfect gas laws when brought together. Thus it seems scarcely gp. 
sistent to assume deviations from Boyle’s law and ignore deviations 
from the unnamed laws of additive volumes—an assumption that hy: 
already served as the thesis in the development of some impressiye 
formulas for correcting gasometric analyses.” 

{t will be worth while to show experimentally the behavior of thoy 
two gases when mixed, in order to illustrate the danger of makins 
the above assumption. : 

C,H,0 will vary from Boyle’s law considerably more than will CO 
For this reason, the samples were prepared by measuring C110 x 
a low pressure, in which state it could be expected to more nearly 
behave itself. Even so, it would not be safe to compute the grayi- 
metric composition of this mixture from such pressure data. Re- 
flection over the following data will serve to support this feeling jy 
the matter. 

Suppose 40 ml of dry air is measured in the burette of the apparatus 
used for these analyses. If it is then stored, and 40 mi of dry C0, 
is measured at the same temperature and pressure—and if, now, the 
40 ml of air is returned to mix with the 40 ml of CO., what would the 
volume be? This experiment was performed, and the resuliing 
volume (average of 10 determinations) was 80.06 (+0.02) ml. Next, 
40 ml of air and 40 ml of C,H,O, measured separately as before, wer 
mixed; the resulting volume (average of 10 determinations) was 80.37 
(+0.04) ml. Now, 40 ml of CO, and 40 ml of C,H,O were measured 
separately, and mixed. If a law of ideal mixing prevailed, the expected 
volume would be 80.00 ml. On the other hand, it might be inferred 
from the above data that the combined volumes of these two gases 
would be either 80.43 or 80.21, depending on the particular assumption 
invoked. But the average of 10 determinations was 80.13 (+0.01) ml. 
(40 ml of air+40 ml of air were measured as above to check burette 
calibration and reproducibility.) It is true these measurements are 
not sufficiently accurate to justify formulation of equations defining 
the exact behavior on mixing; but the order of magnitude of the 
differences is certainly quite significant. The results are convincing 
enough evidence that the volume relationships on mixing in this gas 
analysis (and in others) should be investigated. Without such 
investigation, it is impossible to say whether or not such behavior 
on the part of gas mixtures will cause significant error. i. 

To get a clearer picture of the pressure-volume relationships existilg 
during this analysis, consider the volumes actually observed. These 
are (average values): 


13 Dalton’s law, apparently an extension of the law of Henry, has been interpreted in modern textbooks 
as meaning that the total pressure of a mixture of gases is equal to the sum of the partial pressures of the 
components. The gases are, therefore, still thought of individually, and one is considered to have no physica 
effect on the others. This concept has been extended in some texts, and in the minds of many experimenters, 
to be the practical equivalent of the statement that the total volume of a mixture of gases is equal to et 
sum of the volumes of the components measured separately, and, of course, at the same temperature a0" 
pressure. The latter statement must necessarily involve Boyle’s law, and, in the absence of a foruaal 
sponsor, it is difficult to name any law of additive volumes. The point so far not taken into account by gas 
analysts is that mixed gases may not behave as individuals. Thus, if two gases are each measured seps- 
rately at the same temperature and pressure, the sum of their volumes may not be equal to the total volume 
of the mixed paves Ceasorused at the same temperature and pressure); and the deviation from this sum ™s) 
be significantly different from the deviation that might be expected from the law of Boyle alone. 
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1. The sample is composed of 54 ml of CO,.+6 ml of C,H,O=V,. 
9 24 ml of air, measured separately from V,, and used as an inert 

transfer gas=} ——_ . ; ; 

3 54 ml of CO,+24 ml of air, the residue after absorption of 

C,H,O= I 3° ° . ° ‘ a 
4. 24 ml of air, the residue after absorption of CO,= V,. 

The analysis is computed: 


Vi ha Ys _.0,H.O; a SS CO). 
V 1 \ 1 


Of these volumes, V, and V, are definite; but V; and V; must be 
evaluated, since they cannot be arrived at by simple arithmetic. 

In order to evaluate V, and V3, the following measurements were 
made: 

|. 54 ml of CO, and 6 ml of C,H,O were measured separately and 
when mixed, both gases being dry and at the same condition of 
temperature and pressure (approx. 1 atm, 25° C). The average of 
\) determinations gave this result: 


54 ml of CO,+6 ml of C,H,O=60.07 +0.03 ml= Vj, above. 


2. The above measurement was verified by 10 determinations of 54 
ml of air+-6 ml of air, which served to check the burette calibration 
aid manipulation. 
3. 54 ml of CO,+-24 ml of air, measured separately and mixed as 
ox, ggbelore, gave as the average of 10 determinations: 
vere 
37 54 ml of CG,-+-24 ml of air=78.04 +0.02 ml= V3, above. 
red 
eted 
rred 
“ses 
tion 
mil, 
‘ette 


4. The above measurement was verified by 10 determinations of 
j4 ml of air-+24 ml of air, to check the burette calibration and the 
manipulation. 

5. The measurements made above were further authenticated by 
passing measured volumes of C,H,O and of CO, from the burette to 
the storage reservoir and back, without significant gain or loss of 
pone volume. ’ a 
sae These evaluations of V,; and V3; at once disclose deviations whose 
the Maoder of magnitude (0.1 percent) could produce significant error. 
cing (Out fortunately in this particular instance the computation partially 
oys Aactncels the effect of deviation from ideality. Thus, the deviations, 
wich Mae taken into account, would increase the C,H,O by only 0.05 percent 
vior Maeed the CO, by 0.07 percent. 


ting 2. GRAVIMETRIC DETERMINATION OF ACCURACY 
While the above measurements give assurance that the accuracy 
of the method probably is commensurate with the reproducibility 
obtained, a definite experimental check is in order. The accuracy was, 
therefore, determined by preparing several mixtures of C,H,O+CO, 
tanging In composition from about 7 to 12 percent of C,H,O. These 
mixtures were made up as before by measurement of pressure at 
ves ae Xed volume and temperature, but at the same time the composition 
‘os Mas determined gravimetrically, and by the analytical method here 


olume 


n may reported. 
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(a) APPARATUS 


The apparatus employed is shown in figure 3. It comprise 
5-liter flask, C, to which is attached a barometric manometer B ° ‘ 
small brass cylinder, A, fitted with a needle valve of the vacuum Sra 
connects to C through an interchangeable ground joint. An aurifi : 
flask, E, may be attached at the same joint and serves as a rss 
pipette. The flask and connections may be evacuated by a ate 
through cock 1, and supplied with C,H,O and CO, through the ae 
inlets of cock 4. ° ' 



































ed, 
FiGuRE 3.—Apparatus for preparing known miztures. 


(b) PROCEDURE 


The procedure is as follows: Flask, connections, and cylinder 4 
are thoroughly dried and evacuated. The manometer is read, A and 
2 are closed, A is removed, weighed, replaced, and the connections 
evacuated. Cock 1 is then closed and 2 opened. C,H,O is allowed to 
purge through the atmospheric outlet of cock 3 and the safety seal, F, ' 
immersed in onginy and then slowly admitted to the flask maintaln- 


ing excess flow at Cock 3 is closed when the desired pressure s 
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attained, A is immersed in liquid air, its valve opened, and the 
(,H,0 condensed therein. The pressure is observed, the valve closed, 


liquid air placed around the tube, D, at the bottom of the flask, and 
the pressure again observed. C,H,O and CO, will condense to negli- 
sible pressures at —190° C; but very small traces of air in these gases 
wil prevent their complete transfer to A within reasonable time, by 
accumulating at the opening of the needle valve, and offering a plug 
of uncondensed gas through which the final small portions of C,H,O 
and CO, must diffuse slowly. Most of the gas condensed rapidly in A, 
but the final traces were transferred slowly. Rather than wait for the 
tardy diffusion to occur, a gravimetric correction was made for 0.2- to 
4mm residual pressure. Actually, this introduced the greatest 
uncertainty in the measurements, and would be avoided in more 
exacting work. Cock 2 and A are both closed, and A is removed 
and weighed. Repeated weighings were made to insure reproducible 
temperature and humidity conditions, and against leakage of valve. 

After weighing cylinder A it is replaced, the connections are evacu- 
ated, and the C,H,O is condensed in D. A is closed, and the C,H,0O is 
allowed to evaporate in C. The pressure and temperature are ob- 
served. CO, is now purged through 3 and the desired amount is 
then admitted into C. The pressure and temperature are observed. 
The two pressure observations from which the composition is cal- 
culated are thus made within a relatively short interval of time, and 
temperature conditions remain comparable. The C,H,O+CO, is 
then condensed into A at —190° C. The pressure is observed, A is 
closed, D is immersed in liquid air, and the pressure is again observed. 
Cock 2 is closed, and A is removed and weighed. 

Ais then replaced, the connections are evacuated, and the mixture 
iscondensed in D. Cock 2 is closed, the mixture is allowed to evapo- 
rate and is stirred by thermal convection, E replaces A, the connec- 
tions are evacuated, and a sample taken into Z from C for analysis. 


(c) RESULTS 


The results of the determinations made with this apparatus and in 
this manner are given in table 3. Five mixtures were investigated, 
ance some differences might be expected because of the different 
pressure-volume relationships existing. The mixtures were made up 
it different temperatures (26° to 30° C), and analyzed at different 
temperatures (25° to 31° C), so that any probable effect of tempera- 
ture 1s accounted for in the data. The range studied covers what may 
be expected in commercial use. 
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TABLE 3.—Comparison of the composition of five mixtures of C,H.O+ R 
determined by three methods: fF C1H0+CO, ag 


M=Amount of C2H«O, expressed as mole percent calculated from observed gravimetric data: 
100 No, 4,0 
Total N 
P=Amount of C2H,0O calculated from observed pressures: 
100 Po,H,0 
Total P 
V=Amount of C2H,O, expressed as percent by volume as determined from gasometric analysis: 
100 Voc,H,0 
Total V 
If the ideal gas laws were obeyed, these three would all be expressions of mole percent, and would 


agree within the limits of error of the observations, assuming the reaction involved in the data V were 
quantitative. t 





Percentage of C:H4,O Percentage of C:H,0 





Mixture number 
P M-P | M-V 





_—. 


—0. 04 | +0, 02 
—.07 | 

—.05 | 
—.04 | 
0 

















| 
| 
| 
Algebraic averages —0.04 | 





* Within the significance of the present work, these values may be taken directly as an expression of per- 
centage by weight, since the molecular weights of the two gases are very nearly the same. 


The data given show the differences existing between the compos- 
tion, expressed as percentages of ethylene oxide, as determined (.\/) 
gravimetrically according to the method just described, (P) by the 
measurement of pressure at fixed volume, and (V) by actual gaso- 
metric analysis, according to the method described in this paper. 

A glance at these data discloses the following facts: 

1. The amount of C,H,O, as determined by observed pressures, is 
—— about 0.04 percent higher than that measured gravimetn- 
cally. 

2. The amount of C,H,O as measured by the gasometric analysis, 
is in general about 0.05 percent lower than that measured gravi- 
metrically. 

3. The accuracy attained by the gasometric analytical method, in 
the range of compositions investigated, is commensurate with the 
reproducibility of the method. The average deviations observed over 
the range of composition investigated may be applied as corrections 
of the gasometric results (if desired) in computing mole percent or 
percent by weight. 


VI. NOTE ON SAMPLING 


The fumigant is usually obtained as a mixture of liquid C,H,O and 
CO, under pressure in a cylinder. Some cylinders are equipped with 
a connection leading from the valve to the bottom, so that the contents 
may be discharged from the liquid rather than the gas phase. Even 
so, one must not expect to open the cylinder valve slowly and draw off 
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, sample of gas whose composition approaches that of the liquid 
within the cylmder. The two gases have widely different boiling 
points, and will separate notably by distillation when vapor is drawn 
Jjowly from the cylinder; and vapor will be drawn from the cylinder 
t the seat of the valve is warm enough to vaporize CO,. This vapor 
ill push back liquid trying to escape, and so the C,H,O will be 
afectively separated from the mixture. . 

The tendency of liquified mixtures to separate under such conditions 
has long been known and has been studied experimentally. An inter- 
sting example of such behavior was observed during the gravimetric 
preparation of the mixtures here reported. After weighing one of 
these mixtures, the cylinder containing the liquid (at room tempera- 
ture) was connected to the evacuated flask, the valve was opened, and 
the mixture allowed to vaporize into the flask. The pressure within 
the flask built up immedzately to 10 to 20 mm in excess of the pressure 
of the CO, alone. Thereafter, the pressure increased very slowly to 
qual that originally observed for the mixture of C,H,O+CO,. 

The only recourse is an extremely rapid discharge of the cylinder 
yith valve opened wide) into a long pipe of sufficient diameter so that 
hack pressure or cooling does not cause condensation. Samples may 
be taken from this discharge tube. Better yet, if the cylinder is dis- 
charged completely into a large dry evacuated container used for 
fumigating, samples from this may be treated with more confidence. 

It is also necessary to take the sample in dry containers; and rubber 
tubing should not be used. 


WasHInaton, July 27, 1938. 
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ELASTIC PROPERTIES OF CAST IRON 
By Alexander I. Krynitsky and Charles M. Saeger, Jr. 


ABSTRACT 


An optical method for measuring the deflection of cast-iron transverse-test 
pars during loading and up to the breaking strength has been developed. ‘Trans- 
verse-strength properties were determined on test bars made from three types of 
cast iron heated to the maximum temperatures of 1,400°, 1,500°, 1,600°, and 
1700° C. Test bars were vertically cast bottom-poured in green-sand molds 
at 100°, 150°, 200°, and 250° C above the liquidus temperature. Total, plastic, 
and elastic deflection; modulus of rupture; modulus of relative elasticity; and 
total, plastic, and elastic resilience were determined and the microstructure of 
the test bars was examined. 
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I. INTRODUCTION 


The term “elastic properties” as used in this discussion connotes 
those complex properties that determine the behavior of cast iron 
under load. They may be determined from the results of transverse 
bending tests. 

All of the properties of cast iron are related to and, to a large extent, 
determined by the structure, which consists of a coherent metallic 
matrix throughout which particles of graphite are dispersed. The 
presence of these particles decreases the amount of metallic material 
wailable to resist stressing in a given section and causes nonuniform 
distribution of stresses, but the quantitative interpretation of these 
tects is as yet unsettled. In the consideration of the behavior of 
cast Iron in compression, the question whether the spaces occupied by 
the graphite particles should be considered as voids or as spaces filled 
with compressible, partly compressible, or incompressible material, 


191 





192 Journal of Research of the National Bureau of Standards ya.» 


has been discussed by several authors, for example, MacKenzie (1): 
Thum and Ude [2], Meyersberg [2a], Pearce [3], and Bolton [4]. Aside 
from the general weakening or depreciation of properties, because of 
the presence of the relatively soft graphite flakes, the nonhomogeneoy 
structure frequently leads to erratic results in determinations of the 
properties and renders their interpretation uncertain. Such factors 
as the composition, shape, size, and surface condition of the specimen 
and the details of the testing procedure are of greater importance jn 
determining the properties of cast iron than in similar determinations 
on other metals, steel, for example. ‘ 

The determination of the elastic properties of cast iron is further 
complicated by the fact that even moderate stresses produce plastic 
as well as elastic deformation in this material. Attention was called 
to this fact approximately 50 years ago, by de Segundo [5] and by 
Bach [6]. The observations of these early investigators have been 
confirmed by subsequent workers, and it had been definitely estab. 
lished that plastic deformation of cast iron begins under low loads, 
In the determination of elastic properties by means of transverse tests 
it is necessary, therefore, to define complete stress-deflection curyes 
by means of step-wise loading, and to determine for each load not only 
the total deflection under load but also the plastic deformation or 
permanent set that remains in the bar after the load has been removed, 
The difference between the total and the plastic deformation is the 
elastic deformation. 

Interpretation of the results of transverse tests and of the data 
computed from these results, for example, modulus of elasticity, 
resilience, et cetera, in terms of properties such as brittleness, tough- 
ness, ductility, and energy required to break a bar has been widely 
discussed but the correlations have not been entirely satisfactory 
partly because of the difficulties of conducting transverse tests, with 
resultant uncertainty regarding the accuracy and reliability of the 
data obtained. 

Preliminary results of this investigation, presented in a previous 
publication [7], showed that various properties, including the trans- 
verse strength and running quality during casting as well as the size 
and distribution of the graphite flakes in each of the irons used, were 
affected by the maximum temperature to which the liquid metal had 
been heated. Computations of the modulus of elasticity were made, 
based on the data obtained under loads of 1,200 pounds and also under 
the loads that caused rupture. These results clearly showed the need 
for a more detailed and systematic study of the stress-strain relations 
in cast iron. For the present investigation, it was decided to prepare 
test bars from each of several irons, representing different maximum 
heating and pouring temperatures, and to record, as accurately as 
possible, the complete stress-strain curves up to the breaking point of 
each bar, by means of observations of the successive deflections under 
increasing loads and the permanent set that accompanied the applica- 
tion of each load. 


II. METHOD OF MEASURING DEFLECTION 


_ One of the first requirements for the proposed investigation was al 
improved method of measuring the deflection of the bar under load 


1 Figures in brackets indicate the literature references at the end of this paper. 
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nd the permanent set in the bar after the load was removed. Direct 
spservation of the changing position in space of a reference mark on the 
bar is unsatisfactory because it involves a number of assumptions, for 
sxample, the supports and the testing machine as a whole must remain 
movable and there must be no penetration of the supports into the 
secimen. ‘To avoid these possibilities of error a new method of 
neasuring deflection was devised, as follows: 

In this new method of measuring deflection of bars under transverse 
iyading, a rubber band, B, is stretched tightly along the length of the 
gr, A, and is held in position by two clamps, C, as shown in figure 1. 
\fetal strips on the inner surface of the clamps keep the rubber band 
year to, but not in contact with the surface of the bar at all times. 
The spacing between the two clamps is the same as that between the 
sipports for the bar in the testing machine, 18 inches in these experi- 
neits. Reference mark D is attached directly to the bar, midway 
petween the clamps. ‘The bar is placed in the testing machine and the 
nicrometer telescope, /’, mounted at a distance of 20 inches from 
the bar, is used to measure the distance between the lower edge of 
therubber band and the top edge of the reference mark. When the 
har deflects under load, the reference mark moves away from the rubber 
hand which remains as a straight line connecting the central portions 
of the bar at the two supports. The distance between the rubber 
band and the reference mark therefore indicates the deflection of the 
jaded bar and a similar measurement after the load is removed 
indicates the permanent set. With this apparatus, the maximum 
ror in reading deflections, including errors in the micrometer tele- 
ope and personal errors of the observer, was found to be less than 
+(.002 inch. 

The principal advantage of :this simple method of measuring de- 
fections is that the readings are made entirely independently of the 
esting machine and supporting fixtures. A further advantage is that 
observations can be continued up to and including the breaking load 
without danger to sensitive optical apparatus or to delicate in- 
struments such as extensometers. 


Il. MATERIALS AND PREPARATION OF TEST BARS 


Three types of iron, of the compositions shown in table 1, were used 
in this investigation. Iron A is a stove-plate pig iron remelted with 
the addition of 10 percent of open-hearth ingot iron. Iron B approxi- 
nates the composition of a medium cylinder iron, although the carbon 
sslightly high and the phosphorus is somewhat low. Iron C repre- 
sents a soft iron such as is used for general castings. 


TaBLE 1.—Composition of stock pig irons 





Total Silicon Manga- 


carbon nese Phosphorus}; Sulfur 





Percent Percent Percent Percent Percent 

3. 44 1. 40 0.15 0. 46 0. 020 
3.79 1. 40 - 63 . 181 . 062 
3. 44 2. 43 ote . 395 . 050 
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Cylindrical bars for the transverse tests were cast in groups of fo, 
as shown in figure 2. The pattern was molded in a three-part eylindn. 
cal flask with the full length of the test bars contained in the check. 
Each flask contained two groups, that is, eight bars. The patter, 
were molded in a mixture of eight parts of molding sand and one part of 
sea coal, tempered to approximately 7 percent of moisture. The mold 
cavities were faced with nongraphitic carbonaceous material: ¢/ 
commercial origin. The castings were bottom poured with the bar 
occupying vertical positions in the mold. By this procedure, as do. 
scribed in a previous publication [8], test bars can be produced whic) 
are satisfactorily uniform in diameter throughout their length and free 
from burnt-on sand or other surface defects. The cylindrical sectio, 
of each bar was 21 inches in length and 1.2 inches in diameter. : 

Charges of 230 pounds of iron, in a commercial magnesia crucible 
were melted in a high-frequency induction furnace of the tilting type. 
Each melt was heated to a predetermined maximum temperature o{ 
1,400°, 1,500°, 1,600°, or 1,700° C and was maintained at this tem. 
perature for approximately 1 minute. The metal was then allowed 
to cool until the temperature was 250° C above the liquidus tempers. 
ture ® of the iron used in that particular melt. A set of four bars was 
then cast as a unit by pouring the metal directly from the tilted fur. 
nace into the mold. The metal remaining in the crucible was allowed 
to cool to the next pouring temperature, 200° C above the liquidus, 
and another set of four bars was cast. By continuing this procedure, 
a set of four bars was cast at 150° C above the liquidus, and a final 
set at 100° C above the liquidus. Sixteen test bars, representing four 
pouring temperatures, could be thus obtained from a 230-pound melt 
that had been heated to a definite maximum temperature. The bars 
were allowed to remain in the molds for at least 18 hours after casting. 

Temperature measurements of the molten iron up to 1,600° C were 
made with a platinum—platinum 10-percent rhodium thermocouple 
that was protected by a closed-end, glazed porcelain tube inserted in 
a closed-end graphite tube. The portion of the graphite tube that 
came in contact with the molten iron was protected by a coating of 
aluminum oxide covered with a mixture of 95 percent of zirconium 
silicate and 5 percent of bentonite. 

Temperatures above 1,600° C were measured with an optical pyrom- 
eter which had been previously calibrated with the thermocouple 
in the temperature range of 1,400 to 1,600° C. 


IV. TESTING PROCEDURE 


A universal testing machine of the hydraulic type of 50,000 pounds 
capacity was used for the transverse tests. The rate of loading of 
all bars corresponded to 0.12 inch per minute travel of the free cross 
head of the testing machine. 

The test bar, equipped with the rubber band and reference mark, 
as shown in figure 1, was mounted in the testing machine on two sup- 
ports of the roller type, at a span of 18 inches. The bar was adjusted 
until the rubber band lay in the horizontal plane of the neutral axs 
of the bar. A load of 50 pounds was sufficient to seat the bar firmly 
but was not enough to produce noticeable plastic deflection, was p- 


: Speeninste analysis: 4 percent of volatile; 74 pone of fixed carbon; 22 percent of ash. 


e ao temperature of each iron was estimated from its composition, according to the procedure 


described in a previous publication [7]. 
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Apparatus for measuring the de fle ction of test bars under transverse 


loading. 


* bar: B. rubber band; C, clamps; D. reference mark; and FE, micrometer telescope 
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Figure 2.— Transverse lest bars, as cast. 
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P n: B. down gate; C, semicircular feeder for bottom pouring; and J, verti 
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and the distance between the edge of the rubber band and the 
top of the reference mark was determined as a zero reading. The bar 
was then loaded in a series of steps (100 or 200 pounds each) until the 
preaking load was reached.and the deflection of the bar under each 
of the loads was recorded. After the load reached 800 or 1,000 pounds, 
each application of a load was followed by unloading to the original 
50-pound load, and the permanent set in the unloaded bar was deter- 
mined. Attempts to record the small permanent deformations that 
resulted from loads of less than 800 or 1,000 pounds were seldom made. 
Unloading the bars did not affect the deflection under subsequent loads; 
the breaking loads and the load-deflection curves of duplicate bars 
were practically identical when one of the bars was loaded in the 
manner just described and the other was progressively loaded without 
unloading between successive loads. 


V. PHYSICAL PROPERTIES OF TEST BARS 


plied 


The results of the transverse tests yielded information on the fol- 
lowing properties of the cast iron bars: Modulus of rupture; plastic, 
dastic, and total deflection; relative modulus of elasticity; and 
resilience. 

1. MODULUS OF RUPTURE 


The modulus of rupture of each bar was calculated, according to 
the formula 


2.546LS 
M= ~~ P”? 





+ 


MAXIMUM HEATING TEMPERATURE 














350 1400 1450 15001350. 1400 1450 1500 1300 1350 1400 1450 
POURING TEMPERATURE °C 


Figure 3.—Effect of maximum heating temperature and pouring temperature on the 
modulus of rupture of three irons. 


where M is the modulus for round bars; L, the span of the bar between 
supports; S, the breaking load; and D, the diameter of the bar. 

The modulus of rupture for each of the three irons and the effect on 
the modulus of varying the maximum heating and pouring tempera- 
tures are shown in figure 3. Each of the plotted points is an average 
of values ‘derived from four companion bars. 
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The modulus of rupture of iron A was found to be considerably 
higher than that of either of the others; and that of iron B to be the 
lowest of the three. For all of the irons the modulus increased with 
an increase in the maximum heating temperature and tended to 
decrease with an increase in pouring temperature, the latter effec 
being more pronounced for iron A than for B and C. 


2. PLASTIC, ELASTIC, AND TOTAL DEFLECTION 


‘ 


The total and plastic deflections or ‘set’ of each bar under jp. 
creasing loads, were determined and the elastic deflection was taken 





T 


iRON A 
WEATED TO 
POUREO aT 





iRON B 
HEATED TO 1700 °C 
POURED aT 1460 °C 


1RON C 
MEATEOD TO 1700 
POURED aT 1350 °C 


ee ee | ee | 


10 5 20 








DEFLECTION INCH 


Figure 4.—Typical plastic, elastic, and total deflection curves of irons A, B, and ( 


as the difference between the total and the plastic deflection for each 
load. ; ™ 
The relation between load and total, plastic, and elastic deflections 
of representative bars of irons A, B, and C is shown in figure 4. Tron 
A has the highest breaking strength of the three irons and has about 
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the same deflection at the breaking load as iron C. Iron B has about 
the same breaking load as iron C but considerably more deflection 
at rupeure. 

In all cases the total deflection curve is a continuous curve inclining 
progressively towards the deflection axis. If there is any linear 
portion in the total deflection curve it is confined to the region of 
very light loads. The curves for plastic deformation for loads of 
800 pounds or more, likewise, are continuous curves without straight 
portions. These observ ations are in accord with those of previous 
precenee but the data for elastic deflection are not. Meyersberg 

jconcluded that the upper portions of the elastic elongation curves 
were reslightly bent, whereas Pearce [3] found that the elastic deflection 
curve usually was a continuous straight line for all loads. The results 
this study, illustrated in figure 4, show that the elastic deflection 
wrve is linear only in its lower portion. For iron A, the linear 
portion extends to some load between 1,400 and 1,800 pounds for 
agch of the bars. For irons B and C it extends to the range of loads 
- 1,200 to 1,400 pounds. Above the linear portion, each curve 
nes toward the deflection axis. 

“The departure from linearity in these curves can be demonstrated 
nathematically by the method suggested by Tuckerman in the dis- 
cussion of a paper by Templin [9]. Theoretical values for the elastic 
deflection for each load can be calculated from the equation for a 
straight line and from the modulus of elasticity of the bar at a load 
within the initial, straight portion of the curve. As long as the curve 
conforms to the equation for a straight line the observed and caleu- 
lated values will be in agreement, within the limits of experimental 
eror. Any systematic or appreciable divergence of the observed and 
calculated values indicates a deviation of the curve from the equation 
fora straight line. 


TABLE 2.—Observed and calculated elastic deflection for iron B 


pari f observed values for elastic deflection with values calculated from the relative modulus of 
ty, 12,445,000 Ib/in.’, of this bar under a load of 1,000 pounds. 


Elastic deflection 


Load 7 ee eae Sale ais ee 
Observed | Calculated | Difference 


Pounds Inches | Inches Inches 
700 0. 064 0. 066 —0. 002 
800 .075 .075 . 000 
900 . 083 . 085 . 002 

, 000 . 094 . 094 
, 100 


, 200 
, 300 
, 400 
| , 500 
600 | 
| 
| 1, 700 
i; 1,800) | 
| 
| 
= 





1, 900 | 





These computations were made for a number of the curves. The 
data in table 2 show the results obtained from a bar of iron B heated 
‘0 & maximum temperature of 1,700° C and poured at 1,460° C. 

ads were applied to this bar in increments of 100 pounds, to define 
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the curves more closely. It is evident that the observed and calq) 
lated values for elastic deflection are in agreement, within the limits 
of experimental error, +0.002 percent, for loads of 1,300 pounds a 
less, but not for loads in excess of this figure. The elastic deflection 
curve of this bar, therefore, is a straight line only for loads of | 399 
pounds or less. bia 


3. RELATIVE MODULUS OF ELASTICITY 


In this investigation the term “relative modulus of elasticity” j; 
used in place of “modulus of elasticity,” because it is well known that 
cast iron does not follow Hooke’s law exactly and the term “modulus of 
elasticity,’ as applied to cast iron, generally means the relative 
stiffness of the iron under the particular conditions of loading. 

The relative modulus of elasticity, as determined in the transyerse 
test for a single centrally applied load, is determined by the formulg 


ee 
E=381D" 
where 
E=relative modulus of elasticity in pounds per square inch, 
P=load in pounds, 
J=span in inches, 
IJ=moment of inertia, and 
D=deflection (at load P) in inches. 


Values for D were taken from the elastic-deflection curves. In view 
of the fact that the elastic-deflection curve was found to be linear only 
in the lower portion, two values of the relative modulus of elasticity 
were computed for each bar; E, for the straight portion of the elastic. 
deflection curve, and FE, for the second portion of the curve, at the 
point of rupture. 

The data in table 3 show that the values for H, are consistently 
higher than for Z,. For irons A and B the difference between F, and 
E, is generally greater than 1,000,000 lb./in.?; for iron C the difference 
is usually less than this figure. 

Values for the relative modulus of elasticity, Z,, have been plotted 
in figure 5 to show the variations with different maximum and pouring 
temperatures. In general, values of F, for iron A are higher than those 
for either of the other irons and, under some conditions of heating and 
pouring, the values of , for iron B are higher than those for iron C. 
A tendency for the values of /, to increase with an increase in maxt- 
mum heating temperature is more evident in irons A and B& than m 
iron C. The data for each iron show a slight and somewhat erratic 
tendency for the values of EZ, to decrease with increasing temperature 
of pouring. 
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TaBLE 3.—Relative modulus of elasticity of cast iron 
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‘These data are average values of 2 bars. 


4. 


All other figures are averages of 4 bars. 


RESILIENCE 


The resilience of cast iron, the work absorbed by the specimen 
before rupture, is obtained by measuring the area below the deflec- 
lion curve and between the origin and an ordinate to the deflection 


curve at the point of rupture. 


If the resilience of a specimen could 


¢ simply and adequately expressed there would be no need for 
discussion of “stiffness”, “plasticity”, ‘“‘toughness’’, “brittleness’’, or 
‘miler indefinite but important properties of cast iron. These 
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values measure the total work but do not define the way in which js 
has been absorbed. It is necessary, therefore, to consider not only 
the area below the curve but also its shape, to determine whether ; 
particular resilience value is the result of a high, low, or mediun 
breaking strength coupled with high, low, or medium deflectioy 
And, finally, a complete definition of resilience should differentis, 
between the work absorbed in elastic and in plastic deformation, _ 

Measurements of the total resilience are best determined with , 
planimeter; determinations of the ‘‘triangular’’ resilience, based oy 
the assumption that the total deflection curve is a straight lino 
between the origin and the point of rupture, obviously only approxi- 
mate the true resilience and the degree of approximation decregsps 
as the curvature of the deflection curve increases. 

Values of total resilience of the three irons, measured planimetrically 
are shown in figure 6. The resilience values of irons B and ( gy» 
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Figure 6.—Effect of maximum heating and pouring temperatures on the total resilience 
of cast iron. 


similar. In neither case is there a definite relation to either the maxi 
mum heating temperature or the pouring temperature, except tha 
very low resilience values were obtained from iron C that had bee 
heated only to 1,400° C. The curves for iron A are quite differen 
from those for irons B and C. For iron A there is a pronounce 
decrease in resilience with increasing pouring temperature, at least 
up to 1,450° C. An increase in men, te with increasing maximum 
heating temperature is evident in the bars that were poured at 1,300 
C, but this effect is less evident at higher pouring temperatures. Iron 
A under some conditions of melting and pouring showed the highest 
resilience values of any of the three irons, but under other conditio! 
of melting and pouring its resilience values were lower than any 0 
those of irons B and C. These resilience values illustrate the marked 
effect on the characteristics of cast iron of an addition of 10 percent 
of ingot iron, as in iron A. i. 

A number of attempts to establish more satisfactory criteria fo 
evaluating the properties of cast iron have been made in recent yeals 
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MacKenzie [10] found that there was a definite relation between the 
ratio of the true to the triangular resilience and the ratio of the secant 
modulus of elasticity at half load to the ultimate modulus of elasticity. 
These ratios were computed from the data in this mvestigation, by 
using values of H, for the secant modulus at half load and values of 
F) for the ultimate modulus, but no definite relation between the 
widely scattered plotted points could be established. 

Attempts to correlate resilience values with the carbon and silicon 
contents also have been made. For example, MacKenzie [1] and 
Bolton [4] compared resilience values with those representing the 
carbon content, plus one-fourth of the silicon content. Values of 
('+Si/4, for irons A, B, and C were computed and were plotted against 
total resilience, as shown in figure 7. The plotted points for each 
iron are widely scat- 
tered, the spread of 
the resilience values IRON A 
being greater than ‘RON B 
that of the C-+Si/4 IRON C 
values. However, 
there is some indi- 
cation of an increase 
in total resilience 
with increasing val- 
ues of C+Si/4, as is 
indicated in figure 7 
by the two dotted 
lines that include 
the majority of the 
points for all three 
irons. A similar ene! sa 
relation was found 300 340 380 
for plastic resilience TOTAL RESILIENCE INCH POUNDS 
values, but the elas- 
tic resilience tended 
to decrease with 
increasing values of C+Si/4. On the whole, it must be concluded 
that the relation between resilience and carbon-plus-silicon contents 
is not definite for these three irons. 

Tucker [12] concluded that neither the breaking load, deflection, nor 
resilience values gave reliable indications of the resistance to thermal 
shock, but that the plastic resilience, expressed in percentage of the 
total resilience, furnished a measure of the relative toughness of the 
irons that he investigated. Such values of plastic resilience of the 
bars in the present investigation were obtained by subtracting from 
the total resilience the elastic resilience computed on the assumption 
that the elastic-deflection curve was a straight line between the 
origin and the breaking point. These values, expressed as percentage 
of the total resilience, are given in table 4. The data show that 
the plastic resilience values for each iron are slightly lowered when the 
maximum heating temperature is 1,600° or 1,700° C; the effect of 
varying the pouring temperature, after heating to a particular maxi- 
mum temperature, 1s not definite. 

According to the data in table 4, the ratio of plastic to total resili- 
ence (in percent) of iron B is definitely higher than that of either of 
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Ficure 7.—Relation of total resilience to carbon and silicon 
contents. 
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the other irons, with the values for iron A slightly lower than for iroy 
C. According to Tucker’s data this indicates that the medium 
cylinder iron B is definitely more resistant to thermal shock thay 
either of the other two, and that the soft iron C is somewhat better 
than the high-strength iron A in this respect. It is interesting to 
note the grouping of the three irons classified in this way and accord- 
ing to the total resilience. On the basis of a comparison of the per. 
centage plastic resilience, iron B is superior with irons A and C paired 
whereas, if the comparison is based on the total resilience (fig. 6), irons 
B and C are quite similar and A is different. These results indicate 
the complexity of resilience determinations and the difficulty of their 
interpretation. 


TABLE 4.—Relation of plastic to total resilience of three irons 
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® Average value of 2 bars. Other figures are average values from 4 bars. 


VI. EFFECT OF MAXIMUM HEATING TEMPERATURE ON 
COMPOSITION AND STRUCTURE 


After completing the transverse tests, the bars were sampled for 
chemical analysis and for microscopic examination. The analytical 
data are summarized in table 5, which gives average values, for each 
iron, for aJl the metal that was heated to a particular maximum 
temperature, irrespective of the temperatures at which the bars were 
poured. Although data on this point are not included in table 5, 
variations in pouring temperature did not appreciably affect the 
content of any of the elements determined. ™ 

The data show that there was no appreciable variation in the silicon, 
manganese, phosphorus, and sulfur contents with variations In 
maximum heating temperature. There was, however, a definite 
tendency for the total carbon to decrease with increasing temperature 
of the molten metal. The metal was melted in a magnesia crucible 
with the surface of the melt exposed to the air. Under these condi- 
tions some loss of carbon by oxidation at the higher temperatures 
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would be expected. The total carbon contents of all the bars are 
close to or slightly less than the eutectic values computed from the 
mounts of silicon and phosphorus present. The combined carbon 
contents of bars of iron A are close to the eutectoid composition for 
‘he low-strength bars and exceed the eutectoid composition for the 
hich-strength bars. In all bars of irons B and C the combined carbon 
is less than the eutectoid amount. 


fapLe 5.— Effect of maximum heating temperature on the composition of cast-iron 
bars 
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Corrected to bar diameter of 1.2 inches. 


In general, the strength of cast iron increases with decreasing 
amounts of graphitic (or total) carbon. Data for the breaking loads, 
in table 5, show that the bars of iron A had higher strengths and lower 
graphitic carbon contents than the corresponding bars of either B or C. 
Furthermore, in each of the three irons, an increase in maximum heating 
temperature decreased the amount of graphitic carbon and increased 
the strength. Evidently the strength and associated properties of 
these irons are affected by the amount of graphitic or total carbon 
but examination of the microstructure shows that there are variations 
in the size of the graphite flakes and in their distribution, as well as in 
the amount of graphitic carbon. 

Cross sections of the test bars, from locations close to the fractures, 
were prepared for microscopic examination, according to the pro- 
cedure described by Ellinger and Acken [13]. The specimens were 
ground on a surface grinder, lapped on a lead-tin plate charged with 
ine emery, further lapped on a second plate charged with finer emery, 
and, finally polished with an aqueous suspension of rouge, in an auto- 
matic polishing machine. 

_ The micrographs in figure 8 show that the high-strength bars of each 
iron contained smaller and more uniformly distributed particles of 
graphite than did the low-strength bars of the same iron. The dif- 
erence in particle size is more readily apparent in the micrographs of 
representative areas under higher magnification, in figure 9, and the 
structure of the nongraphitic matrix of each bar, revealed by etching 
the polished specimen, is shown in figure 10. The latter micro- 
graphs show that the matrix of the low-strength bars of iron A (heated 
to 1,400° C) is coarsely pearlitic, whereas in the high-strength bars of 
the same iron (heated to 1,700° C) large amounts of massive cementite 
ure present, which, by stiffening and hardening the matrix, undoubt- 
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edly contributed to the high strength of these bars. Cementite jg 
found only in iron that contains hypereutectoid amounts of combine 
carbon; pearlitic structures are associated with eutectoid and hypo. 
eutectoid compositions. For iron B, the matrix of both high- and 
low-strength bars was found to be pearlitic, with small amounts oj 
ferrite present and the structure of both high- and low-strength by: 
of iron C was found to consist of fine pearlite with areas of ferrite gy, 
of steadite. 

The relations between the strengths, compositions, and structures 
of the three irons may be summarized as follows: 

1. The increase in strength of each of the irons after it has bee 
heated to a high temperature is associated with a decrease in the 
amount of graphitic carbon and in the size of the graphite particles 
and with an increase in the uniformity of their distribution. 

2. The high strength of iron A, as compared to irons B and (. js 
associated with lower graphitic carbon and higher combined carbo) 
contents. 


VII. UNUSUAL STRUCTURAL FEATURES IN CAST IRON 


An unusual structural condition was observed in many of the 
graphite flakes of iron C, in which an appreciable portion of the 
graphite appeared to be subdivided into roughly hexagonal grains, 
This condition is readily apparent in the micrograph of low-strength 
iron C, in figure 9, and is also evident, although less pronounced, in the 
micrograph of a high-strength iron C. Nipper [14] observed similar 
structures in samples of graphite laminas, but it is not entirely clear 
from his paper whether these laminas are natural graphite or have 
been isolated from cast iron. The present example of hexagonal 
graining in graphite refers definitely to the structure of graphite in 
cast iron. 

Several of the specimens of iron C also showed a tendency towards 
a dendritic structure in the arrangement of graphite particles. This 
tendency is evident in portions of the micrograph of the high-strength 
bar of iron C, in figure 8, but in no case was the dendritic structure 
developed to such an extent that the properties of the bar were 
definitely affected. 

Several of the micrographs in figure 9 suggest the existence of « 
laminated structure in some of the graphite particles. [Examination 
under higher magnification confirmed this suggestion by showing that 
laminated areas existed in some of the large graphite particles and 
that the structure of some of the smaller particles was almost entirely 
laminated, as is shown by the micrographs in figure 11. This condi- 
tion was most pronounced in specimens from high-strength bars 0! 
iron A, in which the combined carbon is hypereutectoid, but was also 
observed in hypoeutectoid, low-strength bars of iron A and occasion- 
ally in both high- and low-strength bars of iron B. These laminated 
structures were observed in specimens from the unstrained ends of the 
bars as well as in specimens from locations adjacent to the fractures 
and, therefore, were not caused by stress deformation of graphite 
particles during transverse loading. Furthermore, the laminated ap- 
pearance persisted through repeated repolishings and hence was not 
an accidental or surface phenomenon. 
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Discussion of the mechanism of formation of graphitic particles 
involves the controversial questions of the double iron-carbon dia- 
‘ram, that is, the question whether graphite in cast iron is formed at 
autectic or eutectoid temperatures, or both. : Arguments for and 
against the double diagram were reviewed by Epstein [15], who con- 
cluded that the existence of the iron-graphite eutectic had been con- 
irmed by observations of eutectic-like structures and that the 
existence of the iron-graphite eutectoid had been indicated by the 
results of thermal analyses, although no one had been able to illus- 
rate satisfactorily the appearance of the iron-graphite eutectoid. In 
a recent paper Wells [16] concluded that the iron-graphite eutectoid 
gecurs at 738° C and 0.69 percent of carbon, and that the graphitiza- 
tion observed in high-purity alloys of iron and carbon resulted directly 
fom decomposition of austenite as well as from decomposition of 
cementite. Boyles’ recent paper [17] indicates that the large flakes 
of graphite in cast iron were formed during the solidification at 
eutectic temperatures and that subsequent decomposition of austenite 
or carbide at subeutectic temperatures produced finely divided 
mauphite which had little effect on the flake structure. 
' The authors’ conclusion that the laminated graphite structure, 
illustrated in figure 11, was formed by eutectoid decomposition of 
qustenite or carbide is in agreement with the conclusions of these 
other investigators. Laminated structures have been observed in 
sme examples of eutectic alloys, but the best-known example of 
laminated structure is pearlite, which consists of alternate lamina- 
tions of ferrite and cementite in the proportions of the eutectoid 
composition. ‘The appearance and structure of the laminated graph- 
ite are similar to those of pearlite, and it is reasonable to assume 
that the laminated graphite also is a result of a eutectoid decomposi- 
tion, one in which austenite (or carbide) decomposes to form laminated 
graphite and ferrite. The micrographs in figure 11 therefore are 
ifered as photographic evidence of the existence of the iron-graphite 
eutectoid. 

Etching of specimens that contained laminated graphite revealed 
luther unusual structures. Figure 12 shows the appearance of a 
luke of graphite and its surrounding areas in one of the high-strength 
pars of iron A, etched with a 1 percent solution of nitric acid in ethy] 
alcohol. ‘The graphite flake is surrounded by pearlitic areas the struc- 
ture of which extends into, and at one point entirely across, the graph- 
ite. The appearance suggests that the surface of the pearlitic area has 
been flowed over the graphite in the polishing operation, but there are 
objections to this suggested explanation. The flow of metal surfaces 
durng polishing has been established by Beilby and subsequent 
workers, but it 1s generally agreed that such a flowed film is ‘“amor- 
phous,” that is, structureless, and hence different in properties from 
those of the parent metal. It is difficult to believe that such a layer 
could preserve the patterns that are shown in figure 12. Moreover, a 
reasonable explanation of this structure can be developed if it is 
assumed that both pearlite and the laminated graphite are produced 
by the eutectoid decomposition of austenite or carbide. Apparently, 
in the area shown in figure 12, the same eutectoid decomposition that 
produced pearlite in the surrounding areas formed the laminated 
structure in the graphitic area. The fact that some of the lamina- 
lions are common to both the pearlite and the laminated graphite, 
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indicates that the two structures were formed under the same condi 
tions. These observations further indicate that the laminated ora " 
ite structure shown in figure 11 represents the iron-graphite eutectoid 
mixture. Furthermore, the suggested mechanism of simultaneors 
formation of pearlite and laminated graphite in adjacent areas 9. 
counts for the observation of the existence of the laminated structyy 
in numerous small areas along the boundaries of large grains of graph- 
ite, as in figure 9. The structure of the pearlite is not revealed in 
these polished but unetched specimens. 

Although Bolton [18], in a discussion of a somewhat similar strye. 
ture, suggested that the black area might not be graphite, the authors 
believe that the black area in figure 12 is graphite, containing definit, 
indications of a laminated structure closely related to the structyr 
of pearlite. The present data do not indicate whether the laminate 
graphite is formed by eutectoid decomposition of austenite or of 
carbide. However, Wells [16] has shown that either of these reactions 
can produce graphite in high-purity iron-carbon alloys. 


VIII. SUMMARY 


1. A new method of measuring the deflection of test bars under 
transverse loading has been developed. Measurements of deflectioy 
under interrupted loading can be continued up to the breaking point 
of the bar and the accuracy of these measurements is entirely inde- 
— of the rigidity of the testing machine and the mountings o/ 
the bar. 

2. The properties of three types of cast iron were determined by 
measurements under transverse loading. Iron A was a stove-plate 
pig iron to which 10 percent of basic open-hearth ingot iron had been 
added, iron B represented the medium cylinder type of cast iron, and 
iron C was a soft iron of the type used for general castings. The 
average breaking load of iron A was appreciably greater than the 
breaking load of either irons B or C, and consequently, the elasti 
properties of A were generally superior to those of B or C. 

3. Variations in the maximum temperature to which the liquid 
metal was heated affected the composition, structure, and properties 
of the test bars. An increase in the maximum heating temperature 
decreased the amounts of total and graphitic carbon, decreased the 
size and increased the uniformity of distribution of the graphite par- 
ticles, and in general increased the transverse strength and elastic prop- 
erties. The properties of iron A were more susceptible than those of 
B or C to variations in the maximum heating temperature. 

4. Variations in the temperatures at which the test bars were 
poured had less effect than variations in maximum heating temper- 
ture on the composition, properties, and structure. However, there 
was a tendency for some of the properties to decrease with increasing 
pouring temperature, particularly for iron A. 

5. The results of the transverse tests indicate that the initial por- 
tion of the elastic-deflection curve is a straight line, but at approx 
mately half of the breaking load the curve inclines towards the deflee- 
tion axis. For all bars of iron A the points of inflection occurred at 
loads between 1,400 and 1,800 pounds, for irons B and C the limtt- 
ing value lay between 1,200 and 1,400 pounds. Because of the curva- 
ture of the elastic-deflection curves, a value for the relative modulus 
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Saeger 


of elasticity, computed from the data for any point on the initial linear 
ortion of the curve, is approximately 1,000,000 pounds greater than 
S smilar value computed from the load and deflection at the point of 


, e. 

PS tee unusual structural conditions were observed, indications of 
the existence of hexagonal-like grains within some of the graphite 
particles of iron C, and evidence of the existence of a laminated 
sructure, resembling that of pearlite, in many of the graphite parti- 
es of irons A and B, particularly the former. The laminated graph- 
ite structure is considered to be that of the iron-graphite eutectoid 
composition. 
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RELATION OF CAMERA ERROR TO PHOTOGRAMMETRIC 
MAPPING 


By Irvine C. Gardner 


ABSTRACT 


A portion of the errors in &@ map constructed from airplane photographs arises 
from differences between the actual performance of the airplane camera and the 
ystulated performance as defined by the calibration constants. The different 
sources of errors of this nature are enumerated and equations are derived showing 
the relation between errors in the calibration constants of the camera and the 
resulting errors in the map. 

Equations 3, 4, 10, 11, 18, and 19 give resulting errors in the map arising from 
incorrect V alue of the calibrated foeal length, distortion, and incorrect location of 
the principal point. Other sources of error in the camera, such as failure of the 

film to lie in a plane, the character of performance of the shutter, and imperfect 
filters are considered and their effects upon the image evaluated. 
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I. INTRODUCTION 
1. GENERAL CHARACTERISTICS OF THE AIRPLANE CAMERA 


The airplane mapping camera is a precision measuring instrument. 
Although fundamentally an optical instrument, it is dependent, for 
proper functioning, not only upon the proper design and installation 
of the optical parts but also upon the proper performance of a large 
number of elaborate and complicated mechanical auxiliary parts, 
among which are the mount which prevents the vibration of the air- 
plane from affecting the definition of the photograph, the film mech- 
anism by which the film is advanced and constrained to lie in the focal 
plane, the photographic shutter by which the exposure is timed, and 
the control devices by which the desired overlap is secured between 
the different negatives of a series. In this discussion consideration 
will be restricted to the fundamental requirements that are necessary 
in order that the resulting negatives may have the metrical character- 
istics suitable for the production of a map of the desired accuracy. 
These requirements relate to the precision with which the equivalent 
focal length must be determined, the distortion introduced by the 
objective, the location of the principal point of the focal plane, the 
planeness and shrinkage of the film, the type of shutter, and choice 
of focal length. After the tolerances connected with these character- 
istics have been derived and a suitable photographic objective has 
been chosen, it then becomes the problem of the camera designer to 
design the mechanical parts of the camera so that the photographic 
operations may be satisfactorily performed in a manner consistent 
with the maintenance of these tolerances. The requirements for the 
mechanical parts, except as they affect optical performance, will not 
be considered in this discussion. 

When use is made of an airplane photograph for the construction of 
& map, it is commonly assumed that the photograph is a true perspec- 
tive representation. Usually this is not rigorously true as distortion 
and other defects prevent the realization of a perfect perspective. 
Points A, B, O, C, D, and H in figure 1 are object points that are imaged 
on the photograph at A’, B’, 0’, C’, D’, and H’, respectively. The 
point N’ that brings the points of the photograph as nearly as possible 
into the correct perspective relationship should be chosen as the per- 
spective center. The point 0’, at the foot of the perpendicular from 
N’ to the plane of the photograph, is the principal point! and the 
length N’O’ is the principal distance. These terms are standard 
ones long used in the theory of perspective and are applicable to any 
representation in perspective, whether it is a photograph, painting, oF 


1 This principal point of the photograph or of the focal plane is not to be confused with the principal points 


of the lens, a pair of points coincident with the nodal points in lenses of the types employed for photogram- 


metric purposes. 
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drawing. It will be noted that the definitions are in terms of the 
photograph and bear no reference to the camera. 
The point N is the corresponding perspective center for the terrain 
‘hat is photographed. If B is the angle at N, subtended by O and 
any point A, and if y’ is the distance from 0’ to A’, then since, ideally, 
j=’, it is evident that 


‘di 


= 
tan B= 2 ’ 


(1 

















B oO Cc 
Fiaure 1.—The metrical relations of the camera. 


lhe symbols y’ and p designate the indicated lengths when D’, C’, O’, B’, A’, and IJ’ are assumed to be 
points on a print or diapositive. If they are assumed to be images of points, as formed in the camera (in 
which case they are unaffected by shrinkage of film or subsequent enlargement), the symbols y and f apply. 


where p is the principal distance. This relationship enables the 
photograph to be utilized for the measurement of angles. In practice 
these angles may be determined either by computation, after a direct 
measurement of y’, or, if a plotting machine is employed, the angle 
8is determined and utilized to transfer the desired result to the map 
ina semiautomatic manner without an intermediate reading or record- 
lug of the angle. 
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2. DIFFERENCES BETWEEN IDEAL IMAGERY AND ACTUAL 
CAMERA PERFORMANCE 


Returning to figure 1, let it now be assumed that points A’, 8’, 9 
CO’, D’, and H’ are no longer the points recorded on a photograph, es 
was the case in the discussion of the preceding paragraph, but that 
they are the image points as actually produced in the focal plane by 
the photographic objective. The distinction lies in the eliminatioy 
of possible changes in the spacing of the points as a result of film and 
paper shrinkage, subsequent enlargement from the negative, etc, |) 
order to make this distinction quite clear without making two separate 
drawings, alternate symbols, y’ and y, and p and f, are shown in figure |, 
When D’/7/’ is assumed to represent the trace of a print, made eithe; 
by contact or enlargement, y’ represents a length measured on the 
print and p is the principal distance. When points D’, C’, 0’, BY 
A’, and H’ are assumed to be points of the aerial image as formed by 
the camera lens, the symbols y and f refer to the corresponding di- 
mensions. The point N’ is selected as before to correspond to the 
most nearly correct perspective. The point N remains the corre- 
sponding perspective center for the terrain. We now have the equation 


tan B “ . (2) 


which resembles eg 1, the difference being that f has been written in- 
stead of p and that y’ in eq 1 is measured on the photograph and the 
y of eq 2 is measured on the image in the camera. In other words, 
eq 1 involves the properties of a finished negative, print, or enlarge. 
ment, whereas eq 2 is restricted to the properties of the camera. 

If it be assumed that the camera lens (not shown in the diagram) 
producing the imagery of figure 1 is free from aberrations and the in- 
evitable minor errors of workmanship, and further, if the idealized 
rays of paraxial imagery are assumed, then N and N’ are identical 
with the first and second nodal points o1 the camera objective and the 
quantity f of eq 2 is the equivalent focal length.? 

Inevitable departures from the idealized conditions of paraxial 
imagery, however, cause differences, too large to be neglected, be- 
tween the magnitudes that are actually used in the construction of 
maps and the corresponding quantities for the idealized case. 

Ordinarily the surfaces of a lens are assumed to be surfaces of 
revolution having their axes of symmetry in coincidence. This com- 
mon axis is referred to as the axis of the lens. When a lens is carefully 
examined on the test bench, it is often found, even though the lens 1s 
of excellent quality and repute, that the axes of the different surfaces 
fail to coincide exactly. In such a case it is difficult to define the axis 
of the lens or the positions of the nodal points with precision. Ifa 
line approximately coincident with the axes of the several surfaces 1s 
tentatively accepted as the axis of the lens, it may be found that a ray 
in the object space traveling along this line is deviated several min- 

2 This is based on the assumption that a, the altitude of the camera, may be considered infinite. When 
is 15,000 feet and the focal length of the camera objective is 6 inches (150 mm), the displacement of the image 
from the rear focal plane, because of the finite, though large distance to the object, is approximately two tel- 


thousandths of an inch, an amount negligible in comparison with the other departures from the ideal con- 
dition that are to be considered. 
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utes by transmission through the lens. Consequently, N and N’ of 
ggure 1 have been referred to, not as nodal points of the lens but as 
centers of perspective in the object and image space. The principal 
point of the focal plane is not necessarily the point where the axis of 
the lens intersects the image plane but is at the foot of the perpen- 
dicular dropped from the center of perspective in the image space to 
the image plane. For the pair of conjugate points, N and N’, angles 
3 and 6’ are equal when distortion is absent. These points may, 
therefore, be termed the points of unit angular magnification, al- 
though they satisfy this condition rigorously only in the ideal case. 
4s a consequence of this unit angular magnification, a lens may be 
tipped about N’ through an appreciable angle without seriously 
affecting the metrical characteristics of the negative and, in such a 
case, the point O’ will not coincide with the intersection of the axis 
of the lens and the focal plane. 

The equivalent focal length is usually defined in textbooks on optics 
in terms of paraxial imagery, and it is treated as though it were a 
uniquely defined quantity associated with the lens. Actually, when 
the physical nature of light and other factors, such as the graininess 
of the photographic emulsion, are taken into consideration there is a 
considerable range of positions in which the image plane may be 
established without marked variation in definition. For each of 
these settings there is a different value of the f of eq 2, and this is 
but one reason why the f of eq 2 is not necessarily identical with the 
computed focal length. When distortion is present, eq 2 is not 
valid for all parts of the field unless it is modified by the addition of a 
correction term. If the plotting machine has no provision for the 
introduction of such a correction term, a partial correction can be 
secured by choosing a modified value of f in a manner that is more 
fully considered in a later section. In view of differences of this 
nature between the f of eq 2 and the equivalent focal length, the 
former will be referred to as the calibrated focal length, and it is 
defined as the value that best satisfies eq 2 for all points of the 
image. 

According to this nomenclature, therefore, the equivalent focal 
length is an optical characteristic of a lens defined in terms of paraxial 
imagery; While the calibrated focal length is a property of the lens 
and camera, which not only depends upon the equivalent focal length 
of the lens and the amount of distortion present but is also governed 
by the manner in which the lens is adjusted in the camera. The 
principal distance is a property of the finished negative or print and 
differs from the calibrated focal length in that it is further affected 
by film shrinkage, by the shrinkage of the paper base, if a print has 
been made, and by any enlargement or reduction. If, for example, a 
photograph is enlarged any number of diameters, the principal distance 
is proportionately increased. 


3. SOURCES OF ERROR IN THE MAP 


Errors in the map may result from errors introduced into the photo- 
graph by the camera, from changes such as shrinkage of the film or 
paper, and from errors introduced during the process of making the 
map from the photograph. This last classification includes faulty 
enlargement ; imperfect rectification or tilt elimination by projection; 
‘ystematic errors introduced by the plotting machine as a consequence 
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of inaccurate gears, pantograph linkages, and projectors: erro 
systematic for each print but changing from print to print because 
of incorrect adjustment of photographs in the plotting machine. 
errors arising from lost motion in the plotting machine; and a¢¢j. 
dental errors arising from the inability of the operator to make 
pointings unaffected by accidental variations. Errors of a systematic 
nature may also arise because of errors in the system of control points 

This presentation is restricted to a consideration of the errors that 
have their origin in the camera. They are: a, a difference between the 
assumed and actual value of the calibrated focal length; b, distortion 
produced by the optical system; c, a difference between the assumed 
and actual position of the principal point; d, a departure of the film 
from a plane; e, distortion introduced by the shutter; and f, distortion 
introduced by filter or pressure plate. 


4. GENERAL CONSIDERATION OF ERRORS OF THE CAMERA 


In evaluating the effects of these sources of error upon the map it 
is necessary to consider the manner in which the photographs are used, 
It might at first appear necessary to consider only the displacement or 
apparent displacement of an image point as produced by those varia- 
tions and then to determine the corresponding displacement upon the 
map. This procedure might be followed if the photographs could be 
used for map making without known control points. In airplane 
mapping, however, control points are necessary,’ and consequently the 
introduction of residual errors in the final map is too subtle to be 
evaluated so easily by such a simple direct computation which would 


give estimates larger than the errors actually introduced. 
In the interpretation of aerial photographs one makes direct use of 
the known control points in the picture. If (1) an automatic-plotting 


machine is used to interpret the negative, one pe the negative to fit 


the control points. In this operation it is usually assumed that there 
is no distortion and that the camera is correctly calibrated and colli- 
mated. Implicitly, this adjustment of the negative in the machine 
amounts to a determination of the position of the camera in space at 
the instant of exposure. If (2) the negative is to be interpreted by 
direct measurements and mathematical analysis, one first determines 
explicitly the position and tilt of the camera at the instant of exposure. 
In case (1).or case (2) the choice is made in such a manner that eq 1, with 
the assumed calibration, is valid, or made as nearly valid as possible 
for all control points. If these assumptions are incorrect an incorrect 
position for the camera at the time of exposure will be determined but, 
as has been noted, the error will be such that eq 1, with the incorrect 
values that have been assumed, applies to the control points. Con- 
sequently, this error in the location of the camera is of a compensatory 
nature and tends to correct the errors for other points on the negative. 
3 One may be actually fitting a photograph to control points or, when ‘‘bridging’’ from one set of control 
points to another, there may be no control points on a given print, but it may be adjusted to have the proper 
relation to the preceding photograph of the series. For the purposes of this discussion, however, these t¥0 
operations are essentially similar in character, whether the print is fitted to control points independently sur 
veyed or to the points of the preceding print. e 
‘This compensating effect occurs in all applications of airplane photography to mapping, and it consider: 
ably increases the allowable variations in the performance of the camera. It isinteresting to note, that on the 
contrary when mapping is done by cameras on tripods at the ends of measured base lines the positions of the 


cameras are observed by direct surveying methods rather than by fitting the negatives to control points and 
therefore this compensation is not obtained. 
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The same idea may be expressed, from another viewpoint, by stating 
that the camera serves as a means for determining the positions of 

oints in space by interpolating or extrapolating from the known con- 
trol points. If control points lie at the extreme positions of the 
photographs and at maximum and minimum elevations, then all 
points are “‘bracketed”’ by the control points and the process is 
entirely interpolation. If this is not the case, some of the points to 
be located lie beyond the control points and are located by extrapola- 
tion. The exact functional relationship by which the interpolation 
or extrapolation is to be done is known and the only errors incidental 
to the process are the accidental errors arising from lack of precision 
in the different steps of the operation and the systematic errors (sys- 
tematic in the results obtained from a given pair of photographs) 
arising from the use of incorrect constants, such as calibrated focal 
length and incorrectly located principal point. Consequently, the 
results obtained by extrapolation, although probably less accurate 
than those obtained by interpolation, are not so markedly inferior as 
may be the case when extrapolation is done by an entirely empirical 
method. 

The distribution of the control points is governed by the nature of 
the region to be mapped and the results available from preceding 
surveys. Consequently, the distribution necessarily varies from pho- 
tograph to photograph, with consequent changes in the errors resulting 
from departures of the camera from its calibrated values. In this dis- 
cussion of errors the assumption is made that all the control points lie 
ina datum plane. The systematic errors arising from errors in the 
camera are ea independent of the distribution of the control 
points, but for points located from a given pair of photographs the 
systematic errors resulting from the necessarily imperfect fitting of 
the photographs to the control points increase rapidly if the control 
points are not distributed over the entire region included in the photo- 
graphs. If all control points lie in a plane, the determination of all 
devations is a result of extrapolation, and consequently the errors in 
stereoscopic parallax resulting from errors in the camera may be 
expected to be maximal. If the control points are distributed not 
only in the horizontal datum plane but also over the range of elevation 
presented by the terrain, it is probable that the errors of stereoscopic 
parallax will be less than those embodied in the conclusions of this 
presentation. 


II. CALIBRATED FOCAL LENGTH 


1. PERMISSIBLE RANGE IN SELECTION OF CALIBRATED FOCAL 
LENGTH 


Ithas already been suggested that the value of fin eq 2 that makes 
this equation most nearly valid for all parts of the negative be desig- 
tated as the “calibrated focal length” and that this can be expected 
to differ from the equivalent focal length. For lenses that are not 
assembled in cameras when submitted to the National Bureau of 
Standards for test, the reports give values of the equivalent focal 
length and of the back focal length. If distortion were absent and if, 
ight consisted of rays that could be represented by geometric lines, 
uese values would each be unique. Actually, however, there is con- 
siderable depth of focus, and the values given in a report are self- 
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consistent values lying near the middle of a range of two or more 
tenths of a millimeter, the total variation without marked change in 
definition depending upon the lens and the relative aperture. If aj 
lens that has been tested is mounted in a camera, and if a depth gage 
is used to make the distance from the vertex of the last surface of the 
lens to the plane of the emulsion, equal to the reported back foca] 
length, then, in the absence of distortion or tipping of the lens, the 
calibrated focal length should be equal to the equivalent focal length 
given in the report. If, however, the distance of the lens from the 
focal plane is selected by trial exposures made as this distance jg 
varied, it will be fortuitous if the adopted distance from lens to fijm 
is the same as the back focal Jength given in the report. Consequently 
even in the absence of distortion, the calibrated focal length of a lens 
mounted in a camera may not be the same as the equivalent focal 
length as measured on the same lens submitted unmounted. 

A portion of this depth of focus, which results in an image equally 
good along a considerable range of the optic axis, arises because of the 
diffraction resulting from the wave character of light. This contriby- 
tion to the depth of focus varies as the square of the ‘“f number” 
and, for an ideal lens, has the values given in table 1. . 


TABLE 1.—Depth of focus as a function of aperture 





Relative aper- Depth of focus 





mm 
+9. 004 
. 008 

. 016 

. 032 

. 064 

. 128 

. 256 

. 512 














These limits are considerably extended in actual practice because of 
the graininess of the photographic emulsion, which limits the resolving 
power of the film. For testing purposes at the National Bureau of 
Standards, a special spectrographic plate with an emulsion designated 
by the Eastman Kodak Co. as type V is used. This emulsion has a 
very fine grain with high resolving power but is too slow for use in the 
field. Even with this emulsion, exposures made with a relative aper- 
ture of f/8 and spaced 0.3 mm along the axis show scarcely distinguish- 
able differences in sharpness between consecutive positions. This 
corresponds to a possible variation of +0.15 mm in the value of j, 
depending upon the adopted setting of the lens in the camera. 

The presence of this depth of focus is often very advantageous. 
When a lens has considerable curvature of field, it is possible to locate 
the focal plane in a compromise position for which the definition 1s 
greatly improved at the edge of the photograph without materially 
reducing the quality of definition at the center. 








2. MAP ERRORS RESULTING FROM ERROR IN DETERMINING 
CALIBRATED FOCAL LENGTH 


To determine how precisely the calibrated focal length need be 
determined and the magnitude of allowable variations arising from 
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temperature change of camera, from removal and replacement of 
magazine, etc., it is necessary to consider the effect resulting from using 
an erroneous Value of f in constructing a map. A little consideration 
pads to the conclusion, somewhat startling at first, that so long as the 
samera is vertical and the terrain a horizontal plane, an error in the 
value of f introduces no error at all in the map. All such pictures, 
taken with lenses of different focal lengths or at different heights, are 
rigorously geometrically similar, and consequently the negatives 
differ pa in scale. When the negative is adjusted in the plotting 
machine to fit the control points the scale is adjusted to the correct 
value and the representation 1s correct. 

However, if the terrain is not plane, the relief parallax ® will be 
diferent for photographs taken from different heights because, as the 
devation of the camera is changed, the displacement on the negative, 
jor the points nearer the camera, that is, for the more elevated points, 
changes disproportionately to that for the less elevated points. 
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FicurE 2.—The error arising from an incorrect value of the calibrated focal length, 


The correct centers of perspective are at Niand N3. Because of the error Af in the calibrated focal length, 
points Ny’ and Ny’ are assumed to be the centers of perspective and the elevated point O is incorrectly 
located at O’, a point displaced vertically but not horizontally. 


The effect of an incorrect value of the calibrated focal length is 


illustrated in figure 2, which is an elevation, AB being the trace of 
two overlapping photographs, 1 and 2, assumed to be so placed that the 
common control points coincide. The principal points are at P, and 


‘The relief parallax is the apparent relative lateral displacement of a point because of its elevation or 
epression from an assumed reference plane. In figure 1, the point H, is directly over A, but, in the focal 
plane, H’ is laterally displaced with respect to.4’. The displacement, H’A’, is the relief parallax and always 
extends in a radial direction from the ee. When two overlapping photographs are properly 
arranged for stereoscopic viewing, the relief parallax gives rise to the stereoscopic relief by which the eleva- 
wy of point common to the two photographs can be determined. 

"In this figure and the similar figures that follow, when reference is made to two overlapping photographs, 
this expression is used because of its brevity. Actually, if two photographs were superposed, the overlap- 
ping part of one would prevent the observation of points under the other. The superposition is accomplished 
optically in the plotting machine. The ordinary stereoscope presents a common example of such optical 
Superposition. Futhermore, if photographs are used, the distances NiP; and NP are principal distances 
lustead of calibrated focal lengths. If the photographs are contact prints from the negatives, and if there 
; to shrinkage in negative or -—~ the principal distance of the photograph and calibrated focal length 
hn camera become identical. The lengths are referred to as calibrated focal lengths because the errors 
t the camera are being considered and it seems advisable to exclude, so far as convenient, effects arising 

om the processing of the film or subsequent methods of reproduction. 


4 
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P,, and P,N, and P2N, are the two sage 0 distances, assumed e ual 
to the calibrated focal length. If the photographs are viewed shies 
scopically from N, and Nz, one has a true stereoscopic model, of which 
the scale factor is f/a, where a is the altitude of the camera at the time 
of the exposure. The point O is the correct position in the model, jn 
accordance with the scale, of a point at an elevation h above the refer. 
ence plane. Itisimaged on photographs 1 and 2 at points 0, and 9 

With the correct value of the calibrated focal length, the two photo. 
graphs, when mounted in the plotting machine, are viewed from 
points N, and N;, and O is correctly located at the intersection of 
lines N,O, and N,0,. Because of the error Af in the value of the cali. 
brated focal length, the photographs will be so mounted in the plotting 
machine that the two photographs are viewed from points N,’ and \,’ 
and the point O is incorrectly located at O’, the intersection of lines 
N,’0, and N,’0;. The displacement 00’ is equal to Ahjf/a, where 
Ah, is the consequent error in elevation of the point 0’. From similar 
triangles it is readily seen that 


ah,— "af aos. 3) 


From the similar triangles it is also apparent that 0’ is directly over 0. 
Therefore, if Al, denotes the horizontal displacement arising from the 
error in the calibrated focal length, 


Al,=0. (4) 


In eq 3, if + Af is the probable error of the calibrated focal length, 
+ Ah, will be the resulting contribution to the probable error of the 
elevation of a point on the map. 


3. CAUSES OF VARIATION IN CALIBRATED FOCAL LENGTH 


Although the calibrated focal length of a camera may be correctly 
determined for certain given conditions, it must be remembered that 
changes in the scale of the negative can arise from various causes 
such as failure of the film to always occupy the same position in the 
camera because of poorly designed or poorly constructed devices for 
removing and replacing the magazine, failure of the lens to retum 
precisely to its original position after it is removed for cleaning, expan- 
sion of the camera cone with change of temperature, and other sim- 
ilar causes. All of these can be considered and their effects evalu- 
ated as if they were changes in the calibrated focal length arising 
from changes in the dimensions of the camera. 

For a camera with a lens having a focal length of 150 mm anda 
cone made of brass or aluminum, it is interesting to note that a 
change of 50° F in temperature will introduce a change in the length 
of the cone of 0.08 mm. If the photographic film showed tne same 
percentage change in its dimensions it would largely compensate for 
the change in the cone and the camera would require little or no cor- 
rection for temperature. If the cone and film have different coefii 





imp 
the 
axis 
dist 
cen 
fowl 
que 
eq 2 


aie Camera Errors in Photogrammetric Mapping 219 
cients of expansion the resulting coefficient for the camera will depend 
upon their difference.’ ok is 

However, the great uncertainty in the determination of the tem- 
perature coefficient and its application lies in the lack of information 
js to Whether or not the camera and its spool of film have time to 
attain temperature equilibrium during the period required for ascent 
of the plane and before the exposures are begun. It 1s fortunate that 
preliminary considerations indicate that the temperature correction is 
<osmall that an approximate estimate of its value should be sufficiently 
accurate. 

III. DISTORTION 


1, VARIATION OF DISTORTION WITH SELECTION OF FOCAL 
LENGTH 


A photographie objective is designed to be symmetrical about its 
axis and, although the errors of workmanship and assembly necessarily 
imply departures from this symmetry, it is commonly assumed that 
the distortion of a photographic objective is also symmetrical about the 
axisof the lens. In such a case the displacement of any point from its 
distortion-free position is along a radius, either from or toward the 
center of the plate. In the presence of distortion no value of f can be 
found that will satisfy eq 2 for all parts of the negative. Conse- 
quently, when distortion is to be considered, it is customary to modify 
eq 2 by the addition of a correction term, writing it 


tan pu (5) 


In this equation, y is the measured distance from the principal point 
tothe image point, Ay is the linear distortion, and y—A y is the distor- 
tion-free position of the image point. The angle @ is the angular 
distance of the point from the center of the field. 

The distortion Ay is not a constant but is different for points in 
different parts of the field. If it is assumed that the distortion can be 
satisfactorily represented as symmetrical about the center of the field, 
it will be the same for all points at a given distance from the center and 
therefore will be a function of Bonly. In accordance with this assump- 
tion, in the reports of the National Bureau of Standards, in the table 
headed “Distortion” the values ot Ay are tabulated for intervals of 5°, 
measured from the center of the plate. 

Equation 5 can be written 


Ay=y—f tan B, (6) 


where f tan 6 is the distance from the center of the plate to the distor- 
ton-free position, as given by eq 2. This equation shows that the 
listortion-free position (y—Ay) is a function of the value of f that has 
been chosen and that, in consequence, y being fixed in value, Ay can 
bemodified at will by changing the value of f. After any value of f is 
chosen, the corresponding values of Ay can be obtained from eq 6, and 
‘consistent description of the image is obtained. For lenses of usual 


——— 
‘This argument is not rigorously exact because there will be a change in the position of the second nodal 
Point of the lens resulting from the change in the index of refraction of the lens components with tempera- 
However, in most, if not in all cases, this will be relatively unimportant in comparison with the 
‘pansion of the camera cone and film. 
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construction there will, in general, be a small range of values of f for 
which the values of Ay vary approximately as the cube of tan £8 over 
the field of view that is utilized. When Ay varies as tan’ 8, it follows 
that its value is relatively small for points in the neighborhood of the 
center of the field and over this region, therefore, for these values of 
eq 2 with no correction term applies with fair accuracy. ‘ 

When fis so chosen that Ay varies approximately as tan’ 8, Ay has 
the same sign over the entire field and it has its maximum value at the 
extreme edge of the field. But values of f can be chosen for which 
the distortion has one sign for values of 8 less than some given valye 
and the opposite sign for points beyond. With such an adjusted valye 
of f the distortion is distributed more uniformly, so to speak, over the 
entire negative and never equals in magnitude the value that would 
exist at the edge of the field if the equivalent focal length of geometric 
optics were selected. 

Either system of values describes the imagery correctly and will lead 
to the same final result in the map, provided no quantities are neglected 
when the computations are made. If, however, the plotting machine 
has no provision for introducing the term Ay, or if the distortion is only 
partially compensated by the machine, the two systems of values may 
be expected to lead to different results. To determine the proper 
value of f to be selected as the calibrated focal length, it is necessary to 
consider the method by which a map is constructed. 


(a) CRITERION FOR SELECTING THE CALIBRATED FOCAL LENGTH 


When adjusting the negative in a machine an attempt is made to fit 


certain points on the negative to the corresponding control points 
that have been carefully surveyed. In the presence of distortion an 
exact fit cannot be secured. The calibrated focal length should have 
the value that will lead to the least error in the map as a whole. 
This requirement does not lead directly to a criterion governing the 
choice of the calibrated focal length, because the distribution of the 
control points and the overlap of the photographs govern the manner 
in which the distortion introduces error into the map. Fortunately 
the distortion of the modern lens is so small that the different methods 
of adjusting the focal length lead to results which do not differ greatly. 
A simple and satisfactory basis for adjustment is afforded if the cali- 
brated focal length is chosen to make the maximum value of the distor- 
tion in the photograph as small as possible. 

This criterion can be applied by means of eq 6. If Ay of this equa- 
tion is defined as a function of 8, by differentiation and the usual 
method of determining a maximum or minimum, one can determine 
the value of f for which the maximum value of Ay within a field of 
given angular extent is a minimum. To illustrate this, let it be 
assumed that a lens covers a field having a half-angle of 35° and that 
the distortion, referred to the equivalent focal length, varies as the 
cube of the tangent of the angular distance from the center of the 
field. If Dis the maximum distortion, when referred to the equivalent 
focal length, then the calibrated focal length will be 1.06D greater than 
the equivalent focal length. The maximum absolute value of the 
distortion will be 0.25D, and this value will be attained for the values 
B=19.3° and B=35°. 
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It should perhaps be mentioned that this change from the equivalent 
to the calibrated focal length in order to vary the distribution of the 
distortion over the field has no relation to the depth of focus. In this 
respect it is quite different from an alteration of the back focal length 
because of the presence of curvature. Such an adjustment was men- 
tioned (section II-1) as possible, because the focal plane can be shifted 
over an appreciable range without seriously affecting the definition 
at the center of the field. When, however, the equivalent focal 
length is changed to compensate for distortion, the position of the 
focal plane is not altered. The change relates to computations only, 
the constant of eq 2 being changed in order to obtain a scale factor 
that better fits the photograph, but no change is made in the distance 
from camera lens to photographic film. 


2. ERROR ARISING FROM DISTORTION 


After the calibrated focal length has been adjusted to reduce the 
maximum value of the distortion as much as possible, the effect of 
the remaining distortion upon the map must be considered. If the 
maximum value of this residual distortion is Ar, this will represent 
the maximum error in the horizontal location of any point lying in 
the datum plane. The error, in general, will be less than Av, not only 
because the point in question may be located at a position on the 
photograph where the distortion does not attain this maximum value 
but also because, when the best possible fit is secured between the 
control points and the distorted system of image points on the photo- 
graph, the adjustment that gives the best fit distributes the distortion 
and yields a further compensatory action beyond that secured by 
adjusting the focal length. 

If the point under consideration does not lie in the datum plane 
the determination of the resulting error is more involved. One such 


case is represented in figure 3, which is an elevation, AB being the 
trace of the photographs 1 and 2 assumed to be so placed that the 
common control points coincide. The principal points are at P, and 


P,, and P,;N; and P,N2 are the two principal distances assumed equal 
to the calibrated focal length. If the photographs are viewed stereo- 
scopically from N, and No, one has a true stereoscopic model of which 
the scale is f/a, where a is the altitude of the camera at the time of 
exposure. The point O is the correct position, in the model, in accord- 
ance with this scale, of an elevated point that, as a consequence of 
the relief parallax, is imaged on photographs 1 and 2 at points 0, and 
0, when there is no distortion. 

Actually, because of the distortion, O is imaged at 0,’ and O,’. In 
the present illustration it is assumed that the displacement arising 
from distortion in each photograph is positive, that is, the point 0,’ 
is displaced outward from point P,, and OQ,’ is displaced outward from 
pont P,, the two displacements being Ar; and Ar, respectively. If 
0 were imaged at points 0, and OQ, on the two photographs, observa- 
tions from points N, and N, upon points O, and O, would correctly 
locate O at the intersection of lines N,O; and N02. Actually, how- 
ever, the observations are made on points O,’ and Q,’, and Q is erro- 
neously located at O’. The point 0’ is at the vertex of the triangle 
0’EG, of which the base, by similar triangles, is (Ar; +Ar7.) (a—h)/a, 
and the angles at the two ends of the base may, without important 

118273—39——7 
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error, be assumed to be 90°—§, and 90°—£,. (These values diffe; 
from the true values by the angles 0,N,0,’ and O,N,0,’, angles of , 
few minutes of arc, shown much exaggerated in the drawing.) 
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Figure 3.—The error arising from distortion. 


The elevated point_O is imaged at 0;’ and 0,’ instead of at O, and O; because of the distortions Ar; and Ar. 
Consequently the point O is erroneously determined to be at O’, displaced vertically and horizontally 
from its true position. 


For greater clearness, this triangle is shown to a much larger scale, 
below the main figure. The length O’H is Ah,f/a, where Ah; is the 
error in elevation resulting from the presence of distortion. To obtain 
a convenient expression for Ah;, one writes 


EH _cot E tan fp, 
HG cotG tan B, 
EH+HG _(4r,+Ar.)(a—h) _ tan 6,+tan 6, 
HG HGa tan p; 
TT. (Ar, ae Ar) tan B; (a— h) 
HG= tan 6,+tan p, a 


* Tin —_ (Ar; + Ar.) (a — h) (6) 
Ah HG cot B tan 6,+tan B, a 


If 0,0, (greatly exaggerated in fig. 3) is neglected as small in compan- 
son with P;P;, this equation may be written 
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NP i(a—h) 
P,P, f 

In constructing figure 3, it was assumed that the distortions at 0, 
and QO, are positive in the two photographs. If it should occur that 
the distortion is positive in one photograph and negative in the other, 
the quantity (Ar;+Ar,) becomes the difference of the absolute values 
of the two distortions. In all cases the value of (Ar;+ Ar.) will be 
less than 2A4r, where Ar is the maximum value of the distortion. 
Furthermore, the ratio (a—h)/a for any point lying above the datum 
plane will not exceed 1. If the ratio N,P,/P,P; is denoted by k, eq 7 
may be written 


Aha= (Ar, + Arz). (7) 


Ahi< Petar, (8) 
This inequality will be true for all points that lie above the datum 
plane. io 
The ratio N,P,/P,P2 is equal to the ratio of camera altitude to 
stereoscopic base and is a function of the focal length of the camera 
objective and the amount of overlap. The values for a film 9 by 9 
inches and for different percentages of overlap are given in table 2. 


TaBLE 2.—Values of k (or k’), the ratio of altitude to stereoscopic base for different 
percentages of overlap 





Focal length of objective 


Percentage 
of overlap 





8 in. 
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In a flight strip the customary overlap is 60 percent. The value of 
k for an overlap of this amount or less does not exceed 1.7, 2.0, and 
2.2, for camera objectives with focal lengths of 6, 7, and 8 inches, 
respectively. This illustrates the advantages of the lens of shorter 
focal length, for a given scale factor, f/a, in minimizing the error in 
elevation that arises from the presence of distortion. 

The horizontal displacement arising from the presence of distortion 
will be denoted by Al’,. Referring to figure 3, it 1s evident that a 
vertical through O’ will always cut HG, the base of the triangle O’ EG. 


Furthermore, the length OF or OG will never exceed (a-h)Ar/a. 
Therefore, 
Al’as Ar. (9) 


In eq 8 and 9, Ar is the maximum distortion on a photograph, and, 
consequently, Ah’, and Al’, are the maximum errors resulting from 
the distortion. It is convenient to have these equations in a form 
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suitable for dealing statistically with the resulting errors of the map 
The distortion in a photograph assumes all values from — Ar to + 4, 
with the maximum values in two zones, one near the edge of the field. 
the other at an intermediate distance. For a given object point, the 
distortion on the two photographs will, in general, be different, because 
the images will not fall in homologous positions on the two photo. 
graphs. Given the distortion measurements for a lens, it is possible 
to compute the error of location and elevation for any given object 
point. This error will vary with the elevation of the object point as 
well as with the projected position of the point upon the datum plane 
and will range in value from the maximum positive value, throug) 
zero, to maximum values in the opposite sense, the maximum error 
occurring only when the distortion for the images of a given object 
point has its maximum value in each of the two photographs that are 
combined to produce a stereoscopic model. In a general discussion 
of this nature it is not possible to consider these errors independently, 
but they must be dealt with by a statistical method. 

For errors that are truly ‘‘accidental’’ and follow the normal fre- 
quency curve, the error whose probability is one in one hundred, is 
approximately four times the probable error. Consequently, if the 
magnitudes of the errors resulting from distortion followed such a 
distribution law the probable error, the error that is equal to or larger 
than the errors for 50 percent of the points, could be assumed, with 
satisfactory precision, to be one fourth the maximum error that has 
just been deduced. Actually, however, the errors arising from dis- 
tortion are systematic and this is evidenced by the regular distribution 
of the errors over the photograph. A typical distortion curve for a 
lens has been assumed and, for a large number of points, the resulting 
horizontal and vertical errors arising from distortion have been com- 
puted. A study of these errors indicates that 50 percent of the errors 
may be expected to be less than or equal to 0.4 of the maximum error. 
Accordingly, in order to treat the errors arising from distortion in the 
same manner as those arising from variation of calibrated focal length 
or location of principal point, the maximum error arising from distor- 
tion, multiplied by 0.4, will be considered to be a measure of the error 
arising from distortion and will be applied in the same manner as a 
probable error. This error is analogous to the probable error in 
that it is a “50-percent error’ but differs from it in that it relates to 
errors that are not “accidental” in origin. 
Accordingly, 

ka 
AhaS0.8—Ar (10) 


f 


Al,< 0.4Ar, (11) 
where +Ah, and +Al, are contributions made by the distortion to 
the probable errors of elevation and horizontal displacement of 8 
point on the map. 

It should be noted that eq 10 has the reciprocal of the scale factor 
f/ain the right-hand member and + Ah,, consequently, is the probable 
error of the elevation. Equation 11, on the other hand, gives probable 
error as a displacement on the photograph or on a map having the 
same scale as the photograph. These are the convenient forms of 
these equations because of the manner in which the allowable errors of 
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, map are commonly specified. If one wishes to obtain an inequality 
similar to that of eq 11 giving the error on a map having a scale 
different from that of the photograph, the right-hand member must be 
divided by the scale factor of the photograph and be multiplied by 
the scale factor of the map. 


3, VALUES OF FOCAL LENGTH AND DISTORTION REPORTED FOR 
A LENS NOT MOUNTED IN A CAMERA 


In the reports of the National Bureau of Standards on lenses sub- 
mitted for test unassembled in cameras the equivalent focal length is 
given rather than an adjusted focal length. Several advantages result 
from this. Lenses are frequently submitted for test before acceptance 
and, when the distortions of all lenses are referred to the paraxial 
value of the focal length, the values permit a more direct comparison 
of different lenses, because no question arises as to whether or not the 
difference in distortion may be the consequence of a better choice of 
focal length rather than of greater freedom from distortion. In many 
cases, When a large number of lenses have been submitted for test at 
the same time, a measurement of the distortion at a single point, 
usually 30° from the axis, has been considered sufficient to classify a 
lens. In such cases some convention, of which the choice of the 
paraxial value of the focal length seems to be the simplest, is par- 
ticularly necessary, as otherwise, by a modification of the value of 
the focal length, the distortion for the one point could be made to 
have any desired value. The measurement of the distortion in this 
manner also permits direct comparison of the values with the com- 
puted values of distortion, as given and plotted in many publications 
dealing with photographic objectives.’ 


IV. PRINCIPAL POINT OF THE FOCAL PLANE 
1. INDEX MARKS FOR PRINCIPAL POINT 


Mapping cameras utilize various devices for locating the principal 
or collimation point. When a glass pressure plate is employed, a cross 
may be etched on the side of the plate in contact with the film to 
register on the negative at the time of the exposure. When a pressure 
plate is not used, four projections in the focal plane are placed so that 
they are registered on the negative and are so located that lines 
joining the two pairs of opposite points determine the position of the 
principal point by their intersection. Sometimes the four index 
marks are registered on the film by four small optical projectors 
attached to the interior of the camera cone, this method giving sharper 
and finer lines than can be obtained by the simpler silhouette. 


2. RELATION BETWEEN THE PRINCIPAL POINT AND THE AXIS OF 
THE LENS 


It has already been mentioned that the principal point is located 
at the foot of the perpendicular from the center of perspective to the 
focal plane, and that this will not correspond with the intersection of 
the focal plane and the axis of the lens, if the lens is tipped from its 

'V. Rohr, Theorie und Geschichte des Photographischen Objektivs (Julius Springer, Berlin, 1899). 


Handbuch der Wissenschaftlichen und Angewandten Photographie, vol. I, Das Photographische Objektiv 
(Julius Springer, Wien, 1932). 
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normal position. The errors resulting from such a tipping of the lens 
are not as serious as is perhaps commonly considered. In most cases 
the only error introduced by such a maladjustment arises from the 
possible blurring of the image because the plane of best focus js 
inclined with respect to the plane of the film. If there is a great dea] 
of distortion a further error may be introduced because the distortion 
will not be symmetrical about the principal point. 


3. ERROR ARISING FROM INCORRECT LOCATION OF THE 
PRINCIPAL POINT 


When adjusting ni. ga in a plotting machine an error in the 
location of the principal point, as registered on the negative, causes 
the negative to be located in the machine in a position displaced in 
its own plane from the true position. This is illustrated in figure 4, 
The drawing on the left ; 

is assumed to illustrate of oc 

the actual condition at 
the time of the expo- f 
sure. The true princi- —-+ 
pal point is at O’. The 
image of O, the point 
directly below the cam- 
era, coincides with the 
principal point. If the 
camera is assumed to 
be incorrectly cali- 
brated, so that C’ is in- 
dicated as the principal 
point, the heal. adjust- 
ment in the plotting 
machine leads to the im- 
plicit conclusion that 
the camera was verti- 
cally over the point C 
at the time of the ex- 
posure. This corres- ri 
_— toa erent ran iO” C 
of the camera irom its o 

actual position through h—ac f 4 

the horizontal distance Figure 4.—The error in the assumed position of the 
Aca/f, where Ac is equal camera arising from an incorrect position of the 
to the distance O’C’, ?7"P pee. 

the error in the location The exposure was made with the camera above Q, but it is assumed 
of the principal point. to have been above C because of the error in the indicated location 


As long as the camera is "™® Principal Po task 

vertical and the terrain plane and horizontal, the photographs taken 
with the camera in the two positions will be identical except for the 
displacement of each point of the picture through the distance C0’, 
a characteristic that has no bearing upon the accuracy of the map. 
Since the two negatives are identical the control points, if lying in the 
reference plane, can be fitted accurately when 0’, C’, or any other 
point is assumed to be the principal point. Therefore, just as a 
error in the focal length introduces no error when the terrain 1s plane 
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and horizontal, so an error in the location of the principal point intro- 
duces no error under similar simplified assumptions. 


(s) CASE IN WHICH THE TWO OVERLAPPING PHOTOGRAPHS ARE FROM A SINGLE 
FLIGHT STRIP 


The region to be mapped usually has depressions and elevations 
which produce relief parallax. In order to determine the error that 
results from an incorrectly located principal point, it is necessary to 
consider the interpretation made from two photographs that serve 
as members of a stereoscopic pair when the map is being made. If 
these photographs are successive photographs from a flight, it may 
be assumed, without serious error, that the camera underwent no 
rotation in azimuth between the two exposures. Consequently, not 
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Figure 5.—The error arising from an incorrectly located principal point. 


When the two overlapping photographs are selected from a single flight-strip, the displacements of the 
principal points, denoted by Ac, are similarly oriented in the two photographs and the elevated point at O 
is incorrectly located at O’, a point displaced horizontally but having the correct elevation. Diagrams m 
and n are projections. A projected point that corresponds to two or more original points is designated 
by the symbols of the several points, separated by commas. 


only the magnitude of the displacement of the principal point, but 
also the azimuth of a line joining the true principal point with the 
assumed principal point will have the same value for the two exposures. 
This is illustrated in figure 5 (1) which is a diagrammatic representa- 
tion of the two photographs with the overlapping portions superposed 
to bring the common control points (in the datum plane) into co- 
incidence. The two principal points are at P, and P2, the assumed 
incorrect principal points are at P’, and P’,, The displacements, 
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Ac, may be resolved into components Ap and An, parallel and norma] 
to the stereoscopic base, P,P». 

Figure 5(m) is an elevation in which significant points are projected 
on a vertical plane parallel to the stereoscopic base. AB is the trace 
of the plane of the photographs, and lines P,N, and PN), represen 
the principal distances of the two photographs, assumed equal to the 
calibrated focal length. If the photographs are viewed stereo. 
scopically from N; and N, one has a true stereoscopic model, of which 
the scale is f/a, where a is the altitude of the camera at the time of 
exposure. The point OQ is the position, in the model, in accordance 
with this scale, of an elevated point that, because of relief parallax, js 
imaged at O, and Q, in the two photographs. 

When these two photographs are mounted in a mapping machine. 
if the correct positions of the principal points are utilized, the point 0 
will be correctly located at the intersection of the lines N,O, and 
N,0,. If the principal points are assumed to be at P,’ and P,’ the 
photographs are viewed stereoscopically from N,’ and N,’ and the 
point O will be located at O’, the intersection of lines N,’0, and 
N,/0,. The lengths N,N,’ and N,N,’ are equal to Ap, and from 
similar triangles it is seen that 


Ah,=0, 


where Af, is the error in elevation, and 
h 
Al, = qo?» 


where Al,=O0O’ is the component, parallel to the stereoscopic base, 
of the displacements of the point O on the photographs. 

In figure 5(n) a second elevation is shown in which the significant 
points are projected on a vertical plane normal to the stereoscopic 
base. With this projection the points N, and N;,’ are superposed on 
N, and N,’, respectively, and lines N,’O, and N,’O, are superposed, 
appearing as a single line. These two.superposed lines actually 


intersect at the height hfja above AB. The lengths N,N,’ and N,N,’ 
are equal to An and, from similar triangles it is evident that 


(14) 


where Al,—O0O’ is the component, normal to the stereoscopic base, 
of the displacements of the point O on the photographs. The con- 
ponent An, like Ap, produces no error in elevation. 

From eq 13 and 14 it follows that 


Al.= a (15) 
a 


where Al, may be considered to be the vector displacement of the 
point O from its true position resulting from the vector displacement 
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\c of the assumed position of the principal point from its true position. 
These are actual displacements measured on the stereoscopic model 
constructed from the photographs. 
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FicurE 6.—The error arising from an incorrectly located principal point. 


When the two overlapping photographs are selected from two adjacent flight-strips flown in opposite 
lirections, the displacements of the principal points, denoted by Ac, differ in azimuth 180°, and the elevated 
point at O is incorrectly located at O’, the elevation and the horizontal position both being inerror. Dia- 
grams m and n are projections. A projected point that corresponds to two or more original points is desig- 
nated by the symbols of the several points, separated by commas. 


()) CASE IN WHICH THE TWO OVERLAPPING PHOTOGRAPHS ARE FROM ADJACENT 
FLIGHT STRIPS 


If the two pictures under stereoscopic examination are overlapping 
photographs from two adjacent flight strips, the plane will have been 
flying, in most instances, down one strip and up the other Conse- 
sequently, the azimuths of the displacements of the principal points 
from their true positions in the two photographs will differ by 180°. 
This is illustrated in figure 6(1), which is a diagrammatic representation 
of the two photographs with the overlapping portions superposed to 
bring the common control points (assumed to lie in the datum plane) 
intocoincidence. As in figure 5, the two principal points are at P; and 
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P,, the assumed principal points at P’,; and P’,. The vector displace. 
ments of the points P’; and P’, from P, and P; are Ac and — Ac, respec. 
tively, the negative sign resulting from the 180° difference in sign yi 
the two values. The displacement of P’, may be resolved into the 
components Ap and An, parallel and normal to the stereoscopic base. 
P,P. Similarly, the displacement of P’, can be resolved into the 
components —Ap and —An. 

Figure 6(m) is an elevation, similar to that shown in figure 5(, 
The point O is correctly located at the intersection of lines N,0, and 
N,O,. When points P’; and P’, are assumed to be the principal 
points, NV’; and N’, becomes the assumed centers of perspective and 
O is located at 0’, the intersection of lines N’,O, and N’30,. Also jt 
is evident from the figure that O’ will, in general, be displaced from it: 
true position both in vertical and horizontal directions. The point 
O’ is at the upper vertex of a small triangle with the base horizont,| 
and point 0 located at its midpoint. The length of the base is 2Aph/q, 
By geometrical analysis similar to that used in connection with figure 
3, it can be shown that 


hN,P 
Ah, =2———A 
etl 


kh 
on (16) 


In this equation, Ah, is the error in elevation resulting when the 
principal points of two overlapping photographs are displaced a dis. 
tance, Ap, in opposite directions along the stereoscopic base. Als 
k’ is the ratio N,P,/P,P2, as tabulated in table, 1, and is distinguished 
from k of eq 10 because the customary overlap for successive flight 
strips is 30 percent instead of 60 percent, as in a given strip. Values 
of k’ are 0.95, 1.11, and 1.27, respectively, for focal lengths of 6, 7, 
and 8 inches. 

In the triangle EO’F the horizontal displacement of the point (/ 
will be less than one-half the base, provided O lies vertically over a 
point falling between P; and P,. Consequently, one may write 


h ‘ 
Al, 7 Ap, (17) 


where Al, represents the actual displacement of 0’ from O, as meas- 
ured on the photograph. 

Figure 6(n) is an elevation in which the significant points are pr- 
jected on a plane normal to the stereoscopic base. Unlike figure 5(n), 
the components An and —An are here turned in opposite directions. 


The lines N,’O, and N,’0, do not intersect at all (their projections 1n- 
tersect in figure 6(n), but at this point the lines themselves are separated 


by the distance 0,0,). One line lies on one side of the vertical plane 


containing the stereoscopic base (the trace of this plane is N,0,) and 
the other line lies on the other side. The point of nearest approach 
for these two lines has approximately the same elevation as the poi 
O’ in figure 6(m). It is evident from this figure, where An and —4 
are shown greatly exaggerated, that strict stereoscopic fusion cannot 
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be attained for this condition, but it is attainable in satisfactory de- 

ee even though the two fields in the viewing instrument will be at 
slightly different heights. Reference to figure 6(n) shows that neither 
line departs from the correct position for the intersection more than 
jn h/a, and it appears reasonable to assume that the horizontal posi- 
tion of the point will be obtained with a horizontal error not exceeding 
this amount. The components An and —An do not introduce any 
simnificant error in the elevation. 


(c) EQUATIONS GENERALIZED TQ APPLY TO THE TWO PRECEDING CASES (a) AND (b) 


A consideration of eq 12 and 15 of case (a) and eq 16 and 17 of 
case (b) shows that one may write 
Ahess2 7 Ae (18) 


Al< "ae. (19) 


These equations are applicable to either case (a) or case (b) and cover 
all errors arising from incorrectly located principal points. In eq 
18 and 19, Ah, and Al, are the errors of elevation and horizontal dis- 
placement corresponding to the error Ac in the location of the prin- 
cipal point. The constant k’ is taken from table 1, the overlap being 
that for adjacent flight strips. As for the previous pairs of equations, 
eq 18 gives the actual error in elevation and eq 19 gives the displace- 
ment as measured on the photograph. 

In eq 18 and 19, if + Ac is the probable error of the location of the 
principal point, + Ah, will be the resulting contribution to the prob- 
able errors of the elevation and +Ac the contribution to the probable 
error of the location of a point on the photograph. 


V. DEPARTURE OF FILM FROM THE FOCAL PLANE 


Photographic film, rather than plates, is commonly used for photo- 
grammetric purposes and the assumption, in all the foregoing discus- 
sion, has been that the sensitive photographic emulsion lies in a plane. 
A suction or pressure back, or a pressure plate is commonly used to 
flatten the film. If this does not operate perfectly, it is to be expected 
that the film will be concave toward the objective, and it seems 
probable that the curvature will be different in the direction in which 
the film is rolled and in a direction at right angles to it. If a portion 
of the film is in advance of or behind the position that it should occupy, 
the image falling upon this portion of the film will be displaced the 
distance Ad tan 6 toward or from the center of the negative, when Ad 
is assumed to be the displacement of the film. If the film is curved 
in an approximately symmetrical manner from the focal plane there 
will be a systematic displacement over the field somewhat resembling 
that arising from distortion, and it is a reasonable assumption that 
this will be a maximum at the edge of the field. If the curvature of 
the film were quantitatively known and sufficiently uniform from 
exposure to exposure, a change in the value of the calibrated focal 
length could be made in order to secure compensation, as was done for 
distortion. For distortion, with a field of view of 60°, the correction 
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to the calibrated focal length is approximately equal to the distortion 
at the edge of the field referred to the paraxial focal length, and, gs , 
first assumption, it appears that a somewhat similar relation migh; 
exist between a regular departure of the emulsion from a plane anq 
the change in the calibrated focal length to compensate for it. [jf 
therefore, it is assumed that the film curves from the nominal plane 
surface by 0.1 mm, the resulting error is approximately the same as 
would result from an error of 0.1 mm in the calibrated focal length, 
but it is not the same in all parts of the field and cannot be readily 
compensated. 

By decreasing focal length and film size one can increase the ex- 

ectancy of obtaining planeness of film, but unfortunately the error 
introduced by a given departure from planeness increases with de- 
crease of focal length, so that exact information regarding the manner 
in which the departure from planeness diminishes with size of the 
film is necessary before it can be said whether or not reduction of focal 
length is advantageous in respect to the error arising from curvature 
of film. It is important to note that when a 9- by 9-inch film is re- 
tained and the angular field increased by the use of a lens of shorter 
focal length, the planeness of the film must be improved if the error 
from this source is not to increase. 

All of these considerations indicate that the planeness of the film 
at the time of exposure is an important requirement of the precision 
camera and that any increase in the accuracy of the precision camera 
may require a more complicated or better designed device for holding 
the film flat. This immediately suggests the use of plates rather than 
films in the photogrammetric camera. However, plates as commer- 
cially manufactured are not plane, and, in many cases, the holders are 
of such a nature that the plate is bent by the springs that serve to 
retain them in the focal plane. Arnulf and Perrin® have made 
measurements on photographic plates and report departures from 
planeness as great as 0.26 mm when the plate is in the plateholder of 
an airplane camera. P. Lafouasse " presents data indicating that the 
films measured by him were somewhat more nearly plane than the 
photographic plates that he measured. It is probably advisable that 
additional quantitative data be secured regarding the planeness of 
films and plates in the types of airplane cameras that are used in this 
country, and the curvature of the film or plate may be one of the 
important factors setting a limit to the precision that can be obtained 
with the photogrammetric cameras now in use. 


VI. TYPES OF SHUTTERS 
1. BETWEEN-THE-LENS SHUTTER 


The type of shutter having its leaves mounted within the lens and 
between the components, is the one most favored for airplane photog- 
raphy. With this type of shutter all portions of the film are exposed 
at the same time. Consequently, if the duration of the exposure 1s 
too great, with respect to the translatory or angular velocity of the 
plane, the resulting negative may lack definition because of the move- 
ment of the plane, but in no case will there be distortion analogous 
to that to be considered later in connection with the focal-plane shutter. 


® Arch. Int. Photogrammetry 7, 2d part, 202 (1930-31). 
10 Bul. soc. frang. phot. 24, 29-35(1937). 
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The fundamental requirement for a between-the-lens shutter is 
that it should open as rapidly as possible to its maximum opening, 
remain open for a given period of time, and then close in the shortest 
ossible time. An ideal shutter would consume zero time in its open- 
ing and closing. In such a case the exposure may be considered equal 
to the product of the area of the shutter aperture by the elapsed time 
‘rom the opening to the closing of the shutter. During the actual 
performance of a shutter one does not have the aperture of the shutter 
completely open during the entire period that elapses between the 
beginning of the opening and the end of the closing of the shutter. 
Consequently, one obtains the effective exposure by integrating the 
area of the varying shutter opening with respect to the time, the 
integration being extended to cover the over-all length of the exposure. 
The ratio of the effective exposure to the ideal exposure is defined as 
the efficiency of the shutter. It is evident that the efficiency of the 
shutter, in general, will decrease as the exposure is shortened because 
a proportionately greater time is consumed in opening and closing. 
It is also evident that the efficiency of a given shutter will increase 
as a lens is stopped down because, with the smaller diaphragm, the 
shutter requires less time to open to the full extent of the smaller 
diaphragm. ‘The motion of the plane sets a limit to the maximum 
duration of an exposure that can be made without blurring, and a 
high efficiency is important because it permits more light to reach the 
film for an exposure of given duration. Lutz and Doyle" describe a 
shutter with an efficiency of 93 percent for a one-fiftieth second exposure 
and 80 percent for an exposure of one one-hundred-fiftieth second. 

It is important that successive exposures for a given setting of the 
speed control should be identical, and desirable that the actual 
speeds do not aepart greatly from the nominal speeds. For a precise 
camera the shutter should be so designed that the lenses can be 
removed for cleaning or disassembly of the shutter and replaced 
without necessitating a recalibration of the camera. Some users of 
cameras desire a shutter sufficiently simple and robust to be repaired 
in the field and others consider that a shutter should give the minimum 
of trouble, but that it is desirable to permit no alteration or adjust- 
ment of the shutter except when it is returned to the home office. 


2. FOCAL-PLANE SHUTTER 
(a) EFFICIENCY 


The focal-plane shutter consists of an opaque blind or curtain con- 
taining a transverse slit that travels across the face of the photographic 
film, approximately in the focal plane. The light reaches the film 
through the slit and the exposure can be varied by altering the velocity 
of travel or varying the width of the slit. Adhering to the usual as- 
sumptions of geometric optics, the cone of rays proceeding from an 
object point has its vertex at the film. Therefore, a section of the 
light cone lying in the focal plane shrinks to a point, and if the curtain 
cuts the cone in this plane, the cone is entirely intercepted by the 
curtain or entirely transmitted by the slit at any given instant. Con- 
sequently, any small area of the film is illuminated by the entire aper- 
ture of the lens during the entire period of the exposure and no time 


" Photogrammetric Engineering 3, No. 3, 1-6 (1937). 
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is lost by the process of opening or closing the shutter. Therefore 
if the traveling slit could be placed exactly in the focal plane the eff, 
ciency of the shutter would be 100 percent and this type of shutter 
would be realized in its ideal form. 

Obviously, however, it is impracticable to place the curtain of the 
shutter in the focal plane and it is not unusual to find it a centimete; 
or more from the surface of the film. In this case the section of the 
cone of light in the plane of the curtain is a circle of appreciable diame. 
ter and the aperture of the lens is not entirely unobstructed until the 
curtain has travelled a distance equal to the diameter of this section, 
Consequently, the efficiency of the shutter in its practicable form js 
not 100 percent and, in important cases, may be quite as low as for 
the between-the-lens type of shutter. The efficiency of the shutter 
decreases as the section of the cone of light in the plane of the curtain 
increases, and this increases as the aperture of the lens is increased or as 
the distance from focal plane to plane of shutter curtain increases, For 
the shorter exposures, for which the slit width is narrow, the time re- 
quired for uncovering the entire aperture of the lens is greater with 
respect to the total time of exposure than for a wider slit. Conse- 
quently, the departure from ideal efficiency is greater for the narrower 
slit widths. For typical values, the efficiency of the focal plane shutter 
as a function of the slit width, lens aperture, and distance of curtain 
from film is given in table 3. 


TABLE 3.—Efficiency of focal-plane shutter 





Width of slit, in milli- | Width of slit, in milli- | Width of slit, in milli. 
meters, with //4 aper- meters, with //8 aper- meters, with //16 aper- 

Distance from curtain to| ture ture ture 
focal plane 





10 20 10 30 





% % % % 
89 89 96 

80 : 80 92 

73 73 89 


























A focal-plane shutter with a 20-mm slit, 20 mm from the focal 
plane, may be considered as typical of modern practice for shutters 
on airplane cameras. For such a shutter the efficiencies for lens 
apertures of {/4, f/8, and {/16 are shown in table 3 to be 80, &89, and 94 
percent, respectively. 


(b) DISTORTION INTRODUCED BY FOCAL-PLANE SHUTTER 


In the preceding section the efficiency has been expressed as the 
r y 


ratio of the actual exposure to the exposure obtained with an ideal 
shutter having the same over-all duration of exposure for a given 
point on the photographic film. This over-all duration of the exposure, 
for a given movement of the airplane, has a maximum value that can- 
not be exceeded without the presence of blur and loss of definition 
resulting from the movement of the plane. With the focal-plane 
shutter, unlike the between-the-lens shutter, different portions of the 
film are exposed successively as the slit travels across the field of the 
lens. Consequently, although the exposure for a restricted area of the 
film is sufficiently short to “stop” the motion and eliminate blur, it 1s 
still possible that the plane may move appreciably before the time 
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required for exposing the entire area of the film has elapsed. If this 
js the case, different parts of the photograph are made from different 
positions in ae and the photograph is distorted. 

The typical values for a focal plane shutter that were used in the 
receding section will again be considered. If the slit is 20 mm wide 
and the exposure is 1/75th second, the curtain travels with a speed 
of 1,500 mm per second. The direction of travel is the short dimen- 
sion of the photograph, and it therefore requires somewhat more than 
a tenth of a second for the curtain to travel 175 mm (7 in.), the width 
across the picture. If the plane is traveling at a speed of 200 miles per 
hour, it travels nearly 30 feet in a tenth of a second. Consequently, 
the extreme edges of the film are exposed from positions 30 feet apart. 
There is no true center of perspective for a picture taken under these 
conditions, because the portion of the negative exposed by the initial 
movement of the shutter curtain has a center of perspective 30 feet 
distant from that corresponding to the portion of the negative exposed 
as the shutter slit left the field of view. If the focal length of the lens 
is 6 inches and the altitude of the plane is 15,000 feet, this corresponds 
io a progressive displacement of the principal point of 0.3 mm on the 
photograph. Even if the curtain speed is increased so that the same 
slit width gives a speed of 1/225th second, the displacement of the 
principal point for different parts of the negative is 0.1 mm. 

A similar displacement of the principal point will also arise from the 
angular velocity (rolling or pitching) of the airplane. Estimates, 
based on several intervals of 30 seconds each in smooth air, indicate a 
value of %° per second (0.002 radian per second) as the maximum 
value that was encountered, although large variations from this value 
may be anticipated for different types of planes and for different pilots. 
With this value of the angular velocity, if the elapsed time for an 
exposure is 0.1 second, the airplane will swing through an angle of 
0.0002 radian during its progress. The corresponding displacement 
of the principal point for a lens with a focal length of 150 mm is 0.03 
mm. This indicates that the distortion arising from the angular 
velocity of the plane will not be as great as that arising from the veloc- 
ity of translation when one is flying over smooth and uniform terrain. 
When an airplane is flying over nonuniform territory, where upward 
air currents and bumpiness result, angular accelerations exceeding 
these values are athe a and may introduce serious and unpredictable 
amounts of distortion in the photographs. 

If the direction of travel of the curtain shutter is parallel to the 
direction in which the airplane moves, it will introduce a distortion 
consisting of a lengthening or shortening of the picture in the direction 
of travel. If the curtain travels in a direction at right angles to the 
motion of the airplane, the distortion will be of a different character 
in which a square is photographed as a rhombus. 

The large ratio of the time that elapses during the exposure of the 
entire picture to the length of exposure for a given small area of the 
picture is the fundamental disadvantage of the focal-plane shutter. 
This ratio decreases as the width of the slit is increased but, with the 
wider slit, it is necessary to correspondingly increase the velocity of 
travel of the curtain, a matter of some difficulty, if a short exposure is 
tobe secured. The ratio of total elapsed time to exposure time for a 
series of slit widths follows. In this tabulation it has been assumed 
that the film is 7 by 9 inches and that the direction of travel of the 
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slit is parallel to the short dimension of the film. It will be noted that 
even for a 100-mm (4-in.) slit the elapsed time is three times the oy. 
posure time. With a between-the-lens shutter, for all speeds and 
apertures, the value of this ratio is 1. 


Slit width Ratio 





mm 

20 

40 

60 
80 
100 | 
| 


VII. ACCESSORIES IN OPTICAL PATH OF CAMERA 


1. FILTER 
(a) TYPES OF FILTERS 


In airplane photography the exposures are commonly made through 
filters selected to absorb the light of shorter wave length in order to 
eliminate the loss of contrast produced by the blue haze that results 
from the scattering of light by the atmosphere. These filters are 
light yellow in color and, as the speeds of the available photographic 
emulsions increase, it is probable that filters of deeper color will be 
employed. (For special purposes deep-red filters and infrared sensi- 
tive film are now used.) The filters are of three materials (a) gela- 
tine film, unsupported and unprotected, and usually placed within the 
lens, between the components, (b) gelatine film mounted with canada 
balsam between two plates of plane parallel glass, and (c) filters of 
colored glass, ground and polished with the surfaces plane and parallel. 

The gelatine film is commonly mounted within the lens between the 
components in order to protect it from the effects of varying humidity. 
It is only the extreme thinness of the film that permits its use in this 
manner. Placing the filter between the components has the same 
effect, optically, as changing their spacing. Furthermore, the two 
surfaces of the filter, although approximately parallel, always depart 
greatly from planeness. Even under the most favorable conditions 
the use of a filter in this manner may be expected to affect the defini- 
tion unfavorably. There is also the further danger that moisture will 
eventually reach the filter, even when mounted between the compo- 
nents, and that the filter will develop pronounced wrinkles that may 
ruin a series of exposures before their presence is detected. The neces- 
sary removal of the lens cells from the shutter for replacement of the 
gelatine filter from time to time is an additional undesirable feature. 
In view of these disadvantages, filters of gelatine film should not be 
used for photographic mapping purposes although, for experimental 
work in connection with the selection of a filter shade for a given 
emulsion, the use of gelatine filters is a convenient and inexpensive 
procedure. 

The gelatine filter, mounted between glass plates, has one disadvan- 
tage. The balsam tends to deteriorate with age, the glass plates par- 
tially separate, and the filter assumes a nonuniform shade, or a feathery 
structure resulting from crystallization appears. 





sind Camera Errors in Photogrammetric Mapping 237 

4 filter of colored glass performs exactly as the filter composed of 
, gelatine film cemented between two plates, but offers the advantage 
that it is entirely unaffected by humidity and does not deteriorate 
with age. For some purposes for which filters are required, it is not 
possible to obtain glass of the desired color and a cemented type of 
slter must be used. For eliminating atmospheric haze, satisfactory 
ilters of colored glass are available, and this form of filter is prefer- 
able to the cemented type. 


(b) DISTORTION INTRODUCED BY FILTER 


Filters, either of colored glass or of the cemented type, are several 
millimeters thick and should never be placed between the components 
of a lens unless the manufacturer has designed the lens to be used in 
this manner. If the filter is placed between the lens and the focal 
plane, it will introduce appreciable amounts of negative distortion. 
The resulting distortion is of the same amount as that introduced by a 
pressure plate of the same thickness and this can be ascertained by 
reference to table 5, which gives the distortion at different angular 
distances from the axis for several glass thicknesses. Lenses have 
been built with a filter permanently mounted on the lens barrel, fol- 
lowing the last component with the entire system designed to form a 
unit free from distortion. If a filter of good quality is mounted in 
front of a lens it will not introduce any distortion. If the surfaces of 
the filter are plane but not parallel it will act as a weak prism and 
deviate the rays slightly before they enter the lens. The angular 
deviation will be approximately equal to half the angle between the 
two surfaces. As a result of this deviation, the camera, when in a 
true vertical position, will make a photograph as though it were tilted 
through the angle of deviation produced by the prismatic filter. 
To illustrate this, the difference in thickness at opposite edges of a 
prism is 0.5 mm for a filter 50 mm square (a relatively large deviation 
from parallelism) a photograph taken with the camera in a vertical 
position will show a tilt of approximately 15 minutes. If the resulting 
photograph is properly adjusted in the plotting machine this will 
cause no error, as the effect of this tilt arising from the filter will be 
diminated in the same manner that tilt arising from a nonvertical 
position of the camera is eliminated. If a filter is always to be used 
with a given camera, it is feasible to adjust the level with respect to 
the camera so that a picture without tilt is obtained when the level 
indicates horizontality. 


(c) EFFECT OF FILTER UPON DEFINITION AND CAMERA CALIBRATION 


It should be clearly understood that a filter, used with a photo- 
graphic lens, is just as much a part of the optical system as any com- 


ponent of the lens and that defects in its quality are ~ as detri- 


mental to good definition as are similar defects of the lens. Conse- 
quently, it is not good practice to purchase an expensive photographic 
objective with precisely polished surfaces and equip it with a filter of 
indifferent quality. An interference test, applied to the surfaces of the 
filter, should show plane surfaces of the same quality as lens surfaces, 
ind the glass should be equally free from internal heterogeneities. If 
the filter has slightly spherical surfaces, its use will change the focal 
length and therefore affect the calibration of the camera but this 
118273—39-—-g 
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defect will not be appreciably present in filters of good quality, 4 
wedge-shaped filter in front of the lens does not alter the position of the 
principal point. If a filter is of good quality and is to be used in front 
of the lens it is not necessary to calibrate the camera with the filter ie 
position, although the National Bureau of Standards does calibrate 
cameras in this manner upon request. This method of test is advisable 
when it is suspected that the filter may impair the definition. 


2. PRESSURE PLATE 


The newer types of cameras are equipped with suction or pressux 
backs to hold the film plane during the exposure. Nevertheless. 
seems advisable to discuss the action of the pressure plate of glys 
against which the film is pressed by a spring because many camens 
with this equipment are still in use. There is a loss of light amoun. 
ing to approximately 10 percent of that which would otherwise rege) 
the photographic emulsion because of reflection at the two surfaces 
of the pressure plate. Unfortunately, a portion of this “lost,” light 
will undergo multiple reflection at the two surfaces of the pressur 
plate and ultimately reach the emulsion as stray light, tending t 
reduce the contrast and detract from the brilliance of the image. I 
the pressure plate does not have surfaces of good quality, this will be 
somewhat detrimental to the quality of the image but the pressure 
plate is so close to the focal plane that its surfaces can affect definition 
relatively little as compared with the surfaces of the lens components. 

Probably the most serious objection to the pressure plate has been 
the large amount of negative distortion that it introduces. This dis- 
tortion is proportional to the thickness of the plate and increases from 
the center of the field outward. Values of the distortion for a pressure 
plate having the refractive index 1.52 are given in table 4.” 


TaBLE 4.—Linear distortion (in millimeters) produced by a plane parallel plaig 
behind the lens 





Thickness of plate (in millimeters 





Angular distance from center of field (degrees) 
3 5 10 




















The large umount of distortion introduced by the pressure plate cam 
be entirely and satisfactorily eliminated if the lens with which it i 
used has a corresponding amount of positive distortion. To design 
such a lens does not offer any great difficulty but this method of com 
pensation for a pressure plate has not been employed. Such a pr 
cedure would necessitate the production of a greater variety of lense 
for pressure plates of different thicknesses with consequent increase 
cost and the use of the suction or pressure back has become the mor 
convenient and desirable construction. 

WasHIneTon, October 24, 1938. 

13 For a 7- by 9-inch negative the angular distances to the extreme corners of the field are 55.4°, 44.0°, ~ 


35.9°, respectively, for focal lengths of 100, 150, and 200 mm. For a 9- by 9-inch negative the correspond 
angular distances are 58.2°, 47.1°, and 38.9°, respectively. 
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MUTUAL INDUCTANCE AND FORCE BETWEEN TWO 
COAXIAL HELICAL WIRES 


By Chester Snow 


ABSTRACT 


A formula is found for the mutual inductance and force between two coaxial 
helical wires which, in addition to the well-known current-sheet formula, contains 
sual! correction terms, one of which represents the axial components of current; 
one the finite diameters of the wires; and another wa, which depends upon the 
relative azimuths of the helices, arises, naturally, from the actual helical form of 
the windings. The pitch of the windings may be different in the two, but each is 
considered so small in comparison with the cylindrical radii that terms relatively 
aaller than the square of this ratio may be neglected. The number of turns is not 
necessarily an integer. 


CONTENTS 


. Introduction 
. Formal expressions for the mutual inductance and force between two 
coaxial helices 
. The principal terms ws and w’s 
’, Effect of axial current in helices and lead wires 
’, The azimuthal terms 
I. The helix equivalent to a helical wire 
II. Application to the current balance used in the National Physical 
Laboratory 
. Summary 


I. INTRODUCTION 


Aformula of precision sufficient for absolute electrical measurements 

the mutual inductance or force between two coaxial helical wires 
oes not appear to have been developed to the same degree of precision 
8 in the case of a self inductance. The principal part corresponding 
0 current sheets is well known, but even if the construction were 
perfect, there remain certain small correction terms which must be 
ound by starting with an idealization of the coils which is nearer the 
tual than a current sheet. The procedure here adopted as the most 
natural is to formulate the mutual inductance as Neumann’s double- 

e integral and to expand the integrand in a Fourier’s series as a 
unction of the difference of the angular parameters of the two helices. 
he constant term of this series gives the current sheet formula, pro- 
ded both angular and axial components of current are included in 
his term. The remainder of the series gives a correction depending 
pon azimuth of the two helices, which is relatively small for such 
losely wound coils as are used in practice. By restricting the problem 
0 cases where the windings of the one do not come too close to those 
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of the other, we avoid most of those evaluations of proximity effects, 
troublesome in deriving a formula for self inductance. 
The equations defining the first helix, h,, may be taken as the three 
equations which express the rectangular coordinates, 2,, y,, 2,, of 
any point, P;, on it, in terms of a single independent parameter. The 
most suitable line parameter in this case is the angle 6,, where z,, 
and 6, are the cylindrical coordinates of P,;. If the rectangular en 
are right-handed, a positive pitch, 2rp,, corresponds to a right-handed 
helix. If the plane z=z,, and the azimuth 6=6,, are those of its 


initial point, and r=z,,, 0=8@,, its end point, and if r, is its cylindrical 
radius, the equations of h, are 


t=, + Pi (6-94) 
Wi=n, COs B; , where 6;,5 6, 5 4,,. 


2=7; sin 6; 
Its axial length, /,, is given by 
1, =2e,—Li, =P 1 (Oe, —94,) =2 7M. 
Its element of length has the magnitude 
ds,=-yr?+p;? dé, 


and the direction cosines 


dz _ x Ms_ i sin dz __ 71 cos 4 


ds, yretp? ads = yr2+p,? ds — Vr pe 


If two helices of the type in eq 1 differ only in the values of their two 
constants, r;, and r;, they may be considered as filaments of the same 
helical wire, w,, where, for the purpose of this paper, a “helical wire” 
is defined as follows: 

To specify a helical wire, w,, of pitch 2p, and axial length /, and 
with a wire radius, p,, whose central filament has the cylindrical 
radius r,, the initial plane and azimuth of this central filament being 
z;, and 6,,, the plane and azimuth of its end point, z,, and #,,, we 
define it as the totality of all helices, h,’, represented by eq 1 with 
(z,,’ and 7,’ in place of z;,, 7:), provided that the two constants, 2,/ 


and r,’, lie in the range 
(x,,’—24,)* + (r;’—1;)?S py (9) 


All azimuthal planes cut this wire in circular sections of radius p;, 
the initial face of the wire in the plane 6=6;,, and its end face in the 
plane 6=6,,. All its filaments have the same terminal azimuths and 
the same axial length J;. Giving 6, a variation while holding |, 
constant corresponds to a rigid rotation of the wire about the z axis, 
while a variation in z, with J, constant is a translation parallel to 
that axis. 

The second helical wire, w,, coaxial with w, is specified, as above, 
with subscripts 2. It is assumed that p, and p, are both positive; 
that is, w, and w, are both right-handed helical axes or both left- 
handed, depending upon whiter the coordinate axes are chosel 
right-handed or left-handed. 
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The total “number of turns,” N, and N>, are not necessarily integers, 
js no essential simplification in the final formulas would be obtained 
by placing such a restriction upon the generality of their application. 
“The associated current sheet of a helix is defined as a circular 
cylindrical surface coaxial with it, having the same radius and end 
planes, on which the linear density of current has the angular com- 
ponent j= = 1/2rp and the axial component j,=1/27r. The mutual 
inductance between the sheets of the two helices will be denoted by 
Myysg = MT Me. 

The mutual inductance between one helix and the current sheet of 
the other is identical with that of the two sheets, and this statement 
holds for each component of m,,,, separately, the part mg being due 
to their angular components of current and m, due to their axial 
currents. The latter, m,, may also be interpreted as the mutual 
inductance between one of the sheets and any straight line which is 
a generator of the other. 

‘A definition will be given later of a helix which is equivalent to a 
helical wire, so that the associated current sheet of a helical wire is 
that of its equivalent helix. 

The equivalent helix of a wire has the same pitch, terminal end 
planes, (and hence axial length), as the central filament of the wire, 
but a slightly different cylindrical radius, depending upon the wire 
radius and the nature of the current distribution in the wire, which 
is assumed to flow everywhers in the direction of the generating 
helical filament, its magnitude being a function of r. 

Hence the axial length adopted here for that of the current sheet 
of a helical wire is in harmony with that generally accepted. All 
single-layer solenoids are, in fact, helical wires, but when N is an 
integer they have generally been treated by idealizing them as N 
equal coaxial circular turns of wire, their central planes equally spaced 
at r=2rpn, n=0, 1, 2,3, . . . . N—1, so that the distance between 
central planes of the first and last turn is 2rpN—2zp. If all these 
crcular turns are cut by a plane through the z axis and each given 
a shear, they go exactly into the “helical wire’’ here defined. The 
current sheet associated with this series of circular turns of wire is 
generally taken with length 2zp plus the axial distance 24pN—2rp 
between central planes of the first and last wire, that is, 2rpN, which 
is the length adopted here in general, although we do not restrict N 
to be integral. 


Il. FORMAL EXPRESSIONS FOR THE MUTUAL INDUCT- 
ANCE AND FORCE BETWEEN TWO COAXIAL HELICES 


The mutual inductance, m, between the helix h,, and that hy is 
here defined by the Neumann’s double-line integral 


cos (ds,,ds2) 
—2,)?+ (yo—yi)?+ (Z2.— 2)? 








m= Sia S der 


R?(y) =r,?+7r.2—2r,r2 cos Y, 
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this becomes 


si ee nos W (21,22) +1 
m=] dy) de >t 
Zi Zig V (t2—2)? + RW (21,22) 


where y(z;,72) is a function of (2,22) given by 








z 
¥ (21,22) 24], : (8 


The part of m which is due to the angular components of current in 
the two helices is that part of the integral in eq 7 attributable to the 
first part of the numerator 7,7, cos ¥/pip2; the part coming from the 
item 1 in the numerator is due to their z-components or axial com. 
ponents of current. The latter is a small quantity compared with 
the former (in general, of second order, when p,/r; and p,/rz are both 
small quantities of the first order). 
Writing eq 7 in the form 


Lo Ley 
m= | ds, | dx, f[22—21,9 (21,22) 
Ty Zig 
and changing the variable z, to z’ by the substitution x’ =z,—2, gives 


m={" dx, (" FS lx’ ,W(x1,2’ +2,)]dz’ 


4a—T 


Next, changing the variable z, to z’’ by the substitution 2’’=z,—1, 
gives 


os dx” ° Sle’ W@a—2" 2! —2" + 2a) de’ 


nd . 
Ze r°—(Lo—Tiz 


Leg 


4 dx” {f F(z’ W(tea—2z” %eae—z” +2’) ]dx’ 


J lq_—t 


rig , ” ” y) | 
a SF [x W(Leo—2 Leg — £ +2 )\dx"t 
0 
— ” - , ” ” , , 
= dz” | fz’ ,W(te—2” te—z” +2']dzx 
J 0 
lie af , ” ” , , 
— dx Siz’ WV (ta—2" ta—z" +2']dx 
J 0 0 


Log —Zix Y—(Zog—Zig) 
— {af Sle’ W(te—yta—y +2’) de’ 


Lala 


In the last double integral, the variable of integration z” has merely 
been replaced by y in order that the symbol xz” may be used in another 
sense, so that the four integrals into which m is resolved in eq 9 below 
shall each have the same designation of their two variables, and shall 
have limits of the same general form. In the last double integral we 
next change the variable y to x” by the substitution x” =y— (%.—%4) 
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Snow] 


Lig—Z ey ¥” 
and it becomes +f = ao" { Sie’, Hey —2",, 2 —2" +2’) ee’. 
" ss | 


Writing the 2” integral of this in the form 


(Poem 


m=["" aa" | fiz’, Vous, (2’s x") |dx’ 

— fae" fle’, Vere (2’s x”) |dx’ 
0 

+f as" | fie’, Vigo (2’s x”) |\dx’ 
0 


- = do" {" Sl’, Vy (2's x”) dz’, 
J0 0 


gives 


where 
Pi 


L 


1 1 a 
' tM) — ——— |y”+6.—90, —[{ 2 _— 
AC ’ L ) Aas ae 4 3 ty ( ~—r 


a’ 1 1 te = 
Sol) Se —— — fe" +p —f —f <a —& 
Veet» £ ) 6G. >) TO 9, ( — ) 


z’ 1 1 r,—2z 
gM) 4 (——— e419 —9 —( “4 “a 
Vine, 2 >a ) Po + 1g ej ( Di ) 


, ” es 1 ” Li, x 
Va4,(2 »o a a +0,—0,- (74 ) 


In what follows we shall use x (positive or negative) to represent one 
of the following four z-differences, and @, for the associated constant 
azimuth-difference. 








L=Leq—24) Lqy— Ley L4—Xe,, and t_—2s| 


(11) 
en a a a 


These four differences represent the four possible axial distances 
between a terminal plane of A; and a terminal plane of h.. The cor- 
responding azimuth differences, 0,, are those of the respective pairs of 
terminals (having the same subscripts). Equation 10 may by a 
similar notation be condensed into 


a1 (1h —*r pees 2 
V2(2, 22) Z4(— Do 2+, —1 (12) 
Where ¥e=Vors Vee Ver ANd Yy4, 


From the definitions given in eq 11 and 12 it is found by the use 
ofeq Land 2 that ¥,,,=Wee, 20d Vue = Vag: 
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Hence eq 9 shows that the mutual inductance between two helices js 
of the form 


n= w(x, — r,) — w(r,, —~2 J + w(r,—2,) -- w(a', a.) 13 
Consequently, the z-component, f, of the attraction of h» for h,, whe 
each carries unit current in the same direction, is given by 

, aes 
f= —o! (L_—14) +0! (L——Le,) — 0’ (Ly—Le,) tw’ (t4—2,), (13) 


where primes denote derivatives with respect to the r-argument and 


ir 7 12 4 
ae) dr, Pipe sitelindadlh 
3 1 2 dx, 4 
0 0 vat? (yp) 


where F is a function of y, defined by eq 6 and y, is defined by eq 12, 
It will be found that w is an even function and w’ an odd function of 
the z-argument, both vanishing with it. 

It is necessary to evalute the integral w (x) for positive and negatiye 
values of z in order to treat all possible coaxial helices. To do this 
we note that the integrand in the integral of eq 14 is an even periodic 
function of ¥, with period 27, and may therefore be developed in q 
Fourier’s cosine series. To find this series, consider the function 
k/4./1—k? cos* 6, where the modulus, k, is a positive real in the range 
O0<k<1. This function has the development 

nie — oo(k) 433%, (k) cos 2n8, 
4J1—k cos‘6 «4 n=1 


where the functions ¢,(k) are given by 
k (2 cos 2n8 —~})*, 2 cos 2nd 
dn k = ¢_ (k =— = == d= ~— — _k ee 
) a(®) T Jo y1—k’ cos’ 6 7 Jo y1—k? sin® 9 
They may be expressed in terms of the hypergeometric series 
, 1\,, 1 
jeontil nt+>5 r(n+5) 


: 1 1 ; X, : 
on(k) => ~TQ@n+1) ae F(ntgntgr2n+1, k). (1j 


aa 


l 
=» where 


, bets Bais, 1 
When k->1 every ¢,—~©, its principal part being $n log 


v 


k’ is the complementary modulus. These functions satisfy the dif 
ferential equation 


a . 
ea 6’) |= (4 n?— 1) dn (hk) , 


({s 


where ¢,'= s dn: 


The first two functions of the series, ¢) and ¢,, are expressible in terms 
of the two complete elliptic integrals of the first and second kind with 
modulus k. 
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9 PS Te he »\ 
b(t) = K(k) and d1 @=(44 )- A) ERG . 
T T T 


. (19) 


inv other function, ¢.(k), could be computed in terms of ¢o, and ¢;, 
by successive application of the recurrence relation 


. , 9 
(n re. ons 1 +(x —, Von =2n( — 1)¢, (20) 


The derivative, $n is given by 


1 _ $,,(k) =(n—})] (pe 1 on | 


When & is a function of x, 7;, and 7; defined by 


he. 4ryry 
a+ (r+72)? 

the mutual inductance of two coaxial circles of radii r; and rz (the 
distance between their planes being 7) is M=49?y7r,ro¢,(k) 

Their attraction, /’, with unit currents is F=—D,M, so that 


M rs 
Y, (k) =40°¢) (k) = a and Y2(k) =7°*k*{(k) = vr (23) 


These are convenient for finding ¢,(k) and ¢{(k), because the func- 
tions Y\(k) and Y2(k) are tabulated against k? in table 2 of the Scien- 
tific Papers of the Institute of Physical and Chemical Research 
Komagome, Hongo, Tokyo, 1927) by Nagaoka and Sakurai. 

When & is given by eq 22, the functions ¢,(k) satisfy the partial 
differential equation 


ee B un Qe 
D3+Di, + *D, +S" |(46,) =0 
1 


- 1—2a pew 
= ID: +D;}, + - ———Diy+= 72 n | (ry-4g,) , 


where a is any constant. 
There are also the integral representations 


on(k) 1 (* cos ny — ‘ on 
eas |, ar |, e*!7' J, (78) J, (728)ds, (25) 


where J, is Bessel’s function. The modulus, k, in all that follows, 
will be understood as the function of z, 7,, and r, defined by eq 22. 
unless otherwise stated. When z is replaced by 2; it will be called 
i; and ky means the modulus for z=0. 

Hence, as shown by eq 15, 


1 
Vii + RW.) 





1 @ 
=o) +2 2s on (hs) cos NYzf, (26) 
1 n= 
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so that the required Fourier’s series for the integrand of eq 14 is 


T1"8 00s ¥z+1 


PiP2 1 [rire 
aa 18 (fe k 
Vz} 2 +R} (vr) — TrrslPips ) + dol i) 





+> cos ny] 2 (dni (ky )+dusa(by)) +2¢n(k) | (27) 


Hence eq 14 may be put in the form 
w(x) = we (Zz) +w,(z) +w,(z, 6), 


wane tag VE (da, fT gy (hy) + ee (b) 
cae —a)| 61(ha) + PP 66(hs) fe (29 


Pip 


where 


tag (,02) = Si (2,0,) 


Wa (2,02) <7 i “dit, f, “an cos NW, (2),79) 
TP. 


Fir* (4y_ (he) + desa(b)) +2¢64(k) | (30 


The term w, is a al of the integral of eq 29 involving ¢;; the term 
w, the small part involving ¢p. 

To prove that the wy,2 defined by the integral of eq 29 is the w 
function which, if used in the general formula, eq 13, would give 
the my, previously defined as mutual inductance of the two associated 
current sheets, we may start with the elementary formulation which is 


: Zo Tes Tits * 
m ia d dz "do ab, PiP2 _ (— )+ ] 
3132 Ly 2 
Zi % 0 23 VG— 2)? +R? (6,—6) 
= dx dz. “ay Pipa aie (31) 
: . Qn 2 2 
4 3 q + (@2—2,)?+R?(Y) 


By transformation similar to those used in passing from eq 7 to eq 13 
(but much simpler) one finds that 


ly dx, i) ”" deaf (—2;) =0(ta—2n) 














—w(%2—Za) + w(2i2.—Ze1) —w(ta—Zn) 





fo 
th 


wh 


Nc 
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Snow) 


where 


. =f df sendnm | e—nyfla)da= | “dn ("42d (33) 


provided f is such a function that these integrals converge. This 
shows that w(x) is an even function of z, an odd function, or neither, 
according as f(z) is. Also w(0)=0, and w’(0)=0, where w’ (x)= 
D,w(r). In eq 31 f(x) is an even function of z, so w(x) is an even 
function of z, and m,, is therefore reduced to the general form of 
eq 13, where the w—function is 


t rz. Qe T)lo 
° ——cos y+ 1 
oygzorter= fod fda | SYPm 
4 ; : ™ Vzi+R?(y) 


z Qe 1112 
— jt 
= | (z—z,)dx dy pips YT? (34) 
L J, 2 Veit RW) 


Using the expansion eq 27 for this integrand, the V integral vanishes 
for all terms except n=0, so that eq 34 becomes identical with eq 29, 
thus proving that the latter corresponds to current sheets. 

It is useful to notice the two following partial differential equations 
satisfied by the components a, and w;, which are derivable from the 
definition 29 together with eq 24. 


Mor= (Di,+Dy+7Dn )oe=(Di+Di,+2Dy, 0, (35) 


1172 


(1-20, )oo=(Di+Ds.—z Dr, ” 
r rT; 


2 2 1 rr 
=(Di+Di,—2D,, oye $1 (ko) 


III. THE PRINCIPAL TERMS «, AND «»’ 


The equation 29 gives 


w(t) = nay’ (2) — YEE Lordy 
Pip: J9 
where 


oy’ (e) = 10 f 1(ky) dz, 
Pipa J 
Now, by eq 22 

ti 


dk, 


x,dx2,= —4r,r, E 
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Hence 
ats (a (k,)dx,= — 4(rirs)! 7 di (ky E 


Pipe Pipe k 
By eq 18 


i(k l d —€* Z 
oe —3db| -E oh sdk, (9) — (ia—1 ont) 


ae st 2 aD eee — E(k) 
~ Bardkl ky De »] 


~<a 1 4(r,72)} r K~I 
PiPr2 J, sinteedae i 3mpip2 & —1)(F ~)] 


minus the same function of k), which, i independent of z, may 
be discarded from the w function as it alw ays cancels from eq | 
Hence, placing 


Hence 


4 AnN, N: 


TMP, dale 


9 
“ 


we (a) = we’ (x) 


4m (2 ero tre| K—-(p—-1 (Kk E)| oi 


To find w,(x) we use eq 16. 


2 
w,(r) = Vrirs cos 2a f 


TPiP2 


Kenda _ 
ie 2 cos? 6 


bol 


z 
= 2 | cos 26d6 : Gey __ 
TPiP2J0 ov Ve +R? (28) 


2 


D- 2 
= 212 | cos 20 log ES Feo | dé. 
TPiPoJ 0 RE 


Integrating this by parts gives 





. 


w, (2) — 2(rir2 ed een 
pips Jol? (20) Vz?+R*(26) 


T 


_ayriray, ee 2 sin? 20d0 —_— 
 4rpyp2 o(1 (1—k,? cos* V1 —F cos? 6 


VM _ sin? 9(1—sin? 6)d0 __ 
™P:P2_ (gs =—sin? 0) VIPs sin? 0 
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T 


_wvVnr 72  - ? — (1—k’ sin? 6) 


TAP2 JoV1—! k? sin? 6 k 


eee, se 1 | 
- 1-H | 


_ (ri tr) kok 2B ee {x-n}, 


2TPiP2 





30 that 


where 


a en 
(1 —k,? sin? 0)-V1—F? or 6 


oSn7U 
: fit 2(u, k) “K+ (He zo 


The equation may also be put in the form 


| 2tNN2) ssa r(r;— 
w(t) = i th V2?+ (7, +72)? (K— b+ 1K nl 
For computing, it is sometimes preferable to introduce the Jacobian 
zeta-function, Z. 


k-—-ll= KZ(6, k’)  ——_—  « 


ark (Oo, k’ ) 4 
(38) 


iy’ Jka 


-1 
0<O= sin ae 


13 
Naya? + (ri +12)?(K— E) 


ny —re | xa (80, k’) +- =o, ail | (39) 


Where the plus sign belongs with positive z and the minus with 
negative x, the bracket multiplied by + vanishing when z=0. (Also 
the absolute value |r;?—7r,*| must be observed.) The zeta-function is 
expressible in terms of E(0, k’) and F(6, k’) for 


2(6, k’) = E(6, KN) 2 F(6, k’). (40) 


Using the AGM method of computation,! and finding a,, 6,, and c, 
by the formulas 


L. V. King, On the Direct Numerical Calculation of Elliptic functions aid Integrais, page 8. 
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tan (Ga1— 0) = 2" tan 6, 


and starting with the initial values 
/ 


A=1, bo=k, co=k’, and éy=sin-"9° 
then 


F (0, k’) On - ' , : 
rhe A) ats and Z(, k’)=c, sin 0;+¢. sin 6, . . . . Cy SIN @,. (4])' 


Using the initial values a4=1, b)>=k’, and co=k, 
gives 
K= 3q, and K—E=(@+2ci+4c.... 2" aye 


The formulas ‘ai are the result of successive applications of Landen’s 
transformation which increases the amplitude and decreases the 
modulus. Hence, they will be most suitable when k’ is small. 

If, however, k’ is large it is easier to work the transformation in 
the opposite direction of increasing modulus, in which case the 
formula is 

, 
KO, k’) + 2A Oy EY _ 4% (1—sin Y,)+KDYus Oe), 
2 2 2 
where 
a@=1, 6,= EF’, Co=k 


and 
sin (2¥_41— Va) = os sin y¥, (beginning with y= 64) 
and 


_o, tan (2¥.— yr) tan (2¥3— v2) tan (2¥,— 
2a ¢ Cn) = 2s “cos (2¥3— ppt 6s cos Qy,—v,) 1 14° ‘cos (2Q43— = 


2 — 1 
WF actisics +2(2"*— enya oe rae oem 











King? also gives an alternative to this with the same a,, },, and %, 
and the same recurrence formula between Yn41 and p,, as above, except 
that it starts with 


a v= 1 B= a 


2 There is evidently a misprint in King’s formula7 see eq 75 or 75 defining Bsa). The general tert 


should contain the factor 2(2'+*—1)ca¢s and not 2(2*—1)cn+s. 





Mutual Inductance of Coarial Helices 


Snow] 
in which case 


aa os a, sin vat >l(Ve, cn) | (41) 
2 


ik’ sin 2p 





To interpret the function w,(z) in eq 39 it will be found that if 


the factor N2/l,=1/2mp, is omitted and z replaced by J, in eq. 39 the 
result is the mutual inductance between the helix h, (or its sheet) 
and a circle of radius 7, in its end plane. This is found from eq 13 
by letting 0. The formulas for ws(z) and w,(r) when used in 


ey 13 or 13’ are equivalent to the formulas derived by Jones® for 
ihe mutual inductance of two sheets and the force between them. 
The principal part of the force comes from the w, function and is, 
by eq 13’, 

Sa= [00g (Lig — 2) — 0 (Ltg— Ley) ] — [wy (Leg —7t,) — 5 (Ley—Le,)]. 


By the above interpretation of the function w,(z) it is evident that 
the first bracket is N2/l, times the mutual inductance between helix h, 
or its sheet) and the initial circle of sheet 2 (in the plane T=24,). 

The second bracket is N,/l, times the mutual inductance between h, 
and the end circle of sheet 2 (in the plane z=z,,). 


IV. EFFECT OF AXIAL CURRENT IN HELICES AND 
LEAD WIRES 


From eq 29 


z 1 z 
we (2) = J . bu(lidey—— [ t1$0(k;) dx, =z, (x) 


VP179 Vr,T2J0 


Zz 
= f do (ki )aidaye 
17270 


rite 


Now, by eq 19 and 22 


a y — _4Vriry f K(k) », 
al, do(ky aida, al, kK (ki) a.dz, a A ke dk, 
_4y7 rind = (k) EF eo} 
- k ko 


Hence, we may take 


we(a) au (a) 208 BW), (42) 


dropping the term in kp, which is independent of x, and therefore 
always cancels from eq 13. 


J. V. Jones, Proc. Roy. Soc. (London) [A] 63, 198 (1898). 
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To find w,(r) we may use eq 17 


~~ 1 *) 
e T{ : — r r, 
w,(r) = a : | do(ki)dzy=* ; (s+5 (4r,re)* d. 1 


vita J P(s+i) Jo [ait (in) 


1 l+sny, |; ’ 2) 
= 3 108 {sin yt 2 7 SC, Be) 


where 
x 


r+? 


=\ po 4 1 
str, k=!) 0+ 3) eyH(1—s, 1,5 st +) 
— T?(s+1) 7 | 


y=tan7! 


and 


Now, sin yS is, in general, small compared to the logarithm, except 
when y=0, in which case they both vanish; but as z becomes large, , 
approaches 7/2, and the logarithm becomes infinite, while the series 


») 
approaches the finite limit S($ to) =log Tk’ which vanishes if /,)=0, 
= 0 


and has its greatest value, log 2=0.69, when kj=1. Also S vanishes 
when k)=0, so that, since the term sin yS(7;k) is a small correction 
to a second-order term, a good approximation is obtained (sufficient 
for all cases with which this paper is concerned) by taking 


9 
sin yS(7,ko) =sin ¥ log i TR 
0 
so that 


r+ y2?+ (itr)? x Trtt\ 149 
w(x) =log | — : it i rece log on ) (43) 


M+Pe 








, ; 2 mtr. 
where r,, is the greater of r; and r, ( since att), It may be 
1+ky Tm 


noted that the method used to obtain So oi (ky dar in the preceding 
section leads to 


to Aft, [et iF TRO) 
[jetkndn—? [ log | R() fio 


An approximation similar to this may be obtained from eq 17 in 
the form 





1 . ’ f 
a, on(ki)dz,;=sin yko"S, (y,ko) ’ 
1 





wl 


80 


He 
ob’ 


sO 
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where 





r(stnt5 (1 _ 


ee 
S,(7,ko) == SK T(s+)ll(s+2n+1) N15) >? Sin 7) 


— 


when 


sin? y=1, # i~e~maie, j= ane 
n(e+n+3) 


1 (/ko\" 1 1 : 
k2"S, (Ft0)=5,( 3) P( nn 5p2n ‘ 1,K8) 
1 kK \* 1/n\*. 
=a THR) ann) fn 


—1(n\"; 
=3,(2) if <r}. 


Hence an approximation sufficient for use in second order terms is 
obtained by placing 


so that 


2n 
* sin 1Sa(rsbo) =U) 


so that the following approximation analagous to eq. 43 is obtained: 





zt 1 , 
i (k,)dx,= V2?+ (nr)? i(2) if n>0 and To >T. (44) 


When the lead wire of length 1, for the return current of /, lies parallel 
toa generator of the current sheet but at a slightly greater distance 
from the axis, say 7,,=7:-+Ar,, the mutual inductance between this 


lead wire, 7,, and the helix, hz, is (as far as second-order terms are 
concerned) equal to that between /; and the current sheet of Az, the 
corresponding w-function being 


—w,(z, Tuy, T2) = —w,(z, ‘1, 72) — Ar,D,,w,(2, Thy 12). 


Ifthe lead wire of No. 2 is similarly situated at a distance r;,=r,+ Ar; 
its mutual inductance with h, has the w-function 


—w,(z, "1, Ti») = — (2, 1; r2) —Ar.D, 202 (2, TiP2). 


This arrangement, of course, implies that the number of turns, N, 
" d N., are integers, so that the azimuths of the two lead wires are 
4, and Diy. Their mutual inductance with each other is determined 


by the w-function, «2,2, a8 in eq 13, and their force by w;,, as in 
eq 13’, where 


Wyn (2, a, Ty, Ty) =2w},,— V2?+R?, (a) +R, (a) (45) 
118273—39--9 
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[x2 +R? (a) 4 ‘ 
1, (2s a, Ty, r,) =log Vv . a= f _ dey 
nual) 0 VET R?, (a) 





where 
a=61— 64 


RB}, (a) =77, +7, —2ryry cos a. 


The additive constant, Ri, has been omitted from eq 45, beins 
independent of z, so that it cancels from eq 13. When the radiel 
distances, Ar; and Ar., of the lead wires from the helices are i) 
the axial currents in the lead wires and helices almost compensate ir 
their effects upon inductance and force, the residual being 9 nes 
small term, depending upon the azimuth, a, of the coils. If &,(r. ; 
represents the total w-function due to z-components of current in 
lead wires and helices, then : 


@;(z, a) =w,(2, nr, 12) +on1,(2, a, T+ Ar, 72+ Ar.) 
—w,(z, 71+ Ari, 72) —wz(2, 71, T2+ Ar). 
When Ar; and Ar, are small, this becomes 
@,(2, a) =n (2, & Ty Ty) —2(2, M1, 72)— (AnD, + ArD,,)0;(2, 1, 1) 


Since, however, these w-functions are themselves of second order 
the terms in Ar are smaller than second order and therefore negligible. 
Hence, for this arrangement of lead wires, the w-function ‘and its 
derivative, which takes account of all axial components of current 
in helices and lead wires, are given by 


47m E 
T 


r —JP+Ra) (46) 


Ws (x, a)= ru, (x,a) 4- 


[21 22/,) 
@.(r,a) =log & - an (a) —w,'(r)= ) b, cos na, (46’) 


n=1 





the constant term in this Fourier’s cosine series being zero. Tlie 
coefficients 6, for n>0 are given by 
2 . ; 
= | on(ky)dx,= (46 } 
Vrir2 0 
~ x 1fn\ 

V2?+(rm+r2)?n re 
by the approximation in eq. 45, ifr2>r,. On summing this Fourier’ 
series we find the following finite form as an approximation for eq 46 





(46'"’) 





o’ (2,a) = 


a ee 
Ve+ (+n)? ©? R(a) 


when 73>7}. 
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V. THE AZIMUTHAL TERMS 


To integrate eq 30, we may place (by reference to eq 10’) 


; ; +O 
cos ny.=e™* —¢in(%-5,)+ beth -d ae 
with the understanding that the real part of the result is to be taken. 
Using the abbreviation 


Jin,s) =Fa-1(718) Jo (728) + Tuts 18) Tusa (ra8) +P, (718) Jn (728). 


eq 30 becomes, by use of eq 25, 


w,(t,0;) = "it: , ,in( a5 n) {~ dsJ(n a) fre in( 5 - me de [e pr Hl dy 


PiPo 


Ifa>0, then 2; and 22 are both positive in the range of integration 
of the 242% double integral, so that e~*!*!!' =e, When «<0, then x, 
and zp are both negative in the range of the integral, so that e~#!*! = 
e, and the result of the integration in the latter case is found by 
changing the sign of s in the result of the integration for the first case. 
it is thus found that both results are included in the following: 


inz 
po(e- os — eM) _ prile- 7 —e~#itl) 
w,(t,0:)=—a7-——5 |. dsJ/(n,s De 
" ~ { dss(n,) me. 14s 





ne 


in@z re) 
1 Po€ 
= | dsJ(n,s) 
4 0 


_int _inr 
Pse Pr pre Pi e~#l2| 


»—pn( 122) | (p— pi)( 14:2) OG) 


where the upper sign belongs with positive z and the lower with 
negative x 
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The terms in the first integral which do not contain the exponential 
factor e~*'*! may be discarded, as they will all cancel in eq 13 because, 
gs shown by reference to eq 2 and 11, the values of 6,—2/p, will be the 
same for the first and second z-differences, (xeg—zi,) and (xe.—se;). 
These first and second z-values enter eq 13 with opposite signs. 
Similarly, the values of 6,—z/p, are the same for the third and fourth 
p-values, Zg—Le and Z_2—2y, which enter eq 13 with opposite signs. 
Hence the term in eq 47, cos n(@,—2/p,), disappears by the second 
term canceling the first and the fourth canceling the third. It is 
aso found that cos n(6,—2z/p2) disappears by the fourth canceling 
the first and the third canceling the second. The corresponding 
terms in the second integral, sin n(@,—2z/p,) and sin n(0,—z/p,), do 
not always cancel because of the + factor which has the sign of z. 
If the associated current sheets are wholly external to each other the 
cancelation is complete, but if one lies wholly or partly within the 
other the terms sin n(@,—z/p,) and sin n(@,—2z/ps) of eq 47 do not 
disappear from eq 13. ; 
Hence eq 47 may be written 


@o 


etl (n)(1 —PsPays) 
a TT Ne 
w, (2,02) =— aly cos n6, pst ~io7 
. + oil nl 
0 oe X42) 


o se-"I J (n,s) 74 
<"ita(PitP2) t= Mis ds + w_*(x,02), (48) 
m 4 sin n6, (+85 \a+2s ) 
” 0 nr? n? 


ds 








plus terms which cancel in eq 13, where 


Tl 


(p2—P1) 
Hs Con) —sin mn) [" er 
2 eal 1 oe ea (e 
0 ¥ TD, 0 T3 


We may now show that by limiting the application of the results of 
this section to those cases where the o; (x) are all negligible, we do not 
impose a serious restriction as far as practical mutual inductances are 
concerned, with the exception of a bifilar-wound coil in which one of 
the windings is used as primary, the other as secondary. Such a coil 
lias never been used for absolute measurements, but we see no reason 
why it should not be used in the future. It would possess certain 
advantages from the point of view of accurate construction and 
leasurement and would give the greatest possible mutual! inductance. 
The formula for its mutual inductance has never been evaluated with 
precision as it has for a self inductance, but the problem would present 
no greater difficulty. 

It is easy to foretell that this limitation is that the helices, h, and h, 
shall not through a finite part of their lengths be separated from each 
other by a distance which is small (of the order of pir). 


w, (2,02) ae” 
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If H, denotes the first Hankel’s function then 
eee Cre mr. ; 
fi 2 Pp Pp 
p* 
when r;>7;. In case rz<7,, the two are interchanged in this equation, 


Since r,/p and r2/p are large, the asymptotic expansions of //,, and a 
show that this expression is very close to 





=nlr—nl 
oy ze =i cases ro >T, OF T2<r}. 
ony 1112 


Hence, approximately, 


n . - 
—p,lln—nil+i@—p.,)] 


pie | 


* WE nts 
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where the real part is to be taken. 
When p.=p,=—p, this becomes 


Vrir2 1 ve pliant +ie—pit)) 


rr. _— Re, Se ee 
ertoas) = VEEL A Fanti Brisco 


=" aa 

a |n np 

Since the entire w, function with the factor 7,7, is of second order in 
general, compared with the principal term, ws, which has the large 


Fils 2. 5 - a ; oe 
factor——, it is evident that p p—"! will be utterly negligible when 


1P2 I 
r2—T; 1s finite (not a small quantity of the order of p). Since o* is 
only ditferent from zero when one coil lies wholly or partly within the 
other, this means that the two must not approach very close to each 
other if wi is to be neglected. 

The remaining integrals in eq 48 cannot be similarly integrated in 
finite terms, but it is evident that when |2| is finite (not small of the 
order p/r), we may (on account of the convergence factor e-*"" of the 
integrand) obtain the principal part of w, by placing p,;=p,.=0 in 
the integrand, as the resulting integrals still converge, even if r,=", 
We cannot be content, however, with this restriction on z, as we must 
evaluate w, for cases where x=0. In this case, if we place p,;=p,=0, 
the integrals in eq 48 then converge only if 7,472, but if 7;—7, is a 
finite quantity, the principal part of w, is thus obtained. The limits- 
tion is therefore evident—the helices must be everywhere separated from 
each other by finite distances. Assuming this to be the case, and 
noting that the last term in the definition of J(n, s) is negligible, we 
obtain the result 


Wy (2,02) ae Veto Eg, 1k) +bnu (8) cos n6z, 


so that 
cos né, 


w4(2,02) = — rr, D3 [64-14 burs] 





Mutual Inductance of Coaxial Helices 


and 
tien, Me Nes , COS nbz, - 
we (2,02) — Frnt Ol (49’) 


In finding the self inductance of a helix, the wq rises in rank from a 
second order to a first order quantity, from which originates the most 
important correction term—the evaluation of which was the only 
dificult part of that problem. The same remark would apply to a 
bifilar mutual inductance). Eq 49 holds of course when r;=r2, pro- 
vided the end planes of the two coils are not close together and the 
coils wholly external to each other (so that k& is not close to unity). 
The series 49 and 49’ may be evaluated by finding ¢) and ¢, in 
terms of elliptic integrals by eq 19, and then finding all other ¢, or ¢, 
by the recurrence relation 20 or 21. 
‘Since the series represent second-order terms, it is evident that 
some simpler approximation will be sufficient. To obtain such, con- 
sider the even, continuous function of 6, {(@), which is never negative 
and is a periodic function of 0, with period 2z, defined by the series 


@ 


ye for all real values of 6. 50) 


n=1 


(§) vanishes when 6=2nz, where n is any integer. It has the maxi- 
mum value unity when 6=(2n+1)z. It consists of a succession of 
parabolic ares for 


when 056K 2z, (51a) 


when 275047. (51b) 


cos N(0+6;)_ 


n? 


Also from eq 15 one finds 
k s6 <. 53) 
3 = = ¢1(k) + Dink) + nih] cos 6. (53) 


i—k cos? 2 
< ¥1—k* cos’ 2 


Multiplying eq 52 by eq 53, and integrating the product from ¢=0 
to §=27, gives 


k (2*f(0+40, 16 F 6, 
= ‘f(8+6,) cos 64 = Fe) 2) bert bun) SE 
n=1 








2Jo0 ¥1—k? cos? 2 
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s0 that eq 49 becomes 


Qn 
J (0+ 6, 6d8 
t04(2,0;) = —V rire Tok) {(0+6,) cos 
} aft—# cos cost § 


24 2e 
~~ ae oh ~** $040, ) cos. “ -* J(8+82) cos 9 cos 6 
172 ae te ee” 
, yi —h cost —k? co ¥ 
0 2r—02 


It is sufficient to evaluate this for the range 0 =0,< 7 since its value 
for any other range of 6, is obtainable from ‘this by inspection, remen- 
bering that w, is an even function of 6, with period 27. (This restric- 
tion disappears from the end result. ) With this restriction the 
argument of f in the integral from 0 to 2r—@, lies in the range 
0<0+0, S27, so that by eq 5la 


f0+0,)=(2—-£2) (24%), 


In the integral last written, as eq 51b shows, 


f(0-+0,)=(4—2£%) (282), 


Making these substitutions, and then changing the variable of inte 
gration in the last integral to 6’, where 0=2x—0’, gives 


2n— Oz 


[2- io pett +0 cos 6dé 


. y1 —k cost 5 
rk 
—— [2-2 | = @ cos adi _| 
T cs 6 
yi — k? cos? ; 
0 2) 
or : 


2r—0, 
mon lio [ {l0.+6) cos 040k [-"f(0=—0) cos tater " 
0 


1-—k cost 5 . ye ke cost 2 


where both of the arguments 6,46 and 6,—8 lie in the positive range 
of less than 2x. For this range, instead of the cosine series 50, we find 


a sine series for f(@) 


ae 
w. (0,02) = —a/ryry 5 1(k) ae 


sin a sin — 
V4 22 


56 
4 =az<« 4 


Since 32/x*=1.032, it is evident that a good approximation for the 





integrals in eq 54 is given by placing {(6)=sin 5 for 0S 6S 2r, thatis, 
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j(0.+6)=sin (5) and se.—0) =sin (">"), 


This gives after some simple transformations 


6, 2 sin 6 cos 26 
w, (2,02) = —Jrrde 3 vile (kh) +5 me as’ A) -” "2 vi ——. do 


(3 


ak 6, (5 ) cos 6 cos 26 
+> cos 3 Vie sin? 6 dé A (56) 


a—O2 


Now, 


j (Sat ces 28 dp—7, [log (k cos 6+1—F sin’ 6) 
V1—-F s sin? @ k 
—k cos 6V1—k? sin? 6] 


and 


k sesh Bh lh =A Lip sin= (k sin 0)+k sin 6-¥1—K? sin? 6], 
Vl—k’s sin ~ 


so that 


ge Oe, al Ti? 
w,(2,9,) = — a — IN 1—k? cos? : 


2 





: sin Saf —k? cos? 3 (57) 


Oe 6, I. 0; 
—7z 00s = sin (ik cos =) — pins a P 


Differentiating this with respect to zx gives 


(2,0) = mk py (hk) 


ao 


a 1-3) 
cos & = 5 sin“ ( cos = 3) — ew )1-# cos 4 —k? cos 2 
k sin Sta /1—# cost cos? — 
+k sin % el " 67") 


k’ 


re 





These equations hold for any value of 6, since they are even periodic 
functions of 6,, with period 27. 
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The four azimuthal differences in eq 11 may be written. 


90:4, =A+20N2,Bere,=at2r(N2 —N)), 


ean ee | J =—=py= inn 
B ine, = & 24,9191, = &=01.—6;, (58) 


so that when the number of turns, N, and N2, are integers, every 9. 
may be taken equal to a. Itis evident that no essential simplificatioy 
would have resulted in the formulas eq 57 and 58 had we restricted 
N, and N, to integral values at the beginning. 

To obtain an interpretation of the term w,(z,0,), let M(x,0,,r,7,) 
or, more briefly, 1/(@,), denote the mutual inductance between the 
two incomplete circular arcs which are the parts of the two coaxial 
circles of radii r;r; (the distance between their planes being z) included 
between the azimuth planes 0=0 and @=60,. 

When the circles are complete we have by eq 23 


M (21) =42° yrirodi (k) 
In general, if 0<0,<(27 


ile ™ ” cos (6.—6;) 
M(6,) = nin, de, f dO, Vz? +R?(6,—6,) 





which may be transformed into 


MGs) _ (a (” 208 brah 
2Vrv2 Jo “Jo yi-# cos? 
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Using the expansion in eq 53 for this integrand gives 


M(6.) _& 
rir, 2 
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Equations 49 and 59 are equivalent to 
2 
024202) 0, (2,0) +2 Me.)—(5*) Mer) | when 0<6,S2r, 
Since, by eq 57, 
- — MGs) _nvreek'( se ier ‘) 
w,(2,0) = 12 of ~( 1 I: sin! k ); 


this may be written 
2 
(2,0) = — 422) 4 3| Me. —(#:) Mer) | 
ke! / 50) 
eure FE sint k 


which furnishes an interpretation for part of the term w,(x,92), although 
this form has no advantage in computation since the mutual indu- 





(60) 


ugh 
luce 
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tance \/(8;) of the two incomplete circular arcs is not a well-known 
or tabulated function as in the case M/(27) when the arcs are complete 


circles. In fact, the foregoing formulas give for it 


Se ee cal 
6\? — y! —k* cos* =—k’ 
M »-(s:) Meni 2 





k 
— ‘Tal c0s . sin! (i cos 5)—sin™ k| 


(61) 





~ ie sin > log 
iC ye 


When the series 49 and 49’ are known to converge so rapidly that 
only the first two terms need be retained, there is no need to use 
eq 57 and 57’. 


VI. THE HELIX EQUIVALENT TO A HELICAL WIRE 


Consider the unit current in wire w, to have the vector volume 
density whose magnitude is a function of the distance r; from the z 
axis, say U, (r{). Its direction is that of the generating helical 
filament. 

An axial plane cuts the wire in a circular section of radius p;, as 
specified in eq 5. If dS, is an element of area of this circular section, 
a helical tube whose (oblique) section is dS, carries the current, 


dT, =a 


Ji? 


since 


is the cosine of the angle between the normal to dS, and the direction 
of the current-density vector or tube. The total current carried by 
the wire is unity, so that 


iter 


integrated over the circular section of the wire. 

In place of the rectangular coordinates 2j,—2, and r,—r, of a point 
of this section referred to the center of the circle as origin (as used 
neq 5), we may use the polar coordinates p, and ¢;, where 


Ti, — Zi, =p, COS d, and r;—r,;=p; sin gy, 
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in terms of which 
dS, = pidpidq,. 


Expanding the current density, w(7:)=wW(ritp; sin ¢;), about the 
value r, (corresponding to the central filament of the wire) by Taylor’; 

theorem and applying the condition 62, we obtain (to the second 
order, inclusive, in p;) 


dT, =2tedd E + Ajp; sin $:+B, (v3 sin? 6—f)) 
T pi 4 
where 


2ui(r1) 


U(11) 


Aah) d B= u ‘(r;) 


With similar definitions of A, and Bz, etc., the current, dJ,, in a helical 
current tube of the second wire is 


died Pitot 1+ Aas sin ¢.+B, (0? sin? ¢.— ali 


2 
so that 


dT, I, Pee ool + (Ap sin Gr +Aaph sin ¢) 


2 
+B, (oi? sin? 6) + Ba( oi sin’ dy) + AvAap sin brphsin 61-63) 


The mutual inductance, Myi02, betw een the two helical wires, w, and w,, 
is found by multiplyi ing this expression for dI, dJ, by m and integrating 
over both circular sections of the wire, an operation which is equiva- 
lent to multiplying by 1 when applied to any terms of m except the 
finite part ms, since the operation only alters the subject by a second- 
order fraction of itself. 

Now, by eq 13, replacing the # by 2’ 


fal ffdTmo= ff dl, ffdT,00(x.,—2;,) — ff al, [fdTowe(2.,—21,) 
+ ffl, [fd Toc (a,— 21.) — [fd [fd Tae (x,—21) 


Also 
Ww (7.—%,) =we(zr' ¢ —Z,, "1, r2) 
= w4(%e,—21, + p2 COS $a— pi COB $1,71-+ pr SIN $1,72+ p2 SID dz), (65) 


where x,,— Lay) Tr; and r, refer to the central filaments of the wires. 


The expansion of this by Taylor’s theorem for a function of three 
variables may be written symbolically (to the second order in p, 
and p2) 


cn(0,—24) ={1 + [(p 008 dpi 008 $1) Dat 0} sin $D,,+ p2sin aD, 
+51(p% €08 ¢s— p 008 $s) De-+ o; sin $1D,,+p28in ¢D,P}wo(2,ri7%2), (66) 


where 2=2,,—2,. 
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Snow] 
Multiplying eq 63 by eq 66 and integrating gives 


{faz {fale 1009 (2, — 
ne (D?+ Di +2A,D, +4 (D+ Di, +2AD,,)joo(x,ri72), 


where t=d',—2%,,=2=the distance between end plane of the central 


flament of wire w, and the initial plane of the central filament of 
wires vw. On making use of the partial differential 36, this becomes 


so that 
[fat f [atwrta,—25) -colta—ra) 


+31 a(; +24, )Dnt0i( 2 +24, )Day Juv EAE Vir (ks). (67) 


Pip: 


The expansion of the three remaining integrals of eq 64 is the same as 
this with the appropriate x-arguments. 

The result is that, to the second order, the mutual inductance of 
the two wires is given by 


2 1 
a +5 oi( +24: Dnt 0i(2+24s)D., fre (68) 


where h, and hg are their central helical filaments, with radii r,; and r2. 
This may be written 


Meniwa=Me(7,72) +Me+1Ma, (69) 


where 
7, pi 1 9 d a ps rs 0 
=nt 3 rb PAL an n=rts (1+2A,), (70) 


where all the m’s except the principal one refer to the central helices 
of the two wires, and every m is expressed in terms of its w-function, 
as in eq 13, with corresponding subscripts. 

The sole effect of giving the wires a radius is to increase the effective 
radii of their current sheets from r to r (to the second order), the 
length of the sheets being unaltered. The latter fact is explained in 
part by the particular manner in which their lengths have been 
defined, and in part by the fact that the wire sections are circular, so 
that the terms which would represent an increase in effective leneth— 
being proportional to Diws—have entered the result together with 
log and their combined effects have been expressed in terms of 
D, io and D,,.we by means of the partial differential eq 36. 

For distributions of current in both wires whose r-derivative 


vanishes at the central filament, the constants A, and A, are both 
zero, 
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The uniform current distribution is an example of this class, 4 
correction equal to the above but opposite in sign is obtained for djs. 
tributions in which A;=—1/r;, and A,;=—1/r:, an example of this 
class being the “natural” distribution (B,=1/r}), where the currey; 
density varies inversely as the linear length ¥r?+p? of the helical 
current filaments. A distribution for which A,;=—1/2r,, A,=—}), 
would have no correction of this order of magnitude. 

Whatever assumption be made as to the current distribution, it wi 
still be necessary to evaluate D,,m and D,,m roughly in order to 
estimate the effects of errors in determination of the mean radij 
r,; and 7». 

At present we are ignorant of the current distribution in such wires 
so that for simplicity it would seem better to ignore altogether this 
correction for finite thickness of the wires. The assumption usually 
made of uniform current density carries with it a spurious precision, 

We can, however, allow for future increase of knowledge as to 
current distribution, with no less simplicity of treatment, by defining 
the helix which is equivalent to a helical wire as one having the same 
pitch and end planes as its central filament, but a cylindrical radiys 


2 
raat 4(5+24)p 


where a is the mean cylindrical radius of the helical wire and p is the 
wire-radius—the interpretation of the constant A being the value at 
the central filament of the r-derivative of the current density, which 
was assumed to be a function of r only and in the direction of the 
helical filament. 

As far as we are concerned here, this disposes of the difference 
between helices and helical wires. All the formulas in the preceding 
sections regarding mutual inductance of two helices and the force be- 
tween them are valid for the helical wires to which they are equivalent, 


VII. APPLICATION TO THE CURRENT BALANCE USED 
IN THE NATIONAL PHYSICAL LABORATORY 


If the x axis be taken vertically upward, the lower (or suspended 
helix, called No. 1, has a mean radius r,=10 cm and length /;=15.2 
em. The upper helix has radius r,=16 cm and length /,=11.0 em. 
Both helices have the same pitch, 27p,=22p.=0.2 cm and an integral 
number of turns, N;=76 and N.=55 turns. Hence, every 6,=a= 
6;.—6;, the azimuth of the helices. The upper end of No. 1 projects 


into No. 2, so that the third z-difference of eq 11 is negative. 
2 =L.—24, = + 20.1 cm—k,?=0.5926 
L2 =Leg—Ley = +4.9 cm—k,?=0.9143 
23 =21.—Le, = — 6.1 cem—k,?=0.8974 


4 =2.—L,= + 9.1 em—>k,?=0.8434 


The force, f, acting on the suspended helix when both carry wit 


current is 


S=SotF2(a) +ha(a), (71) 
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Snore) 
where fy is the principal part due to angular components of current in 
the helices and by eq 13’. 

fo=— es (01) +09" (12) — 00! (as) +e (20), (72) 


where the w0’(x) are given by eq 39. The small part f,(a) denotes the 
residue due to z-components of current in helices and lead wires 
together and is given by 


Fa(ce) = —%, (111,02) + 0; (222,02) — %, (213,02) + ©, (24,0), (73) 


where the w,(z,@) are given by eq 46’. 
The last part, f,(a), is given similarly by 


te (a) ae — We (x1,a) + We (22,a) ey W (x5,a) + Wa (24,@), (74) 


where the w,(z,a) may be computed by eq 49’ or 57’. Both w;(z,a) 
and w,(r,a) are even periodic functions of a, with period 27, so the 
same is true of f(a) and fa(a). 

In the actual coils, however, there are two helices on each cylinder 


differing in azimuth by z, so that when each carries unit current the 
force on the suspended coil is F’, where 


F=F,+F,(a) +F,(a), (75) 


F where 


Iy=4fo 
F(a) =2Uf2(a) +f2(a+7)] (76) 
F(a) =2[fa(a) +fa(a+m)] 


This makes the azimuthal terms, F(z) and F,(a), each even periodic 
functions of a, with period z, so that 


F,(a)-+F,(a)=> Cr, cos 2na. (77) 
n=1 


?. Vigoureux of the National Physical Laboratory of England has 
recently reported to this Bureau that in the absence of a formula for 
the coefficients C;, he has evaluated experimentally the first coefficient, 
(,, assuming the first term to be the only one of the series of practical 
importance. This he did by measuring the force for three different 
azimuths a=3°, 15°, and 85°. 

In order to see what part of C; is due to axial components of current 
and what part is due to the azimuth terms «@,%, let 


Con=CS,+C3,. (78) 


Reference to eq 72, 76, and 46’’ shows that the effect of axial currents 
is given by 


1 2n 
n= —2(*) {—ayk,+ak.—agk3+a4ky}, (79) 
2 


Vrire 
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and eq 49’ shows that 


1 1 ai , : 
[i Viits an? { —axyki[don4s (ky) + gon—1 (ky)] + rok3[ don 11(ke) -t. Pan—1 (hk) 
1/2 


—xrskildon+ (ks) + 2n-1(k3)]+ tskil bons (ks) + G2n- i(hy)] } (80) 


Placing n=1 


1 2 
Cz= *') [—ayk, + 2k.—axghs+ 24h] = +0.103 dynes (81) 


vrir\P2 
ae SF ey +p SE) FOE) 


Vinita 








From the recurrence relations 20 and 21, it is found that 


kh) +45) ]_1 (01, 8k] op) : : 
pl sOreO|-H 5 tar Tv ho, (k) ors (i- 1 tx $, (| 
1 1, 8k’™1,, = ‘ ~— 


where Y, and Y, are the functions defined in eq 23, which may be 
taken from tables. 


Y,(k) =2.086, 9.100, 8.165, and 6.072 for ky, k., ks, and k, 
Y,(k)=1.490, 25.663, 20.049, and 10.796 
¥(k)=0.401, 9.04, 6.57, and 3.41 


C= = —a2, Y (ky) +22 Y (ko) —23Y (k3) +24 ¥ (hy) } =0.86 dyne. 
Wyre 
Hence 


C3+-Cz=0.86+0.10=0.96 dyne, which differs from the observed 
value 1 dyne by less than the experimental error. 
The force /,=94,527.22 dynes, so that 


F=94,527.22+0.96 cos 2a 


For azimuth a=0 or z, this azimuthal term amounts to 10 parts 


in a million. 
VIII. SUMMARY 


The self inductance between two coaxial helices and the force 
between them are given in terms of four w-functions, and their 
derivatives, respectively, as in eq 13 and 13’, the four z-distances 
being defined in eq 11 with their respective azimuthal differences ?.. 
Resolving these functions into three components, w(x) =ws(r) +e-(:) 
+-w,(2,0,), the principal part corresponds to the mutual inductance 
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between the two associated current sheets, with angular components 
of current, and w,(z) with axial components. The last term, w,(z,6,), 
is like w,(z), @ Sec ond- order correction term applied to the finite term, 
w(t), Where the ratio of pitch to cylindrical radius is considered a 
first-order infinitesimal. 

The principal term, we(x), is given by eq 37 and w,(x) by eq 39. 
The first 1s equivalent to Jones’ formula for current sheets. The 
sw, ir) and w,(2) are given by eq 42 and 48. 

iwinnuthal terms, w,(2,0,) and w,(z,6,), are given by the series 

and 49’ and in finite terms by eq 57 and 57’. 

: liese p Somes are not restricted to the case where the helices have 
all ernl number of turns but are quite general. They need not 
have the same pitch. For the case most common in practice where 
the numbers of turns are integers and the lead wires are straight lines 
parallel and close to the generator of the cylindrical current sheet, 
the total effect of all axial components of current in helices and Jead 
wires together is represented by a function @,(7) which is given by 
ey 46, and its a-derivative by eq 46’. 

is shown in section VII that the azimuthal variation of the force 
computed by these equations agrees (within experimental errors) with 
that observed at the National Physical Laboratory. 

The correction terms given in section VI take account of the finite 
diameter of helical wires. They require a knowledge of the current 
distribution in the wires. 

These formulas are valid to the second order, inclusive, provided 
no parts of the two helices lie very close to each other, that is, at 
distances of the order of magnitude of the pitch of the windings. 
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WasHinaton, December 21, 1938. 





